Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6260GA27.7% 124.7 / 16.1 66F180F (TTC→TTT) ycfTinner membrane protein YcfT
Reads supporting (aligned to +/- strand):  ref base G (37/10);  new base A (0/18);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 2.63e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGC  >  NC_000913/1174395‑1174866
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
ttGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTc                                                                                                                                                                                                                                >  2:248321/1‑250 (MQ=255)
  gATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTcagc                                                                                                                                                                                                                             >  2:481106/1‑251 (MQ=255)
   aTAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTcagc                                                                                                                                                                                                                             >  1:28848/1‑250 (MQ=255)
    tAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTcagc                                                                                                                                                                                                                             >  2:349939/1‑249 (MQ=255)
       aaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTac                                                                                                                                                                                                                        >  1:104391/1‑251 (MQ=255)
           gCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTACGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTacaaac                                                                                                                                                                                                                    >  2:142551/1‑251 (MQ=255)
                    ccaccaaAACCGATACCAGCGACAACAGCAGGGAGATAGAAAACACCAAGC‑CCCCGCCCGCAAAAAGGGAAAGAAAAGAA‑CCCATCAGCAGATGGCGG‑CGAAGC‑GGCACCTCTTAAACACAGGTCAT‑TTTTGTGGCGCAAACCCATGCGCCAA‑GGTTGAAAAAAGGCAAA‑TTGCGGAC‑CCCACTG‑TTCATCCCCCCCCACGGCGTGGGAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                                                                                                                                           <  1:195755/250‑1 (MQ=255)
                     caccacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑ACCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGAACTAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCAACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                                                                                                                                           >  2:504340/1‑250 (MQ=255)
                     caccacAATCGATAACAGCGACAACAGCAGCGAGATAGAAAAAAACCAGG‑CACAGACCGCCAGAACGGTAAACAAAAAA‑ACCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGACAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGGGGGAA‑CGTACTGAA‑CCGACACACTCAGCAGTACAAAACAGGCGAAt                                                                                                                                                                                                           >  2:504326/1‑250 (MQ=255)
                         acaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAgagg                                                                                                                                                                                                      >  1:91524/1‑248 (MQ=255)
                             ccGATACAAGCGAAAACAGCAGCGAGAGAGTAAACACCCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGC‑CGCAAC‑GGCACCTCTTTACCACGGGTCAG‑TATTGTGGCGCCAACCCAGACGCCAA‑GGCTATAATAAGGAAAA‑TTG‑GGATGCACACTC‑TTCCTCCCCCACCGCGGCGTGGGAA‑CAAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                                                                                                                                                  <  1:120787/250‑1 (MQ=255)
                                    cAGCGACAAAAGCAGCGAAATAGAAACCACCCAGG‑CCCCGACCGCCAGAACGTTAAGCAAAAAA‑GCCTCCAGCAGATGGCGG‑CGAACC‑GCCCCCTCTTTAACACAGTTCAT‑TATTTTTGCGCCAACCAATGCTCCAA‑GCCTGTAAAAAGGAAAA‑TGGCGAAT‑CACCCTG‑TTCATTCCCCACCGCGGGGTGGAAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                                                                                                                                           <  1:109655/251‑1 (MQ=255)
                                     aGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                                                                                                                                          >  1:53387/1‑251 (MQ=255)
                                       cGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCACACCACGGCGTGGGAA‑CGCAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGCGGCAGCGCCAGACGGc                                                                                                                                                                                        >  1:397776/1‑251 (MQ=255)
                                        gACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCAAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGAGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAACCAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                                                                                                                       >  1:416784/1‑251 (MQ=255)
                                        gACAACAGAAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGACAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTCAAAACACGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                                                                                                                       >  2:460659/1‑251 (MQ=255)
                                        aacaacaGCAGCGAGAAAGTAA‑CCACCAGGGCCCCGACCGCCAGAACGATAAGAAAAGAA‑GCCACCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTGGCCCCAACCCATGCGCCAA‑GGTGGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATCCCCCCCCACGGCGGGGGAA‑CAAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                                                                                                                       <  1:147797/250‑1 (MQ=255)
                                            acggaggcGAGATAGAAAACAACCAGG‑CCCCGCCCGCCAGAACGGTAATAAAAGAA‑GCCAGCAGCAGAGGGGGG‑GGCAAC‑GCCCCCTCTTAAACAAAGGTTAT‑TTTTGTGGCGCCAACCCATGCGCAAA‑GGCTGAAAAAAGGCAAT‑TTGCGGAT‑CCCTCTG‑TTAATCCCCCCCCACGGGGGGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                                                                                                                                   <  2:494762/245‑1 (MQ=255)
                                              agcagcGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                                                                                                                                                 >  2:343/1‑251 (MQ=255)
                                                     agaTAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                                                                                                                          >  1:266645/1‑251 (MQ=255)
                                                     aaaaagaaaaCAACCAGG‑CACCGCCCCCAAGAACG‑AAAACAAAAAACGCCATAACCAGA‑GGCGGGCGAACC‑GGCACCTCTTAACCACAGGTCTT‑TATTGTTGCGCAAACCCATGCGCAAA‑GGCTGTAA‑AAGGAAAAATGGCGAAT‑CACCCTG‑TCCATCCCCCCCCACGGCGGGGGAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                                                                                                                          <  2:266645/244‑1 (MQ=255)
                                                      aatagaaaccacccAGG‑CCCCGCCCGCCAGAGCGGAAAGAAAAGAA‑TCCATCAGCAGGTGGCGG‑CGCAAC‑GGCCCCTTTTTAACAAAGTTCAT‑TTTTTTGGCGCCAACCCATGCGCCAG‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CCCCCTT‑TTCTTCCCCCCCCGCGGCGTGGAAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                                                                                                                                         <  1:11422/239‑1 (MQ=255)
                                                        aaatAAAAAACCAGG‑CACCGACCGCAAGAACGGTAAGCAAAGAA‑GCCATCGGCAGATGGCGG‑CGCACC‑GGCCCCTCTTAAACAAGGGTCTT‑TATTTTTGCGCCAACCGACGCGCAAA‑GGCTGAAATAAGGAAAA‑TTGCGGAT‑CACACTG‑TTCTTTCCCCCCCGCGGGGTGGGAA‑CAAAATTAT‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                                                                                                                                       <  1:187588/248‑1 (MQ=255)
                                                           ataaccaacaggcccccgaccgCCAGAAGGAAAAGCAAAAAC‑GCCAAAAGCAAAGGGGGG‑CGCAAC‑GGAACTCCTTAAACACTGGTTTT‑TATTTTTGCGCCAACCCCTCCGCAAG‑GGCTAAAAAAAGGCAAA‑TGGCGGAT‑CACCCTT‑TTCAT‑CCCCCCCCCCGCGGGGGAAACGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                                                                                                                                                     <  1:232966/236‑1 (MQ=255)
                                                               aaaacgaggcccctgccgCCAGAACGGGAAGCAAAGAA‑GCCCACAGCAGGGGGCG‑‑CGGAACCGGCCCCTCTT‑AACCCAAGTCCTATCTTGTAGCCCAAAACCATGCGCTAA‑GGCTATAAGAAGGCAAA‑TTGCGTCC‑CCCACTC‑ATCATTCCCCCCCACGCTGTGGGAA‑C‑AAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                                                                                                                                 <  2:70318/236‑1 (MQ=255)
                                                                    cggcccccaacccccAGACCGGTAAGAAAAGAA‑GCCAACAGCAGATGGCGG‑CGCAAC‑GCCCCCTCTTTACCACAGGTAAT‑TATTGTTGCGCCAAACCATGCGCCAG‑GGCTGTAATAAGGCAAA‑GTGCGGAT‑CACACTG‑TTCATTCCCCCCCCCGGCGTGGGAA‑CGAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGAc                                                                                                                                                            <  1:241572/238‑1 (MQ=255)
                                                                    aGG‑CACCGACCGCCAGAACGGTAAGCAAAGCA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGAc                                                                                                                                                            >  1:162206/1‑250 (MQ=255)
                                                                         accgaccgCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                                                                                                                                        >  2:261398/1‑250 (MQ=255)
                                                                          ccgacccCCAAAGCGGGAAGCAAAAAA‑GCCATCAGCAGGGGGGGG‑CAAACC‑GGCCCCTTTTAAACACAGGTAAT‑TTTTGTGGCGCAAACCCAGGCCCAAA‑GGTTGTAAAAAGGAAAA‑TGGCGTAC‑CCCCCTT‑TTCTTCCCCCCCCAGGGCGTGGGAA‑GAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                                                                                                                                      <  2:436393/251‑1 (MQ=255)
                                                                           cgaccgCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAAc                                                                                                                                                     >  2:250775/1‑251 (MQ=255)
                                                                            gaccgCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACg                                                                                                                                                    >  1:402092/1‑251 (MQ=255)
                                                                                        ggTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                                                                                                                                        >  2:442969/1‑251 (MQ=255)
                                                                                         gTAAGAAAAGAA‑GCCATCAGCAGAGGGCGG‑CGCAAC‑GGCACTTTGTTAACACAGTACAT‑TATTGTTGCGCCAACCCTGGCGCCAA‑GGCTTTAATAAGGCAAA‑TTGGGGAC‑CACACTG‑TTCATCCCCCCCCAGGGCGGGGGAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                                                                                                                                       <  2:61866/251‑1 (MQ=255)
                                                                                            aGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                                                                                                                                     >  1:20750/1‑250 (MQ=255)
                                                                                                aagaa‑gCAATCGGCAGAGGGCGC‑CGAACC‑GGCACCTCTTTAACATAGGTATT‑TATTGTGGCGCAAACCCAGGCGCCAA‑TGCTGTAATAGGGCAAA‑TGGCGGA‑‑CCCCCTGGTCCATCCCCCCCCCCGGCGGGGAAA‑CGAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCATACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg                                                                                                                                <  1:149072/251‑1 (MQ=255)
                                                                                                                   gccgg‑caCAAC‑GGCACCTCTTAAACACAGCTCTT‑TTTTGTTGCGCCAACCCATCCGCAAA‑GGGTGAAAAAAGGCAAT‑TTGCGTAT‑CACCCTT‑TTCATACCCCCCCACGGCGGGGGAA‑CGAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatcatc                                                                                                              <  2:104391/249‑1 (MQ=255)
                                                                                                                    gggg‑cgCAAC‑GGCCCCTTTTTAACACAGGTCAT‑TATTGTTGCGCAAAACCAGGCGCCAA‑GGCTGTAATAGGGAAAA‑TTGCGGAT‑CACACTG‑TTCTTCCCCCCCCCCGGCGTGGGAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATcc                                                                                                             <  1:435128/249‑1 (MQ=255)
                                                                                                                          gaaaC‑GGCCCCTTTTAAACACAGTTCAT‑TTTTGTTGCCCAAAACCATGCCCCAA‑GGCTGTAATAAGGAAAA‑GTGGGGAT‑CACCCTG‑TTCATCCCCCCCCCCGGCGTGGAAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGc                                                                                                        <  1:395781/249‑1 (MQ=255)
                                                                                                                           cAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCGGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGttc                                                                                                        >  1:353817/1‑248 (MQ=255)
                                                                                                                                   aCCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGAGATCAACCCCAGCAGGAAcac                                                                                                >  2:474589/1‑249 (MQ=255)
                                                                                                                                       cttaaaaacaGGTCAT‑TATTGTGGCGCAAAACCATGCCCCAA‑GGCTGTAAAAAGGAAAA‑TGGCGGAT‑CACCCTG‑TTCATCCCCCCCCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAc                                                                                              <  1:481106/242‑1 (MQ=255)
                                                                                                                                                    aat‑tattGTTGCGC‑AAACCCAGGGCCAAAGGCTTAAATAAGGAAAA‑TTGCGGAT‑CACACTT‑TTCATCCCCCCCCAGGGCGGGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCa                                                                                <  1:168952/249‑1 (MQ=255)
                                                                                                                                                                     aaCCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGTACATCCCGCGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCCTTGGAGGc                                                               >  2:110660/1‑251 (MQ=255)
                                                                                                                                                                     aaCCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATAATCCCGGGCAGGAACTCACCGGGGGGAACGGGATAAGCGGCATTgggggg                                                               >  2:169304/1‑247 (MQ=255)
                                                                                                                                                                        cATGCGCCAG‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTAATTCCCCCCCACGGCGGGGGAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCAtt                                                            <  1:343/251‑1 (MQ=255)
                                                                                                                                                                         gggccccAA‑GGCTGAAATAAGGCAAA‑TGGCGGAC‑CACACTT‑TTCATCCCCCCCCACGGCGGGGGAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTg                                                           <  2:20750/246‑1 (MQ=255)
                                                                                                                                                                            ccccAA‑GGCTGTAAAAAGCCAAT‑TGGCGGAT‑CACACTG‑TCCATCCCCCCCCACGGGGGGGGAA‑CAAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCtt                                                        <  1:136534/249‑1 (MQ=255)
                                                                                                                                                                              ccAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTa                                                       >  1:81460/1‑250 (MQ=255)
                                                                                                                                                                                         aataaGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGGAACGGCATAAGCGGCATTGGAGGCAATGCTTAAATCGCGCTc                                             >  1:39518/1‑250 (MQ=255)
                                                                                                                                                                                          ataagcaaaa‑TTGCGGAT‑CCCACTA‑TCCTGTCCCCACCACGGCGTGGAAA‑CAAAATAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAg                                           <  1:220239/244‑1 (MQ=255)
                                                                                                                                                                                          ataaGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGGGCAGGAACTCACCGGTGGAATCGGCATAAACGGCATTggggggaatggttaaaacccgctcac                                           >  2:225992/1‑227 (MQ=255)
                                                                                                                                                                                           acagtaaaa‑TTGTGGAT‑CAC‑CTGGTTCATCCCCCACCACGGCATGGTAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAgg                                          <  1:221072/245‑1 (MQ=255)
                                                                                                                                                                                             aGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGGGCAGGAACTCACCCGTGGAATCGGCATAAGCGGCATTGGAGGCATTTGTTAAATCGCGCTCAGGTg                                        >  1:261621/1‑251 (MQ=255)
                                                                                                                                                                                                gaaatgggCGGAT‑CACCCTG‑TTCATCCCCCCCCGCGGCGGGGAAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgcc                                      <  1:169304/244‑1 (MQ=255)
                                                                                                                                                                                                 aaa‑TGGCGGAT‑CACACTG‑TTCATCCCCCACCACGGGGGGGAAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccag                                    <  1:214845/251‑1 (MQ=255)
                                                                                                                                                                                                       ggcggaacccacTG‑TTCATTCCCCCCCACGGCGTGGAAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccagcca                                 <  1:248321/241‑1 (MQ=255)
                                                                                                                                                                                                         ggAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAACGGCATTGGAGGCCTTGCTTAAATCCCGCTCAGGAGCCAGCCCCTgg                             >  1:378940/1‑251 (MQ=255)
                                                                                                                                                                                                          gAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAAGTGCCAGCCACTggg                            >  1:423843/1‑250 (MQ=255)
                                                                                                                                                                                                                 acTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCCACCCGCGCCGGGACACACCGGGGGAAACGGCATAAACGGCAATGGAGGCATTGCTTAAATCGCGCTCAGGGGCCAGCCCCTgggttagcg                      >  2:426087/1‑244 (MQ=255)
                                                                                                                                                                                                                 acTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCCTAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGGGCCAGCCACTGGTTCAgcg                      >  1:284899/1‑251 (MQ=255)
                                                                                                                                                                                                                   tG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCCAGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGGGCCAGCCCCTAGTTCAgcgc                     >  1:73337/1‑250 (MQ=255)
                                                                                                                                                                                                                     tttCATCCCCCACCCGGGGGTGGAAA‑CAAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCa                    <  1:349939/249‑1 (MQ=255)
                                                                                                                                                                                                                        cATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACCGGGTCAGCGCCCtt                 >  1:203136/1‑250 (MQ=255)
                                                                                                                                                                                                                          ttCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCCCAGGTGCCCGCCACTGGGTCAGCGCACTTAAc              >  1:453114/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cacGGCGTGGGAA‑CGAAATTAG‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACgccagcca      <  1:419939/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAAATCCCCGGTGGAATCGGCATAAGAGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCCGCGCACTTAACGCCAGCCACTGc  >  1:327078/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
TTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGG‑CACCGACCGCCAGAACGGTAAGCAAAGAA‑GCCATCAGCAGATGGCGG‑CGCAAC‑GGCACCTCTTTAACACAGGTCAT‑TATTGTTGCGCCAAACCATGCGCCAA‑GGCTGTAATAAGGCAAA‑TTGCGGAT‑CACACTG‑TTCATTCCCCACCACGGCGTGGGAA‑CGAAATTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGC  >  NC_000913/1174395‑1174866

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