Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,112,209	C→T	D449D (GAC→GAT)	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,112,209	0	C	T	96.8%	186.0 / ‑2.5	62	D449D (GAC→GAT)	opgH	osmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (34/26);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

GAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGACTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATG  >  NC_000913/1111963‑1112449                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   gAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCg                                                                                                                                                                                                                                                >  1:290477/1‑251 (MQ=255)  aaCGCCGGGATCATTCAGTCTTCGCCGAAAGCGTCCGGTATGGATTCGCTGTATGCGCGCTGTCAGCAGTTCGCTACCCGCGTGTATGGGCCACTGTTTACAGTCTGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCTTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTc                                                                                                                                                                                                                                                 <  1:304782/249‑1 (MQ=255)              agtCAGTCGTCGCCGAAAGCGCCCGGTATGGATACGCTGTATGCGCGCTGCCAGCAGTGCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGtt                                                                                                                                                                                                                                   <  2:234464/249‑1 (MQ=255)               ttCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGttt                                                                                                                                                                                                                                  >  2:274813/1‑250 (MQ=255)                 cAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGCAGCGGCGTTGAt                                                                                                                                                                                                                                >  1:224833/1‑251 (MQ=255)                  aGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATg                                                                                                                                                                                                                               >  1:12857/1‑251 (MQ=255)                   gtcgtcGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGTAAGCGGCGTTGATgc                                                                                                                                                                                                                              <  1:415664/251‑1 (MQ=255)                        cgCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCgg                                                                                                                                                                                                                         >  2:160347/1‑250 (MQ=255)                        cgCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGt                                                                                                                                                                                                                         >  1:364198/1‑251 (MQ=255)                              aaGCGTCCGGTACGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGt                                                                                                                                                                                                                   <  1:345595/251‑1 (MQ=255)                                       gTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑Tc                                                                                                                                                                                                         >  2:89574/1‑251 (MQ=255)                                       gTATGGATACGCTGTATGCGCGATGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATACCGCGATTATCCGCGTGAAACCGTTTATCGAGCAATGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGACGGTTCAATCCTGTCACATGATTTCGTGGAAGCGACGTTGATGCGCCGTGAGGGTTGGGGGG‑Tc                                                                                                                                                                                                         >  1:484245/1‑251 (MQ=255)                                        tATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCCCATGATTTCTTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCt                                                                                                                                                                                                        >  1:34996/1‑251 (MQ=255)                                         aTGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGACTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTg                                                                                                                                                                                                       <  2:170334/251‑1 (MQ=255)                                             aTACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGAtt                                                                                                                                                                                                   <  1:454703/251‑1 (MQ=255)                                               aCGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCg                                                                                                                                                                                                                                                >  1:82758/1‑205 (MQ=255)                                               aCGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCg                                                                                                                                                                                                                                                <  2:82758/205‑1 (MQ=255)                                                        gcgcgcTGTCAGCAGTTGGCGACCCGCGTGTATGGGCCCCTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCAGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGAtct                                                                                                                                                                                       <  2:162563/251‑1 (MQ=255)                                                           cgcTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTccc                                                                                                                                                                                    >  1:352047/1‑251 (MQ=255)                                                                 cagcagTTCGCGCCCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGtt                                                                                                                                                                               <  1:410702/250‑1 (MQ=255)                                                                 cagcagTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGTCAACTTAGCGAGTCGCACTACTGGGGGCATAACGCGATTATACGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCACCGCCGCCGGGCGAAGGATC‑CTT‑TGCAGGTTCAATCCTATCACCTGATTTCGTGGAAACGGCGTTGTTGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTc                                                                                                                                                                              >  1:435093/1‑251 (MQ=255)                                                                  ggcggTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCTCTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCt                                                                                                                                                                             <  2:80736/247‑1 (MQ=255)                                                                         cgcgACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATga                                                                                                                                                                        >  1:242425/1‑249 (MQ=255)                                                                                cgcgTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTgcc                                                                                                                                                               <  2:137724/251‑1 (MQ=255)                                                                                    tgtATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCAAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGCGGAATTGCCGCCt                                                                                                                                                           >  2:138886/1‑251 (MQ=255)                                                                                       aTGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTaa                                                                                                                                                         <  1:334599/250‑1 (MQ=255)                                                                                          ggCCACTGTTTTCAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAActtg                                                                                                                                                     <  2:377812/251‑1 (MQ=255)                                                                                          ggCCACTGTTTTCAGCCGGTTGGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAATAATTGCCGCCTAActtg                                                                                                                                                     <  2:54938/251‑1 (MQ=255)                                                                                          ggCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAActtg                                                                                                                                                     <  1:89574/251‑1 (MQ=255)                                                                                           gCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTACGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAActtgc                                                                                                                                                    >  1:342623/1‑251 (MQ=255)                                                                                                      cAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCt                                                                                                                                          >  2:135691/1‑250 (MQ=255)                                                                                                       aGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCGCTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAActtgct                                                                                                                                                   <  2:465976/240‑1 (MQ=255)                                                                                                       aGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAActtgct                                                                                                                                                   >  1:465976/1‑240 (MQ=255)                                                                                                           ggTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACg                                                                                                                                    >  1:52042/1‑251 (MQ=255)                                                                                                            gTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCTTAACTTGCTTGATGAGCTAAAACGt                                                                                                                                   <  2:409448/251‑1 (MQ=255)                                                                                                            gTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGt                                                                                                                                   >  2:186435/1‑251 (MQ=255)                                                                                                                gccattCTGGCAACTTGGCGAGTCGCACTACCGGGGACATCACGCGATTATCCGCGTGAAACCGTTTATCGAGCCCTGCGCACGGGCTCCGCTGCCGGGCGAAGGTTC‑CTG‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGAc                                                                                                                                <  2:302284/246‑1 (MQ=255)                                                                                                                gCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACGGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGAc                                                                                                                                <  2:412492/250‑1 (MQ=255)                                                                                                                           aaCTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACGGCGCACTGGCTCCGCTGCCGAGCGAAGGTTC‑CTT‑TGCCGGTTCACACCTGTCACATGATTTCGTGGAAg                                                                                                                                                                                                                                          >  1:200861/1‑135 (MQ=255)                                                                                                                           aaCTTGGCGAGGCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGGTTATCGAGCTCTGCGTACTTGCTCCGCTGCCGGGCGAAGGTTT‑CTT‑TGCCGGTGTAATCCTGTCACATGATTTAGTGGAAg                                                                                                                                                                                                                                          <  2:200861/135‑1 (MQ=255)                                                                                                                                ggCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACgg                                                                                                               >  1:182088/1‑251 (MQ=255)                                                                                                                                               tGGGGACATAACGCGATTATCCGCGTGAAACAGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCGAACTTGATTGATGAGCTAACACGTGACCGCCGATGGTGCCACGGTAACCTGATGaaatt                                                                                                >  2:394295/1‑248 (MQ=255)                                                                                                                                                 gggACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTcc                                                                                              >  1:383722/1‑251 (MQ=255)                                                                                                                                                   gACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGt                                                                                            >  1:142797/1‑251 (MQ=255)                                                                                                                                                    aCATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTc                                                                                           >  2:424426/1‑251 (MQ=255)                                                                                                                                                         ccgtgATTACCCGGGTGAACCCGTTTATGGAGCCCTGCGCACGGGCCCCGCGGCCGGGGGAGGTTTC‑TTT‑TGCTGGTTCACCCCTGCCTCATGATTTCGGGAAAGCGGCGTGGAGGCGCCGTGCGGGTGGGGGTG‑TTTGGATGTCTT‑CCGATCTCCGGGTTTTTTAGGAAGATTTGCCGCCTAACTTGCTTGATGAGATAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTc                                                                                      <  1:481661/247‑1 (MQ=255)                                                                                                                                                                       gTGAAACCGTTTATCGAGCACTGCGCACTGGCTCAGCTGCCGGGCGAAGGGTA‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑ACTGGATTGCTT‑ACGATCTACCGGGATCTTATGCCGAATTGCCGCCTAAATTGCTTGATGAGCTAAAACGGGACCGCAGATGGGGCCACGGTCACCTGGTGAACTTCCGTCTGTTCTTGGGGAAGGGTAt                                                                        >  2:270446/1‑251 (MQ=255)                                                                                                                                                                        ggaacccGTTTACCGAGCCCTGCGCACTGGCCCGGCTGCCGGGCGAAGTTTC‑TTT‑TGCCGTTTCAATCCTGCCAAATGTTTTCGGGGAGGCGGCGTGGATGCGCGGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCCCCCGGGTTTTAATGAAGAATTGCCGCCAAACTTGCTTGAGGAGCAAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATg                                                                       <  1:274813/246‑1 (MQ=255)                                                                                                                                                                           aacggTTTATCGAGCCCTGCGACCTGGTTCCGCGGGGGGGCGAAGGTTC‑CTGGTGC‑GGTTCAATCCTGAAACATGATTCCGTGGAAGCGGTGTGGTTTCGCCGGGCAG‑TTTGGGGG‑TCCGGTTTTCTTTACAATCTCCCGGTTTCTAATGAAGAATGCCCGCCTAACTGGCTAGTTGAGTAAAAACGTGTCCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCAc                                                                    <  1:54662/247‑1 (MQ=255)                                                                                                                                                                                tttATCGAGCACTGCGCACTGGCTCCGCTGCTGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGACTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGt                                                               >  1:194918/1‑251 (MQ=255)                                                                                                                                                                                tttATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGt                                                               >  1:463358/1‑251 (MQ=255)                                                                                                                                                                                       agcactctgcACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATATCGTGGAAGCGGCGTTGATGCGCCGTGCCGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATAGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGCTGGTGCCACGGGCACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATACCCCCGGTTAACCGt                                                        >  2:190981/9‑251 (MQ=255)                                                                                                                                                                                       aGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGCTGAGCTAAAACGTg                                                                                                                                  >  2:74923/1‑177 (MQ=255)                                                                                                                                                                                       aGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTg                                                                                                                                  <  1:74923/177‑1 (MQ=255)                                                                                                                                                                                              cgcACTGGCTCCGCTGCCGGGCGACGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACAGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGgtgt                                                 >  1:191787/1‑251 (MQ=255)                                                                                                                                                                                                cACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCACTCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGAACCCGGTTCACCGTGCGGTGTTc                                               >  1:216886/1‑251 (MQ=255)                                                                                                                                                                                                        ccgcggcggggCGAAG‑TTCTCTT‑TGCCGGATCAATCCTTTCATATGTTTGCGGGGAGGCGTCGTTGTCGCGCCGTGTAGGTTGGGGGG‑TCTGGATTGCTT‑ACGACCTCCCGGGTTTTTATGAAGAATTGCCGCCTAATTTGTTTGATGAGCTAAACGGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTCCCTGGTGAAGGGTATGCACCCGGTTCATCGTGCGGTGTTCCTGACg                                         <  1:138886/241‑1 (MQ=255)                                                                                                                                                                                                         cgcgcccGGGTGAAGGTTC‑CTT‑TGCGGGTGCAGTTCTGTCACAGTATTGCGTGTAAGGGGTGTTGATGCGCCGTGC‑GGTTGGGGGGGTCTGGTGTGCTT‑ACGCCCTCCTGGGTTCTTAAGAAGAATTGCCGCCTAACTTGCTGGAGGAGCTAAACCGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGc                                      <  2:224833/246‑1 (MQ=255)                                                                                                                                                                                                           cTGCAGGGCGAAGGTTC‑ATT‑TGCCGTTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTCGATGAGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTACTTATGAAGAATTGCAAACTAACTTGCTTGATGAGCTAAACCGTGACCGCAGATGGCGCCACGGGAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCAGCATCACCGTGCGGAGTTCCTGACggccgg                                    >  1:335046/1‑247 (MQ=255)                                                                                                                                                                                                                   gcgagttttcgtt‑tGCCGTTTCTATCCTGTCACAGGATTTGGTGGAGGCGGCTTTGACGCGCCGGGCTGGTTGGGGGG‑TCTGGATTGCTT‑CCGATCCCCCGGGTTTTTATGAAGAATTCCCGCCTAACTGTTTTGAGGAGCTAAACCTTTCCCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCTCCGTGCGGTGTTCCTGACGGGCGTGATGttt                             <  1:368879/240‑3 (MQ=255)                                                                                                                                                                                                                                     ggTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGa                                           >  1:406227/1‑220 (MQ=255)                                                                                                                                                                                                                                               tGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTTTTATCTCTCCGCTCCGCTGTGGTTTATg  >  1:476898/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   GAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTC‑CTT‑TGCCGGTTCAATCCTGTCACATGACTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGG‑TCTGGATTGCTT‑ACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATG  >  NC_000913/1111963‑1112449

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