Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	25.9%	intergenic (+112/+261)	yohP → / ← dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,229,177	0	A	C	25.9%	133.3 / 16.4	90	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (13/45);  new base C (17/2);  total (37/53)
Fisher's exact test for biased strand distribution p-value = 2.63e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326                                                                                                                                                        |                                                                                                                                                        tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑a                                                                                                                                                         <  1:1353786/150‑1 (MQ=255)   cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑gg                                                                                                                                                      <  1:1090965/150‑1 (MQ=255)    tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑ggg                                                                                                                                                     >  1:332130/1‑150 (MQ=255)     ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑gggg                                                                                                                                                    <  1:1346347/150‑1 (MQ=255)      gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATAGAATCACTGGCGCAAGGAGCGC‑GC‑C‑cgggg                                                                                                                                                   >  2:4083011/1‑150 (MQ=255)       tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑C‑gggggg                                                                                                                                                  >  1:2056292/1‑149 (MQ=255)        ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGggcg                                                                                                                                                 >  1:1359933/1‑150 (MQ=255)        ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGggcg                                                                                                                                                 >  2:235697/1‑150 (MQ=255)        ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATCATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGggcg                                                                                                                                                 >  2:235800/1‑150 (MQ=255)         gttaCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGggcgg                                                                                                                                                <  2:2909629/147‑1 (MQ=255)         gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGggcgg                                                                                                                                                <  1:3825341/150‑1 (MQ=255)         gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGggcgg                                                                                                                                                <  1:3026418/150‑1 (MQ=255)          tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑GC‑C‑ggggggggg                                                                                                                                               >  1:3209060/1‑149 (MQ=255)          tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑CGGggcggc                                                                                                                                               >  1:3209212/1‑150 (MQ=255)           gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Gg‑g‑gggggggggc                                                                                                                                              >  2:1197719/1‑148 (MQ=255)           gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑gggggggggc                                                                                                                                              >  2:114571/1‑148 (MQ=255)               ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGTAGCGC‑GC‑A‑GGGGGCGGCCATTc                                                                                                                                          <  2:1122937/150‑1 (MQ=255)               ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATc                                                                                                                                          <  1:2243889/150‑1 (MQ=255)               ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATc                                                                                                                                          <  1:2243883/150‑1 (MQ=255)                gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTGTCTCATCCTGTTGCGGTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Gg‑g‑ggggggggggcggtg                                                                                                                                         >  1:463317/1‑143 (MQ=255)                gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGcccaaag                                                                                                                                         >  2:4161140/1‑145 (MQ=255)                      ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATcgccgccg                                                                                                                                   <  1:2320405/150‑1 (MQ=255)                      ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATcgccgccg                                                                                                                                   <  1:2531969/150‑1 (MQ=255)                       tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGGCCAAcgccccccc                                                                                                                                  >  1:2716597/1‑148 (MQ=255)                         aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCAAAAGCCGcccccc                                                                                                                                >  1:2513095/1‑150 (MQ=255)                            ttatatTTTAAAATTAATTAATGTCATCAGGTCCTAAAATAACGAGAATATTTCAGTCTCTAATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                             <  1:2800082/148‑1 (MQ=255)                            tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                             <  2:2145389/150‑1 (MQ=255)                             gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCccccccccccag                                                                                                                            >  2:3563832/1‑148 (MQ=255)                              atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGc                                                                                                                           <  2:3815194/150‑1 (MQ=255)                                 ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑CGGGGGGGCCCAACgcccccgcccccaccccc                                                                                                                        >  2:2261187/1‑144 (MQ=255)                                  tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑G‑GGGGGGGGCCAA‑CGCccccccccccccccccc                                                                                                                       >  1:2907562/1‑143 (MQ=255)                                  tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑ACGGGGGGGGCCCAACccccccccccccccaccc                                                                                                                        >  2:1580423/1‑144 (MQ=255)                                  tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑G‑GGGGGCGGCCCAAACCCCCcgccccccgccccc                                                                                                                       >  2:1993700/1‑150 (MQ=38)                                  tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGGGGGCCACAGCCGccccccccccccccc                                                                                                                       >  1:3671054/1‑143 (MQ=255)                                  tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAAAAGCCGCCCCCCCCACCAccc                                                                                                                       >  1:713574/1‑150 (MQ=38)                                   ttAATATTATTTATTGTATTCTGGTCCTAAAATAAGGATAATATTTCAGTCTCTATTCTTGTTGCGCTCTTTTCTTGTGCTTTGCTTAATATAATCCCTGGCGCATGTTGCGG‑GC‑T‑GGGGGGGGCAATTCGCCGCCGCCCCCTGCAcccc                                                                                                                      <  1:2559390/150‑1 (MQ=255)                                   ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑C‑GGGGGGGGCCAAAAGACCCCCCCCCCAGcccccc                                                                                                                      >  1:1710326/1‑150 (MQ=255)                                    tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                     <  2:1727546/150‑1 (MQ=255)                                     aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                    <  2:44329/150‑1 (MQ=255)                                       aattaattAATGTCATCATGTCCTAAAATAACTAGAATATTTCATTCTATCATCCTTTTGCGATCATGTAATGTGAATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:1593741/150‑1 (MQ=25)                                        attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGc                                                                                                                 <  2:1359933/150‑1 (MQ=255)                                        attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGc                                                                                                                 >  2:4270436/1‑150 (MQ=255)                                          taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCGGGctc                                                                                                               >  2:1429668/1‑150 (MQ=255)                                           aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                              <  2:1711149/150‑1 (MQ=255)                                             ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                            <  2:2557756/150‑1 (MQ=255)                                                  gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑A‑GGGGGGGGCCCAGAGCCCCCCCCCCCCCCACCCCCCGGCACAGG‑CGaaa                                                                                                      >  1:3565999/1‑149 (MQ=37)                                                   tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGG‑GGGGGGGGCCAAAC‑CCCCCCCCCCCCGCCCCCCCCGGCGCTGG‑Cgccca                                                                                                     >  1:2166465/1‑146 (MQ=37)                                                     atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑GGGGGGGGGGCATAGGCGCCGCCCCCAGCACCCCCCGGCGCTGG‑GGGACaaa                                                                                                   >  1:968278/1‑150 (MQ=255)                                                      tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑G‑GGGGGCGGCCCATCGCCGCCGCCCCCCGCACCCCCCGGCGCAGG‑CGAACaaaa                                                                                                  >  2:2294264/1‑150 (MQ=25)                                                      tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACaaaa                                                                                                  <  2:1143985/150‑1 (MQ=255)                                                      tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTGGCGCTCCTGTCATGTGCATTGCTTCATATAAGCACTGGCGCACGGAGCGC‑GC‑C‑GGGGGCGGCCCATCGCCGCCGCCCCCAGCCCCCCCCGGCGCGGG‑CGCACcaaa                                                                                                  >  2:4232063/1‑150 (MQ=255)                                                       atgttCCTAAAATAACGAGAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGATTTGCTTCATATAACCACTGGCGCATGGCGCGC‑GC‑T‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAAt                                                                                                 <  1:2607061/146‑1 (MQ=18)                                                        tggTCCGAAAATAACGATATTTTTTCTGCCCCTCCTCGTGTGGGGGTCCTGTCTTGTGCTTTGCTTCATATACTCCGGGGCGCTCGGTGCGG‑GC‑G‑GGGGGGGGCCAAGCGCCGCGGCCCCCGGCACCCCCGGGCTCTGG‑CGAACAAAATc                                                                                                <  1:4053383/149‑1 (MQ=12)                                                           tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑GGGGGGGGCCCAACGCCCCCCCCCCCCCCCCCCCCCCCCGCGGG‑CGCACaaaaaacac                                                                                             >  1:2016939/1‑145 (MQ=255)                                                           tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑GGGGGCGGGCCAAAGCCCCCCCCCCCCCCCCCCCCCGGCCCAGGCCGACCaaaacacc                                                                                              >  1:2485677/1‑146 (MQ=255)                                                               aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCCAAAGCCGCCGCCCCCCGCACCCCCCGGCGCAGG‑CGAACCAAATcgccgcgg                                                                                         >  2:712144/1‑148 (MQ=25)                                                                  aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTgcgc                                                                                      <  2:3647427/150‑1 (MQ=255)                                                                        agaATATTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGcc                                                                                <  2:3592345/150‑1 (MQ=255)                                                                             atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATACTCCCTGGCGCACGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCTCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:3211912/150‑1 (MQ=255)                                                                              tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGCTCTGG‑CGAACAAAATAGCCGCTGCGCGGTGCCCACGGc                                                                          >  2:4286267/1‑150 (MQ=18)                                                                                    tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:4215956/150‑1 (MQ=255)                                                                                        tcACCCTGTTGCGCTCTTGTCTTGTGATTTGATTAATATAATCCGTGGCGCATGGCGGGC‑GC‑G‑GGGGGGGCCCACGCGCCGCCGCCCCCGGCACCCCCGGGCTCGGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:873141/150‑1 (MQ=14)                                                                                         cagccTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                               <  1:2728833/147‑1 (MQ=255)                                                                                         cagccTGTTGCGCTCCTGTCATGGGCAGTGCTTCATATAATCCCGGGCGCAAGGCGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCGGCGCGGTGCCCTCGGCTTATCCCTTAc                                                               <  2:1396446/147‑1 (MQ=17)                                                                                                  gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  2:84962/150‑1 (MQ=255)                                                                                                      tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑CGGGGCGGCCCATCGCCCCCCCCCCCAGCACCCCCCGGCGCTGG‑CGAACAAAAACGCCCCTGCGCGGTGCCCACCGCGTATCCCATACGGCGACCCGGTCg                                                  >  2:701321/1‑150 (MQ=17)                                                                                                        cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCACCCCCAGCACCCCCGGGCACTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                >  2:3536235/1‑150 (MQ=255)                                                                                                             aTGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:3764991/150‑1 (MQ=255)                                                                                                             aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCCACCGGGTCGGGcgcg                                            >  2:1526059/1‑150 (MQ=255)                                                                                                             aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                            >  2:1526129/1‑150 (MQ=255)                                                                                                             aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAAACGCCGCCGCCCCCAGCCCCCCCCGGCGCAGG‑CGAACAAAAAAACCCCCGCGCGGTGCCCGCCGCATATACCATAAGGCAAACCGGTCGGGcgcgc                                           >  2:3795276/1‑149 (MQ=12)                                                                                                              tgtgCATTGCTTCATATACTCACTGGCGCAAGGCGCGC‑GC‑T‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:1615903/150‑1 (MQ=255)                                                                                                               gtgCATTGCTTCATATAATCCCGGGCGCACGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:502846/150‑1 (MQ=255)                                                                                                               gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:502906/150‑1 (MQ=255)                                                                                                               gtgAATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  2:814049/150‑1 (MQ=255)                                                                                                                 gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                       <  2:3619999/150‑1 (MQ=255)                                                                                                                    ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:463317/150‑1 (MQ=255)                                                                                                                       cTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:3908032/150‑1 (MQ=40)                                                                                                                        ttCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAAAAGCCCCCCCCCCCAGCACCCCCCGGCGCAGG‑CGGACAAAAAAGCCGCGGCGCCGTGCCCGCCGCGTATCCCCTACGGCGACCCGGGCGGGCGCGCGGGTACAAccc                                >  2:246915/1‑150 (MQ=9)                                                                                                                         tCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:3563832/150‑1 (MQ=39)                                                                                                                         tCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:4161140/150‑1 (MQ=39)                                                                                                                         tCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:332130/150‑1 (MQ=39)                                                                                                                                aTCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:2261187/150‑1 (MQ=33)                                                                                                                                aTCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:1558536/150‑1 (MQ=33)                                                                                                                                     tGGCGCATGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                   <  1:4286267/150‑1 (MQ=17)                                                                                                                                          cAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:1197719/150‑1 (MQ=32)                                                                                                                                            aGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:569882/150‑1 (MQ=32)                                                                                                                                            aGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:569855/150‑1 (MQ=32)                                                                                                                                            aGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  2:2694458/150‑1 (MQ=32)                                                                                                                                             ggAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  1:2713122/150‑1 (MQ=32)                                                                                                                                                gcgc‑gc‑T‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca        <  1:1833797/150‑1 (MQ=17)                                                                                                                                                     gc‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa    <  1:2759127/150‑1 (MQ=32)                                                                                                                                                     gc‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa    <  1:2759122/150‑1 (MQ=32)                                                                                                                                                        a‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:1528776/1‑150 (MQ=32)                                                                                                                                                        |                                                                                                                                                        TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

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