Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,287,040 G→A intergenic (+56/+26) ychS → / ← rttR putative uncharacterized protein YchS/small RNA RttR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0400GA93.6% 123.4 / ‑4.1 49intergenic (+56/+26)ychS/rttRputative uncharacterized protein YchS/small RNA RttR
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (21/23);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAG  >  NC_000913/1286745‑1287330
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                   
ttCGCTCTTCGGGGCGTTGTCTGAGGCAACGCTCGCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCGCACCCTTCCCCGCTTGGTCAGAATACTTAATTGCTGATTCGTTTTAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAagag                                                                                                                                                                                                                                                                                                <  1:469642/300‑1 (MQ=17)
ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTGGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAagag                                                                                                                                                                                                                                                                                                <  1:219174/300‑1 (MQ=255)
ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATGTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAagag                                                                                                                                                                                                                                                                                                <  1:221366/300‑1 (MQ=255)
ttCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTTGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAGTATTTGATTGCGGATTCGTTTTAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAagag                                                                                                                                                                                                                                                                                                <  1:469815/300‑1 (MQ=255)
                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAg                                                                                                                                                                                                                                                                              <  1:217026/300‑1 (MQ=255)
                           aaCGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGAtggtg                                                                                                                                                                                                                                                                      <  1:233039/299‑1 (MQ=255)
                                tctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGACGGTGGTgggg                                                                                                                                                                                                                                                                <  1:633613/300‑1 (MQ=255)
                                 ctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTAGGGTCGTTGCCTGCGACCCCGCTCTCTCGCTGACACCCGAGTCGAACATTAGTCGAAGCTTCTCATCCTTCCACGACTGGGTAGGATATTTACTTGCGGATTCGCTTAGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                                                                                                                                                                                                    >  2:727147/1‑295 (MQ=14)
                                                 aGTCGAATCTTGGTCGAAGTTTCTCATCCTTCCCCGCTTGTGCAGAATATTTGATTGCTGATTCGTTTGAGAATTCCGTGGCTTTTGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                               <  2:59073/300‑1 (MQ=255)
                                                      acacTTGGTCGAAGCTTCTCATCCTTCCCCGTGGGGGCAGAGTATGTGATTGCGGATTAGTGTTAGGGTTCCGTTGCTTGTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                                                                                                                                                                                                          <  2:676816/297‑1 (MQ=17)
                                                      aaCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTAGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCAGAACCTTAGTCGAAGCTTCTAATCCTTCCCCGAATGGGCAGAATATTTGATTGCGGAGTCGCTTAAGAGTTCAGGGAc                                                                                                                                                                                                                                                                                       >  2:50431/1‑255 (MQ=255)
                                                                    cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                                                                                                                                                                              >  1:469370/1‑298 (MQ=255)
                                                                             cTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCTGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCCTCGCTACGCGCTTCGCCCTTCGGGTCGtttcc                                                                                                                                                                                                                   >  1:665182/1‑297 (MQ=255)
                                                                               tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACCCAGCGCTGCGCGCCTCGCCCTTCGGGTCGTTGCCTGCGACCACGCTCTCTCGCTAACGCTCGAGCCGAACCTTAGTCGAAGCTTCTAATCCTTCACAGCGTGGGCAGAGTATTTGCTTGCGGATTCGCTTCATAATCCAGGAACTTTTGCACGTGGTGGTGGTGGGGCAAGACTTCCTCAGCGCTGCGCGCTTAGCACTTCGTGTCGTTGCCt                                                                                                                                                                                                                  >  2:306152/1‑299 (MQ=255)
                                                                                                            gATTCGTTGGAGAATTCCGGGACTTGTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTACGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                                                                                                                                                                                                                                                              <  2:64374/226‑1 (MQ=255)
                                                                                                                      agacttccggtacTTTTGAAAGCGATGGTGGTGGGGGAAGGATTACGCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc                                                                                                                                                                                                                                  <  2:45409/232‑1 (MQ=17)
                                                                                                                                      tGAAAGTGATGGTGGTGGGGGAAGGACTACTTAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCACCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                                                                          <  1:639225/300‑1 (MQ=255)
                                                                                                                                             gATGGGGGTGGGGGAAGGACTACTCAGCGCTGCGCGCTTCGCCTTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCGAAGTTTCTTATCCGTCCCCGTGTGGGCAGAATATGTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                               <  2:59067/208‑1 (MQ=11)
                                                                                                                                               tgggcgtggGTGAAGAANTACTCACCGCTGCGCGCTNCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATGTGGGATTGCGGATTCGTGTAAGGGTTCAGGGATTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTNCTCATCCTTCCCCGCa                                                                                                                                                 <  1:611955/295‑1 (MQ=14)
                                                                                                                                                ggtggtggtGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCgatg                                                                                                                                                <  1:360809/300‑4 (MQ=14)
                                                                                                                                                             ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                                                                   >  1:90166/1‑300 (MQ=255)
                                                                                                                                                              gATTACTCAACGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTATTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGCGTAAGAGTTCAGGGACTTTGGAAAGTGATGGTGGGGGGGGAAGGATGACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                                                                                <  2:103855/238‑1 (MQ=14)
                                                                                                                                                                  aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAANAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTAGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGTCGCAGCTTCTAATCCTTcc                                                                                                                                                      >  2:533288/1‑276 (MQ=255)
                                                                                                                                                                       gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAGTATTTGCTTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCgatg                                                                                                                                                >  2:688911/1‑274 (MQ=255)
                                                                                                                                                                       gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                                                                                         >  1:655172/1‑300 (MQ=255)
                                                                                                                                                                       gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCgatg                                                                                                                                                <  1:688911/277‑4 (MQ=255)
                                                                                                                                                                               cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                                                                                    >  1:94744/1‑297 (MQ=255)
                                                                                                                                                                                   tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                                                                             >  1:91929/1‑300 (MQ=255)
                                                                                                                                                                                        cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGTGGCAACGCTCTCTCGATGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCAGCATGGGAAAAACATTTAATTGCGGAtctgtgggggagttca                                                                                                          >  2:445319/1‑283 (MQ=17)
                                                                                                                                                                                          tCGGGTCTTTGCCTGCGGCAACGCTTTCTCGCTTACGCTCGAGTCGAGCCTTAGTCGAAGCTTCTCACCCTTCCCCGCGTGGGCAGAATGTGTGATTGCGGATTCGCTGAAGAGTTCAGGGACTTGGGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGa                                                                                                      <  2:295710/300‑1 (MQ=17)
                                                                                                                                                                                                   tGCCTGCGTCAACACTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCTTTCCCCGCATGGGCAGAATATTTGGGTGCGGATTCGCTTAAGATTTCATGGACTTTGGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                                                                             <  2:655172/300‑1 (MQ=17)
                                                                                                                                                                                                         cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTAGCCATTCGGGTCGTTGCCTGTGGCAACGCTCTCTCGATGGCGCTCGCGTCGAACCTTGGTCGAAGCTTCTCATCCTTCACCGCATGGGCAGa                                                                                                                                         >  2:355500/1‑250 (MQ=17)
                                                                                                                                                                                                         cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTAATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCCTAAGAGTTCCGGGACTTTTGAAGGTGATGGTGGTGGGGGAAGGATTACTCGCCGCTTCGCTCCTCGCCCTTCGGGCAGTTGCCTGTGGCAACGCTCTCTCGCTTTCGCTCGAATCGAAcatt                                                                                                                                                                            >  2:409589/1‑212 (MQ=255)
                                                                                                                                                                                                                 ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                  >  1:257568/1‑297 (MQ=255)
                                                                                                                                                                                                                                tCGAGTCGAACCTGGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATGCGATTGCGGATTCGCTTAAGAGTNCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc                                                                 <  1:533551/299‑1 (MQ=255)
                                                                                                                                                                                                                                       gAACCTTGGGCGAAACTTCTCACCCTTCCCCGCATGGGCAGAATACTTGATTGCGGATTCGCTTAATAGTTCGTGGACGTGTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAg                                                         <  2:108881/300‑1 (MQ=17)
                                                                                                                                                                                                                                         aCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGCACCTTAGAAGTc                                                       >  1:67197/1‑300 (MQ=255)
                                                                                                                                                                                                                                         aCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGAAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGCAGTc                                                       >  1:67073/1‑300 (MQ=255)
                                                                                                                                                                                                                                         aCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGAAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTc                                                       >  1:66855/1‑300 (MQ=255)
                                                                                                                                                                                                                                              ggTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTCCGAAGTCGATGa                                                  >  2:399636/2‑300 (MQ=255)
                                                                                                                                                                                                                                                            aTCCTTCCCCGCATGGGCAGAAGACTGAATTACGGATGAGCTGGAGAGTGCAGGTATTTGTGAAAGTGAGGGTGTTGGGGGAAGGATTACTCAGCGCTGCTCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAg                                    <  2:205529/300‑1 (MQ=18)
                                                                                                                                                                                                                                                                ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCt                                 >  1:227925/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTCTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCTCTTAGCACTCCGGGTCGTTGCCTGCGGCACCGCTCTCTCGCTGACGCCCGAGCAGAACCTTGATCGAAGCTTCTCATCCTTCCCCGCATTTGCAGAATCTTCAATTGCAGATTTGTTGGAAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGAAACTTCGAAGTCGATGACGGCAGATTTACAGTCt                                 >  1:228975/1‑299 (MQ=17)
                                                                                                                                                                                                                                                                           ggCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGc                     <  1:727147/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ggCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGc                     <  1:727138/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCt                 >  1:712372/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                 atatTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTc                >  2:51337/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                              gATTCGCTTGAGAATTCAGGGACTTTTGAAAGTGATGGTGGTTGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAg  <  1:6947/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   gCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAATCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGAATGGTTAGAATCTTTa                                                                                                                                  >  1:570446/1‑167 (MQ=14)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                   
TTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAG  >  NC_000913/1286745‑1287330

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: