Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,4380TG31.8% 98.4 / 41.4 63pseudogene (1047/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base G (13/7);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAA  >  NC_000913/1466149‑1466709
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                
aGGACTGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                                                                                                                                                                                                                                                                         <  2:304289/299‑1 (MQ=255)
             ggATGTCAGTGGCGGCGGTCACGGTATTGATGTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGGGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                                                                                                                                                                                                            <  2:178644/299‑1 (MQ=255)
                             ggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                                                                                                                                                                                                                                           >  1:626316/1‑300 (MQ=255)
                                   ggTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGTCGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCAGCACCAACTGTGGCACCGAAACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt                                                                                                                                                                                                                                        <  1:555341/297‑1 (MQ=255)
                                      aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGCCCGACGGCGACAAAGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCAAACAGATTAATGGTGATGAAG‑CAAAGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGAAAGGATTCCAACGGTAAGGAAATCGATGGTAATAACGGGAAAGTGATTCAGGACGGCGACCTGGATGTCAGCGGCGGCGGTCGCAGTATTgct                                                                                                                                                                                                                                  >  2:352211/1‑298 (MQ=14)
                                      aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGCAATAACGGAAAAGTGATTCAGGATGGCGATCTGGGTGTCGGCGGCGGCGGTAACGGTATTGat                                                                                                                                                                                                                                  >  2:140863/1‑300 (MQ=255)
                                       ttGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGTTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                                                                                                                                                 <  1:686347/300‑1 (MQ=255)
                                       ttGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                                                                                                                                                 <  2:626316/300‑1 (MQ=255)
                                       ttGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                                                                                                                                                 <  1:685931/300‑1 (MQ=255)
                                                    tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGATGGCGATCTGGCTGTCAGCGGCGGTGGTCACGGTATTGattccaaatacaaca                                                                                                                                                                                                                     >  1:372476/1‑286 (MQ=11)
                                                             gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGTTGATGACG‑CCACGGCAAACAACAACTGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACgg                                                                                                                                                                                                            <  1:166729/299‑1 (MQ=255)
                                                              aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCCAATCGACGGTGACAAGACCATCGTCAATAATGAAGGTGAGAGCACCATCACTAACGGTGGCACCGACACGCAAGTTAACGGTGATGACG‑CCGCGGCAAACAACAACGGCAAA‑ACCA‑CTGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAGGGTGATTCAGGAAGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa                                                                                                                                                                                                                              >  2:101811/1‑280 (MQ=11)
                                                              aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGCGCCAAGGCCATCGTCAATAACGACGGCGAGAGCGCCATCACCACCGGTGGCACCTGCACACCGATTACCGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGCAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCCCAGGCGACAGCGAAACGGt                                                                                                                                                                                                           >  2:100615/1‑299 (MQ=14)
                                                              aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCAACAACGGTGGCACCGGCACGCAGATTAACGGTGATGACG‑ACAAGGCGAACAACAACGGCACA‑AACA‑CCGTTGACGGCAAGGATTCCAACGGTACGGAAATTGCCGGTAATAACGGGGAGGTGATTAAGGACGGCGATCTGGGTGTCAGCGGCGGCGGTCACGGTATTGATATCCACGGCGACAGCGCAACGGt                                                                                                                                                                                                           >  2:100507/1‑299 (MQ=11)
                                                                 ggggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                                                                                                                                                       <  1:360316/298‑1 (MQ=255)
                                                                 gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCGCCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                                                                                                                                                       <  1:360312/300‑1 (MQ=255)
                                                                 gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGGTTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                                                                                                                                                       <  1:360322/300‑1 (MQ=255)
                                                                         cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAATAACAACGGCAACCACC‑‑CCGTCTACACCCCTAATCACCCCTGTGCGAGCTTCGACTATAGTGGCGGAAAAGTTATTCAGGACGGCGATATGGATGTAGGCGGCGGAGGTCACGGTATTGATATTAACGGCGACAGCGCAACGGTGGATaataa                                                                                                                                                                                                  >  1:189804/1‑294 (MQ=12)
                                                                         cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCNGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATAACaa                                                                                                                                                                                                  >  1:516217/1‑297 (MQ=255)
                                                                                 ccATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                                                                                                       >  1:125244/1‑300 (MQ=255)
                                                                                               gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCNGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATc                                                                                                                                                                          <  2:516217/298‑1 (MQ=255)
                                                                                               gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAAAGGTGGATAAAAAGGGCACCATGACCGTCAA‑CGATcc                                                                                                                                                                         >  1:21375/1‑299 (MQ=255)
                                                                                                       gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGCTATCCCTCCCACCCACGTCGACACCGACAATCAGATTAATATTACTACCG‑CCCCAACCCACCACAACTCCAAA‑CCCC‑CCATTTACGGCAAGGATTCAAATGGTACGGAAATCGACGGTAATAACGGAAAAGTGAGTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCGACGGCGACAGCGCAACGGTGGATAACAAGGGTACCATGACCGT‑ACACGATCCGGAGTCCa                                                                                                                                                                 >  2:356890/1‑299 (MQ=2)
                                                                                                       gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGCTAGCGCCCCCACCCACGTCGAAAAAGACAACCAGATGAAGAGTAAT‑ACGACCCCAACAAACCACAAATGCAAACACC‑‑CCACTGACAGCAAAAACTCCAACGGTAAGGAAATCAACGGTAATAACGGAAAAGTTATTAAGGAAGGAGATCTGGAGGTCAGCGGCGGCGGTAAAGGTATTGACATCAAAGGCGAAAGCGCAACGGTGGATAACAAGGGCAccctgcacagcaccaatccagagtccgt                                                                                                                                                                 >  2:356853/1‑274 (MQ=255)
                                                                                                       gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGTAATAACGGGAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                                                                                                  >  2:653831/1‑267 (MQ=255)
                                                                                                       gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                                                                                                                                                               <  1:416593/238‑1 (MQ=14)
                                                                                                       gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                                                                                                                                                               >  2:416593/1‑238 (MQ=14)
                                                                                                         ccATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTCATGATACTGACG‑CAACGGCAAACAACCACGG‑AAA‑ACCAACCGTTGCATACAAGTATTCCAAAGGTAGGGAAATAAAAGGTAATAACGGAAAAGTGGTTCAGGATGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTCTTGATATCAACGGCGACAGCGCAACGGTGGATAACAAGGGCACCCTGACTGTCGC‑CGATCCGGAGTCCATc                                                                                                                                                               >  2:397745/1‑299 (MQ=17)
                                                                                                           aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCgg                                                                                                                                                             <  1:101811/299‑1 (MQ=11)
                                                                                                            tGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGGAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTa                                                                                                                                                           >  2:692605/1‑300 (MQ=11)
                                                                                                              ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                                                                                                                                                       <  2:699776/207‑1 (MQ=11)
                                                                                                              ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                                                                                                                                                       >  1:699776/1‑207 (MQ=11)
                                                                                                               gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAccctggcaccatacacgtccccaat                                                                                                                                                                            >  2:267349/1‑257 (MQ=255)
                                                                                                               gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATAAc                                                                                                                                                                                                    >  2:622324/1‑257 (MQ=11)
                                                                                                               gTATCCAGATCGACGGTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                                                                                                                                                    <  1:622324/257‑1 (MQ=11)
                                                                                                                     aGATCGACGTTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACTGTGGCACCGGCACCTAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTg                                                                                                                                                  <  1:377496/300‑1 (MQ=255)
                                                                                                                       aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATNGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGAc                                                                                                                                                <  1:137619/300‑1 (MQ=255)
                                                                                                                                           cGTCAATAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAAAGGTGGATAACAAGGGCACCATGAACGTCAC‑CGATCCGGAGTCCATCGGTATCCAGATTGAAGGCGAACAGGccattattaa                                                                                                                            >  2:117345/1‑291 (MQ=11)
                                                                                                                                                                      cACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCANCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAATCGTCAACAACGAAGGCGAGAGCACCATCACCaa                                                                                                 >  2:571579/1‑300 (MQ=255)
                                                                                                                                                                       aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGTTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAc                                                                                                <  2:125244/300‑1 (MQ=255)
                                                                                                                                                                            cGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGt                                                                                             >  1:583644/1‑298 (MQ=255)
                                                                                                                                                                            cGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCAGGAg                                                                                                                                                                     <  1:185181/226‑1 (MQ=255)
                                                                                                                                                                            cGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCAGGAg                                                                                                                                                                     >  2:185181/1‑226 (MQ=255)
                                                                                                                                                                                         atttggATTAATGGTGATAATG‑CCACGGTAAACAACAATGGCAAA‑ACCA‑CCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGTCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaac                                                                                                                           <  2:319483/250‑1 (MQ=14)
                                                                                                                                                                                             ggATTAAGGGTGATGATG‑CCACGGTAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTc                                                                                                                                                                                                                                                          <  2:191931/124‑1 (MQ=11)
                                                                                                                                                                                              gATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAgtttct                                                                                                                              >  2:687677/1‑244 (MQ=255)
                                                                                                                                                                                              gATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCAACGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACC‑‑GCGGTTGTTCACACAGAAAGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATCAAc                                                                        >  1:258114/1‑299 (MQ=18)
                                                                                                                                                                                              gATTAAGGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                                                                                                                                                            <  2:429440/154‑1 (MQ=255)
                                                                                                                                                                                               aTTAATGGTGATGACG‑CCACGGTAAACAACAATGGCAAA‑ACCA‑CCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                                                                                                                                                 <  1:685844/148‑1 (MQ=255)
                                                                                                                                                                                               aTTAAGGGTGATGATG‑CCACAGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCgg                                                                                                                                                             <  2:117399/215‑1 (MQ=255)
                                                                                                                                                                                                     ggTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCAATCGTCAACAACGAAGGCGAGAgcacacc                                                                                                       >  1:563874/1‑258 (MQ=255)
                                                                                                                                                                                                     ggTGATAATG‑CCACAGTAAACAACAATGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACgg                                                                                              >  1:697028/1‑272 (MQ=255)
                                                                                                                                                                                                     ggTGATAATG‑CCACAGTAAACAACAATGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACgg                                                                                              <  2:697028/272‑1 (MQ=255)
                                                                                                                                                                                                                     gtaaaCAACAATGGCAAA‑ACCA‑CCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATAGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac                                                   >  1:196345/3‑300 (MQ=255)
                                                                                                                                                                                                                     gCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAA‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAACAACGGTGGCACCGGCa                                                                                   >  2:178286/1‑268 (MQ=255)
                                                                                                                                                                                                                                 ggCAAA‑ACCA‑CCGTTGACGGTAAAGATTCCAACGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGAACAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAACAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaa                                       >  1:289282/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        aaagATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGTTCACGGTATTGATATCACCTGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCTGTTGGTAACAACGAAGGCGAGAGCGCCATCACAAAAGGGGACAACGGCACGCGGATGAACGGTGATGATGTCACGGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                   <  2:196345/296‑1 (MQ=17)
                                                                                                                                                                                                                                                           gATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACCATCACCAATGGTGGCACCGGCACGCatgttca                                                                         >  2:136209/1‑236 (MQ=255)
                                                                                                                                                                                                                                                           gATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCGATTGTTAACAACGAAGGCGAGAGCACCATCACCAATGGCGGCACCGGCACTCAGATTAACGGTGATGACGCCACAGCAAACAACAACGGCaaaaaca                                   >  2:275571/1‑277 (MQ=14)
                                                                                                                                                                                                                                                                          gAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAgg                                                                                                                   >  1:636286/1‑185 (MQ=255)
                                                                                                                                                                                                                                                                          gAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAgg                                                                                                                   <  2:636286/185‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCCCCGTTGACGGCAAGGATTCCACCGGTACGGaaa  >  2:244481/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                            aaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa  <  1:244481/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                
AGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACG‑CCACGGCAAACAACAACGGCAAA‑ACCA‑CCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAC‑CGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAA  >  NC_000913/1466149‑1466709

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑