Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 559,176 T→A 17.8% V68E (GTA→GAA)  sfmC → putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913559,1760TA17.8% 118.4 / 6.0 52V68E (GTA→GAA) sfmCputative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base T (36/6);  new base A (2/8);  total (38/14)
Fisher's exact test for biased strand distribution p-value = 1.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAA  >  NC_000913/558945‑559414
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
ccGTCAAGAATTCACCAGGCTTTTTAACTATTATGACTAAAATAAAGTTAGTGATTCTCATTTTATTTTATTTAATCATTTCGGCCAGCGCCCAAGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑Taaaagaa                                                                                                                                                                                                                                         <  2:878277/239‑1 (MQ=255)
         aTACCCAAGGCGTTTTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGtcatt                                                                                                                                                                                                                                 <  2:642786/238‑1 (MQ=255)
              caAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATTCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTcatcat                                                                                                                                                                                                                           <  1:935376/239‑1 (MQ=255)
                      aTTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGtcatt                                                                                                                                                                                                                                 <  1:234665/225‑1 (MQ=255)
                                atgaCTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAACTAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGGG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACtgtttgt                                                                                                                                                                                                         >  1:510862/1‑239 (MQ=255)
                                atgaCTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATTAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACACAAGCGGGG‑TAAAAGAAAAGTCATTCATCATTACAcctgtctctaaaaa                                                                                                                                                                                                         >  2:234665/1‑227 (MQ=255)
                                     tAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCATGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTAAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TCAAAGAAAAGcaattca                                                                                                                                                                                                                               >  1:337849/1‑212 (MQ=255)
                                        aaTAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCCTGCGGTATTGATTAGACATAGCCATACCAATGAGCGCTTTCTGGCCAA‑TTCGTGGATTGAAAACAGCAGCGCTG‑TAAAAGAAAAGTCATTCCTCATTACCCCGCCACTGTTTGTTAGAGaa                                                                                                                                                                                                  >  2:218412/1‑238 (MQ=255)
                                             aGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTAATTAATCATTACACCGCCACTGTTTGTTAGGAAACCCaaa                                                                                                                                                                                            >  2:580641/1‑239 (MQ=255)
                                                                    tatttaATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtattat                                                                                                                                                                     >  1:512416/1‑239 (MQ=255)
                                                                         aaTCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAATCGAAAAAACTTTGCGAATTATTTac                                                                                                                                                                 >  1:1012316/1‑238 (MQ=255)
                                                                           tCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTAGACCGCCACTGTTTGTTAGTGAACACAAAAACGAAAAAAATTTACTTATAATTTACAcc                                                                                                                                                              >  2:462499/1‑239 (MQ=255)
                                                                           tCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCCCCACCGTTGTTTAGTGAACCCAAAAACAAAAAAACATTCGGAATTATTTAAAcc                                                                                                                                                              >  2:422718/1‑239 (MQ=255)
                                                                                tCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACAGCACCGTTTTGTTAGTGAACCCAAAAGCAAAAATACTTTGCGCATTATTTAAACAGGTcc                                                                                                                                                         >  2:396844/1‑239 (MQ=255)
                                                                                     cccccccctggctgcGGGAGGTGTCGTATGGGGTGCCGGGTTTATTATTTCCCCCCATGATGAAAACCATGCGGCGGTATGGATTAAAAATAGCCATCCAAAGGAGCCTTTTCTGGTAAA‑TTCGTGGATTGAAAACAGCAGCGGTA‑AAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGc                                                                                                                                                    <  2:541189/229‑1 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTGAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCCAGTTCTGGTCAA‑TTCGAGGATTGAAAACAGCAGCGGTG‑AAAAAGAAAAGTCATTCATCATTACACCTCAACTGTTTGTTAGTGAACCCAAAAGCGAAAAACCTTTGGGAATACTTTAAACAGGCCCACCCCGGgaagca                                                                                                                                            >  2:549988/1‑239 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TACAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCCCTGgcagca                                                                                                                                            >  1:685708/1‑239 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATCATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGGGAAAATACTTTGCGaattat                                                                                                                                                                     >  1:723445/1‑214 (MQ=255)
                                                                                                    cGGAGGCACCGCATGAGGGACCACGCTTATTATTTACCCCGCTGATGCAAAACAGACTGCGGTAGGGATTAGAAATAGCAATACAAATGAGCGCTTTCTGGACAA‑TTCGGGGATTGAAAACAGCAGCGGTT‑AAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTg                                                                                                                                     <  1:218412/239‑1 (MQ=255)
                                                                                                           aTCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCACCGGGG‑ACAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCAAAAAGCGAAAATACTTTTGGTATTTTTTAAACCGGTCCCCCGCTGgcagcaactccaggatc                                                                                                                                 >  2:946343/1‑225 (MQ=255)
                                                                                                             cGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTAAACCGGTCCCCCCCTGgcaagaaaacggagggtcttgt                                                                                                                            >  2:751638/1‑220 (MQ=255)
                                                                                                                    ggTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATCCTTGGCGTATTATTTAAACCGGTCCACCCCTGGCAGCAGATCGGGAGTTTTTGTTTTGGAt                                                                                                                     >  2:581827/1‑239 (MQ=255)
                                                                                                                     gTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTTGGTATTATTTACACCGGTCCACCGCGGGCGGCCAACCGTGAGTCTCTGTTTTGGGTg                                                                                                                    >  2:23258/1‑239 (MQ=255)
                                                                                                                           ggggttttttttacCCCGTTGTTGCTAAACAGATTGCGGTATGGATTAAAAAAAGCCATACAAATGGGCGTTTTCTGTTCAA‑TTCGGGGATTGAAAACAGCAGCGGTG‑AAAAAGAAAAGTCATCAATATTTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGtt                                                                                                              <  1:422718/230‑1 (MQ=255)
                                                                                                                                attattTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAgtc                                                                                                         >  1:539917/1‑237 (MQ=255)
                                                                                                                                      tatCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTTAACCCAAAAGCGAAAATACTTTGCGACTTATTTACACCGGTCACCCGCTGGCGGAAGATCGTGGGTTCCTGTTCTGGATGAATGTTAAAAGGATc                                                                                                     >  2:663237/1‑237 (MQ=255)
                                                                                                                                       atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                                  >  1:721949/1‑239 (MQ=255)
                                                                                                                                       atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                                  >  2:310216/1‑239 (MQ=255)
                                                                                                                                         cccGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTc                                                                                                >  1:176396/1‑239 (MQ=255)
                                                                                                                                             cTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTa                                                                                            >  1:1059627/1‑239 (MQ=255)
                                                                                                                                              tGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAg                                                                                           >  1:288935/1‑239 (MQ=255)
                                                                                                                                                    aaaaCAGACGGGGGTATGGATTAAAAATAGCCATCCCAATGAGCGTTTTCTGGTCAA‑TGCGTGGATTAAAAACAGCAGCGTTG‑TAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATaaaa                                                                                     <  1:23258/238‑1 (MQ=255)
                                                                                                                                                          gACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGGATTATTTACACCGGTCCACCGCGTGCAGAAGATCGGGAGTCTATGTTCTGGATGAATGTTAAAACGATCCCTGCGGGAGATAAAAAAGCAt                                                                               >  2:474820/1‑239 (MQ=255)
                                                                                                                                                             tGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTAAACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTTTGGCTGAATGTTAAGACGACCCCTCCGGGAGATAAAAATGCAttaa                                                                            >  2:563696/1‑237 (MQ=255)
                                                                                                                                                                  tttgtataaaaaaTAGCCAACCAAATGGGTGTTTTTGTGTAAA‑TTCGGGAATTAAAACCAGAAGGGGTA‑AAAAAAAAAGTTCATTCACCATTACCCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                                        <  2:288935/228‑1 (MQ=255)
                                                                                                                                                                  tttGGATTAAAAAGAGCCAAACAATGGAGCGCTTTCT‑GTCAATTTCGTGGATTGAAAACGGCAGGGGTG‑AAAAAAAAAAGTAATACATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                                        <  1:751638/236‑1 (MQ=255)
                                                                                                                                                                       aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGaa                                                                  >  1:997208/1‑239 (MQ=255)
                                                                                                                                                                       aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCAACCGCTGGCAGCAGATCGTGAGTCTCTATTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGAATTGAAAGGCCGGaa                                                                  >  1:389613/1‑239 (MQ=255)
                                                                                                                                                                                aGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCa                                                         >  1:1015822/1‑239 (MQ=255)
                                                                                                                                                                                 gCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTGGATAAAAATGCATTGAACGACAGGAATGTTTTGCaa                                                        >  1:104855/1‑239 (MQ=255)
                                                                                                                                                                                 gCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCAACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGAACCCTTCGGTCGATAAAAATGCATTGAACGGCAGGGAATTTTTGCaa                                                        >  2:154567/1‑239 (MQ=255)
                                                                                                                                                                                           tGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGTTAGGTCAAAATGCATTGAACGGCAGGAATGTTTTGCAAAGGGCGAttt                                              >  2:657496/1‑239 (MQ=255)
                                                                                                                                                                                               cgcTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAAGGCATTTAACGGCAGGAATGTTTTGCAAATGGGGATTTTATc                                          >  2:395473/1‑239 (MQ=255)
                                                                                                                                                                                                  tttCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGAATTGAACGGCAGGAAATTTTTGCAACTGGGGTTTTTATcccc                                       >  2:575809/1‑236 (MQ=255)
                                                                                                                                                                                                   ttgttgttaattccggggATGGAAAACGGCAG‑GGTGAAAAAAGAAAAGTCATTCATCATTCCACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATcgcg                                       <  1:154567/223‑1 (MQ=255)
                                                                                                                                                                                                    tgtggaaaa‑TTTGTGCATGAAAAACATCAGGGGTA‑AAAAAAAAAAGTCATTCATCTTTCCTCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCa                                     <  1:563696/232‑1 (MQ=255)
                                                                                                                                                                                                        gttttttttgtgtATTAAAAACAGCGGGGGTG‑AAAAAAAAAATGTAGTCACCATCACCCCGCCCCTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGa                                  <  2:510862/230‑1 (MQ=255)
                                                                                                                                                                                                        gTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGaa                                 >  2:876431/1‑239 (MQ=255)
                                                                                                                                                                                                            a‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGCAAGGAATGTTTTGCAACTGGCGATTTTACCGCGCATGAAatta                             >  1:121291/1‑239 (MQ=255)
                                                                                                                                                                                                                   ggATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGCTCCCTTCGGTAGATAAAACTCCATTGAACGGCAGGAATGTTTTGCAACTGGGCATTTTAACGCGCATGAAATTAtttat                        >  1:439891/1‑236 (MQ=255)
                                                                                                                                                                                                                                   cGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAAc                                                                          >  2:353194/1‑172 (MQ=255)
                                                                                                                                                                                                                                   cGGTG‑TAAAAAAAAATTCATTCACCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAAc                                                                          <  1:353194/172‑1 (MQ=255)
                                                                                                                                                                                                                                         tAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGaa  >  1:665045/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
CCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTG‑TAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAA  >  NC_000913/558945‑559414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑