Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,111,727	2 bp→TT	coding (865‑866/2544 nt)	opgH →	osmoregulated periplasmic glucans biosynthesis protein H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,111,727	0	C	T	88.9%	147.3 / 9.6	54	P289M (CCG→TTG)	opgH	osmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (29/19);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01
*	NC_000913	1,111,728	0	C	T	88.7%	141.4 / 10.0	54	P289M (CCG→TTG)	opgH	osmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (29/18);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

GATTCTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCC  >  NC_000913/1111485‑1111969                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  gATTCTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATc                                                                                                                                                                                                                                            >  2:154853/1‑251 (MQ=255)              gTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGc                                                                                                                                                                                                                               >  2:307341/1‑251 (MQ=255)               ttAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGc                                                                                                                                                                                                                               >  1:291826/1‑250 (MQ=255)                tAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCa                                                                                                                                                                                                                              <  2:186061/250‑1 (MQ=255)                             gCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCtggatgg                                                                                                                                                                                                                >  1:179015/1‑251 (MQ=255)                                       aTTGGGCGCTATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGGTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCg                                                                                                                                                                                                       <  1:260750/250‑1 (MQ=255)                                           gTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGaa                                                                                                                                                                                                  <  1:540320/251‑1 (MQ=255)                                               cgATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCg                                                                                                                                                                                              >  1:326089/1‑251 (MQ=255)                                                gATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCgg                                                                                                                                                                                             >  1:528901/1‑251 (MQ=255)                                                 aTAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACCCTTTGATGTCTACATTCTTAGTGACAGTTATACCTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGt                                                                                                                                                                                            >  2:260313/1‑251 (MQ=255)                                                         aGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAgg                                                                                                                                                                                    >  1:548399/1‑251 (MQ=255)                                                            atatCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCa                                                                                                                                                                                 <  2:261455/251‑1 (MQ=255)                                                               tCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGTAGGTc                                                                                                                                                                                  >  2:315005/1‑247 (MQ=255)                                                                 tGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGGCAGAttt                                                                                                                                                                            >  2:451427/1‑251 (MQ=255)                                                                          aGTTGGCGATGAACCATTAAACCCGGGGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGACAGATTTTCTATcgcc                                                                                                                                                                   >  2:190066/1‑251 (MQ=255)                                                                           gTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATcgccg                                                                                                                                                                  <  1:442186/251‑1 (MQ=255)                                                                           gTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCATGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATcgccg                                                                                                                                                                  >  1:350286/1‑251 (MQ=255)                                                                                     aaCCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTcgc                                                                                                                                                        >  1:313534/1‑251 (MQ=255)                                                                                      aCCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTcgcc                                                                                                                                                       >  2:545445/1‑251 (MQ=255)                                                                                      aCCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTcgcc                                                                                                                                                       <  1:452202/251‑1 (MQ=255)                                                                                         aTTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCcgcg                                                                                                                                                    >  1:147524/1‑251 (MQ=255)                                                                                                  ggAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGt                                                                                                                                            >  1:188395/1‑250 (MQ=255)                                                                                                    aGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAACCACTTTGATGTCTACATTCTTAGTGACAGTTATAATTTGGATATCTGCGTCGCAGAGAAAAAGGCCTGAATGGAGCTTATCGCAGAAGTAGGGGGGGAGGGTAGGATTTTCTTACGACGCCGCCGTGGCCGAGTGAAGCgaaa                                                                                                                                          >  2:375282/1‑247 (MQ=255)                                                                                                          gCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTAGCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGt                                                                                                                                   <  2:73112/251‑1 (MQ=255)                                                                                                           cACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGt                                                                                                                                   >  1:23830/1‑250 (MQ=255)                                                                                                             cGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAat                                                                                                                                <  1:154853/251‑1 (MQ=255)                                                                                                                         tGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTc                                                                                                                    >  2:316340/1‑251 (MQ=255)                                                                                                                           ccTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTg                                                                                                                  >  1:208601/1‑251 (MQ=255)                                                                                                                                  gTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCg                                                                                                               <  1:396152/247‑1 (MQ=255)                                                                                                                                  gTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGc                                                                                                           >  1:159499/1‑251 (MQ=255)                                                                                                                                     aCGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCGGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTgg                                                                                                        <  2:368900/251‑1 (MQ=255)                                                                                                                                        aaGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGCTGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGc                                                                                                     >  1:489872/1‑251 (MQ=255)                                                                                                                                        aaGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGc                                                                                                     <  1:316340/251‑1 (MQ=255)                                                                                                                                                              tGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGCAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtggtgg                                                                               >  2:155266/1‑251 (MQ=255)                                                                                                                                                              tGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtggtgg                                                                               >  1:416005/1‑251 (MQ=255)                                                                                                                                                                          aaCGGGGGAACCAGTAAAAGCCACCGGTAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTAGAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCGGGATGGAGCGTATCGCTGAAGTCGGTGGCGAAGGTCATATGTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGTCTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTg                                                                   <  2:40866/251‑1 (MQ=255)                                                                                                                                                                             gTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTAGCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACt                                                                <  2:326089/251‑1 (MQ=255)                                                                                                                                                                                   aTCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGAGGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTaa                                                          <  2:528901/251‑1 (MQ=255)                                                                                                                                                                                     aaGTAAAAGCCACCGGGGATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCAGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTCATg                                                        >  2:116450/2‑251 (MQ=255)                                                                                                                                                                                      aGTAAAAGCCACCGGGAATGCCAAACACATTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGa                                                       <  1:104812/251‑1 (MQ=255)                                                                                                                                                                                        tAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATAACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGAcc                                                     >  1:384155/1‑251 (MQ=255)                                                                                                                                                                                        gaaaaGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGAcc                                                     <  1:180181/250‑1 (MQ=255)                                                                                                                                                                                              ccACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAt                                               <  1:556179/251‑1 (MQ=255)                                                                                                                                                                                                aCCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTAc                                                                                      <  2:225509/210‑1 (MQ=255)                                                                                                                                                                                                aCCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTAc                                                                                      >  1:225509/1‑210 (MQ=255)                                                                                                                                                                                                     gAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAttgtttgt                                        <  2:159499/251‑1 (MQ=255)                                                                                                                                                                                                        tGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCg                                     >  2:42074/1‑251 (MQ=255)                                                                                                                                                                                                            aaaCACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCt                                 <  1:540914/251‑1 (MQ=255)                                                                                                                                                                                                                acTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGt                                                                                              <  2:155575/186‑1 (MQ=255)                                                                                                                                                                                                                acTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGt                                                                                              >  1:155575/1‑186 (MQ=255)                                                                                                                                                                                                                  tgtGATGTCGAAATTCTGAGTGAAAGTGATAATTGGGATAGCTGGGTAGCAGAGCAAAAAGCATGGAGGGAGCTTAGCGATGAAGTAGGTGGCGAAGGTCAGATTGTCTACCGCTGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGAGGACTGCTGCCGTCGCTGGGGCAGGCAGTATAGTTAAATGGTGGTGCTGGATGCTGACTCGGTAATGATCGGTGATTGTTTGCGCGGGCTGGTgc                            <  2:440941/248‑1 (MQ=255)                                                                                                                                                                                                                                   tAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAAc           <  1:463579/250‑1 (MQ=255)                                                                                                                                                                                                                                         cAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGCTGGAAGCCACCCCGAACg    >  2:168997/1‑251 (MQ=255)                                                                                                                                                                                                                                           gTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGcc  >  1:142262/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  GATTCTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCC  >  NC_000913/1111485‑1111969

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