Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,222,656	A→G	20.8%	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
*	CP009273	2,222,656	0	A	G	20.8%	93.3 / 17.4	85	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (4/46);  new base G (10/3);  total (19/66)
Fisher's exact test for biased strand distribution p-value = 1.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCA  >  CP009273/2222512‑2222803                                                                                                                                                   |                                                                                                                                                      gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑ggg                                                                                                                                                   <  2:1540132/148‑1 (MQ=255) gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑ggg                                                                                                                                                   <  2:601562/148‑1 (MQ=255) gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑ggg                                                                                                                                                   <  2:1540108/148‑1 (MQ=255) gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑ggg                                                                                                                                                   <  2:2793404/148‑1 (MQ=255) gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑ggg                                                                                                                                                   <  2:1540101/148‑1 (MQ=255)     tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑gg‑g‑ggggggg                                                                                                                                               >  2:2263836/1‑146 (MQ=255)     tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑gg‑g‑ggggggc                                                                                                                                               >  2:2263843/1‑146 (MQ=255)        cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Gg‑g‑gggggggggg                                                                                                                                            >  1:989099/1‑146 (MQ=255)        cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Gg‑g‑gggggggggg                                                                                                                                            >  1:989128/1‑146 (MQ=255)          ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCaa                                                                                                                                          <  1:2563044/148‑1 (MQ=255)              tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGG‑G‑GGGGGCGGcccaaca                                                                                                                                       >  2:3584338/1‑143 (MQ=255)              tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑C‑GGGGGGGCCCaaccc                                                                                                                                       >  2:3583785/1‑145 (MQ=255)               aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑G‑GGGGGGGGcccaagcc                                                                                                                                      >  2:2921606/1‑142 (MQ=255)                    aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATcgccgccgc                                                                                                                                <  2:3537590/148‑1 (MQ=255)                     gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑G‑GGGGGGGGCCCA‑Cgccgccccc                                                                                                                               >  2:2333839/1‑148 (MQ=255)                     gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑C‑GGGGGGGGCCCCACgccgccccc                                                                                                                               >  2:2333818/1‑148 (MQ=255)                     gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑G‑GGGGGGGGCCCAAAgccgccccc                                                                                                                               >  2:2462109/1‑148 (MQ=39)                     gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGATCCTGTCATGTGCATTGCTTCATATAATCAATGGCGCAAGGAGCGC‑GC‑A‑GGTGTCTGCCAATCgcctccgcc                                                                                                                               <  2:45535/148‑1 (MQ=255)                       tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑CGGGGCGGCCCCGACCCGCCGcccc                                                                                                                             >  2:2476341/1‑148 (MQ=255)                       tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGcccc                                                                                                                             <  2:345885/148‑1 (MQ=255)                        gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑C‑GGGGGGGGCGACTCcccccccccccc                                                                                                                            >  2:2350742/1‑148 (MQ=255)                        gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGccccc                                                                                                                            <  2:1618911/148‑1 (MQ=255)                            ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGG‑GGGGGGGGCCAAAGCCCGccccccccccc                                                                                                                         >  2:60578/1‑145 (MQ=255)                                   attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑G‑GGGGGGGGCCAA‑CCCCGCccccccccc‑ccccccccg                                                                                                                >  2:3299159/1‑146 (MQ=255)                                   attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑GGGGGGGGGCCCCCGccccccccccccc‑ccccccccc                                                                                                                >  2:2286688/1‑146 (MQ=255)                                   attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑C‑GGGGGGGGGCACGAGCCGCccccccccc‑ccccccccg                                                                                                                >  2:2110308/1‑146 (MQ=255)                                    ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑GGGGGGGGGCCGCCGccccccccccccc‑cccccccccc                                                                                                               >  1:2269583/1‑145 (MQ=255)                                     taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGc                                                                                                              <  2:1260669/148‑1 (MQ=255)                                      aattaatGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGct                                                                                                             <  2:185472/148‑1 (MQ=255)                                                  caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑G‑GGGGGGGGGCCCACGCCCCCCCCCCCCC‑CACCCCCCCCCCCGGGCGAACaaa                                                                                                 >  2:30471/1‑148 (MQ=255)                                                  caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACaaa                                                                                                 <  2:545227/148‑1 (MQ=255)                                                   tggTCCGAAAATAACTAGAATATTTCAGTCTCTCATCCTGTTGGGCTCCTGTAATGTGCTTTGCTTCATATAATCAGTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCAAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACaaaa                                                                                                <  1:2525767/147‑1 (MQ=21)                                                   aGGTCCTAAAATAACTATAATATTTCAGTCTCTCATCGTGTTGGGCTCCTGTCTTGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGC‑GC‑T‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACaaaa                                                                                                <  1:2817000/148‑1 (MQ=21)                                                   aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACaaaa                                                                                                <  2:3735792/148‑1 (MQ=255)                                                    ggTCCTAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAAt                                                                                               <  2:2309177/148‑1 (MQ=255)                                                      tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑AGGGGGGGGGCCAAAC‑CCCCCCCCCCCCG‑CACCCCCCCGCCCCGGCGCACaaaaaac                                                                                             >  2:948546/1‑145 (MQ=255)                                                       ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑C‑GGGGGGGGGCCCGCGCCCCCCCCCCCCG‑CCCCCCCCCGCCCGGGCGGACaaaaacac                                                                                            >  2:2584115/1‑144 (MQ=255)                                                       ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                            <  2:731258/148‑1 (MQ=255)                                                         gAAAATAACGAGAATATTTCAGTCTCTCTTCCTTTTGCGCGCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGC‑GC‑T‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATagccg                                                                                          <  2:334766/148‑1 (MQ=25)                                                         gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATagccg                                                                                          <  2:334782/148‑1 (MQ=255)                                                         gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATagccg                                                                                          <  2:79657/148‑1 (MQ=255)                                                            aaTAACGATATTTTTTCAGTCTCTCATCCTGTTGCGCTCGTGTAATGTGCATTGATTAATATAATCCCTGGCGCACGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                       <  2:1018908/148‑1 (MQ=18)                                                            aaTAACGATAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTAATATAAACACTGGCGCAAGGAGGGC‑GC‑T‑GGGGGGGCCCAATCGCCGCGGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                       <  2:1019068/148‑1 (MQ=18)                                                               aacgtcaaTATTTCAGTCTCTCCTCCTTTTGCTCTCCTGTCATGTGCATTGCTTCATATAATCAGTGGCGATAGGAGCGC‑GC‑A‑GGGGGAGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcgc                                                                                    <  2:3596492/142‑1 (MQ=21)                                                                 cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCACGGGCTCTGGCGAACAAAATCGCCGCTGCGCgg                                                                                  >  2:2491174/1‑148 (MQ=255)                                                                           tCAGTCTCTCATCCTGTTGCGCTCGTGTCATGTGCTTGGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAACCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                        <  2:1352244/148‑1 (MQ=255)                                                                           tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                        <  2:1352283/148‑1 (MQ=255)                                                                           tCAGTCTCTCATCCTGTGGCGCTCCTGTCTTGTGCTTTGCTTCATATAATCCCGGGCGCACGGAGCGC‑GC‑A‑GGGGGCGGCAAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                        <  2:1333866/148‑1 (MQ=18)                                                                           tCAGTCGCTCAGCCTGTTGCGCGCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                        <  2:2623804/148‑1 (MQ=255)                                                                                  ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                 <  2:3068142/148‑1 (MQ=255)                                                                                  ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                 <  2:3057597/148‑1 (MQ=255)                                                                                  ctcATCCTGTGGCGCTCCTGTCTTGTGATTTGATTAATATAATCCGTGGCGCAAGGCGCGC‑GC‑T‑GGGGGCGGCCAAGCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                 <  1:951892/148‑1 (MQ=17)                                                                                  ctcATCCGGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                 <  2:839682/148‑1 (MQ=255)                                                                                  ctcAGCCTGGTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                 <  2:3068036/148‑1 (MQ=255)                                                                                   tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                <  2:1516494/148‑1 (MQ=255)                                                                                    cATCGTGTTGCGCTCGTGTCTTGTGCTTTGCTTCATATAAGCAGGGGCGCCCGGCGCGC‑GC‑C‑GGGGGCGCCCAAGCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                               <  1:210660/148‑1 (MQ=17)                                                                                    cATCGTGTTGCGCTCCTGTCTTGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑C‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                               <  1:3400078/148‑1 (MQ=17)                                                                                    cATCGTGTGGCGCTCCTTTCTTGTGATTTGTTTAATATACTCAGGGGCGCACGGGGCGC‑GC‑G‑GGGGGCGCCCAACCCCCCCCGCCCCCGG‑CACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                               <  1:1798437/148‑1 (MQ=9)                                                                                    cATCCTGTGGCGCTCCTGTCATGTGCAGTGCTAAATATAACCAGGGGCCCACGGCGCGC‑GC‑C‑GGGGGCGCCAAACCGCCCCCCCCCCCCC‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                               <  1:1408352/148‑1 (MQ=14)                                                                                       ccTGTTGCGGTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGGCAATAGCCGCCCCCCCCCGCCCCCCCCGGGCGCTGGGGAACAAAAACCCCGCGGCGCGGTGCCCACGGCGTATCCGTTAc                                                             >  2:842888/1‑148 (MQ=18)                                                                                        gtGTTGGGCTCTTGTCTTGTGATTTTATTAATATAACCAGTGGCGCATGGTGCGC‑GC‑T‑GGGGGGGGCCAATCGCCGCCGCCCCGTG‑CACCCCCGGGATCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                           <  1:2406934/147‑1 (MQ=11)                                                                                               gcTCGTGTCATGTGATTTGCTACATATAATCACTGGCGCCTGGAGCGC‑GC‑T‑GGGGGCGCCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                    <  1:1506745/148‑1 (MQ=14)                                                                                                cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                   <  2:3691934/148‑1 (MQ=255)                                                                                                cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                   <  2:3691910/148‑1 (MQ=255)                                                                                                 tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                  <  2:2684245/148‑1 (MQ=255)                                                                                                        aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcgc                                           <  2:2843602/148‑1 (MQ=255)                                                                                                         tgtgAATTGCTAAATATAACCACGGGCCCATGGTGCGC‑GC‑T‑GGGGGCGGCCAAGCGCCGCCCCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                          <  1:2095134/148‑1 (MQ=11)                                                                                                           tgatttGCTTCATAAAACCCCGGGCCCCGGGCGCGC‑GC‑G‑GGGGGCGCCCAAGCGCCGCCCCCCCCGC‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                        <  1:2921606/144‑1 (MQ=11)                                                                                                            gctttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                       <  2:494889/145‑1 (MQ=255)                                                                                                            gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                       <  2:3748065/148‑1 (MQ=255)                                                                                                             aaTTGATTCATATAATCACTGGCGCATGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                      <  2:2022684/147‑1 (MQ=255)                                                                                                              aTTGCTTCATATAATCACTGGCGCACGGCGCGC‑GC‑T‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                     <  1:292528/148‑1 (MQ=255)                                                                                                               tggCTTAATATAACCCCGGGCGCACGGCGGGC‑GC‑G‑GGGGGCGCCCAAGCCCCGCCCCCCCCGC‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                    <  1:1095219/146‑1 (MQ=11)                                                                                                                tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                   <  2:212686/148‑1 (MQ=255)                                                                                                                     cATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                              <  2:1202193/148‑1 (MQ=38)                                                                                                                              aggggCCCATGGCGCGC‑GC‑T‑GGGGGCGCCAACCCGCCGCCGCCCCCGC‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                     <  1:3551142/145‑1 (MQ=11)                                                                                                                                   cgcAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                <  1:934701/148‑1 (MQ=32)                                                                                                                                    gcACGGAGCGC‑GC‑T‑GGGGGCGGCCACGCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt               <  1:1511267/148‑1 (MQ=11)                                                                                                                                    gcAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt               <  2:758181/148‑1 (MQ=32)                                                                                                                                         gAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc          <  2:2012486/148‑1 (MQ=32)                                                                                                                                         gAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc          <  2:3141558/148‑1 (MQ=32)                                                                                                                                          aGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc         <  2:645478/148‑1 (MQ=32)                                                                                                                                          aGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc         <  2:2052140/148‑1 (MQ=32)                                                                                                                                             gc‑gc‑T‑GGGGGGGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca      <  1:544663/148‑1 (MQ=17)                                                                                                                                             gc‑gc‑T‑GGGGGCGGCCAAGCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca      <  1:2462109/148‑1 (MQ=17)                                                                                                                                             gc‑gc‑C‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca      <  1:2462596/148‑1 (MQ=17)                                                                                                                                             gc‑gc‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca      <  2:1341636/148‑1 (MQ=32)                                                                                                                                              c‑gc‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt     >  2:607759/1‑148 (MQ=32)                                                                                                                                                gc‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc    <  2:510105/148‑1 (MQ=32)                                                                                                                                                   a‑GGGGGCGGCCAAGCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa  <  1:2476349/148‑1 (MQ=17)                                                                                                                                                   |                                                                                                                                                      GTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑A‑GGGGGCGGCCAATCGCCGCCGCCCCCTG‑CACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCA  >  CP009273/2222512‑2222803

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