Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	766,694	T→A	16.2%	N237K (AAT→AAA)	mngA →	fused 2‑O‑a‑mannosyl‑D‑glycerate specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	766,694	0	T	A	16.2%	79.3 / 2.4	52	N237K (AAT→AAA)	mngA	fused 2‑O‑a‑mannosyl‑D‑glycerate specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base T (34/7);  new base A (4/4);  total (39/13)
Fisher's exact test for biased strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGA  >  NC_000913/766462‑766922                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         cgCCTGATGAGTTGCTCTCGGCGCTGGATGCCAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTAATGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCCGACGCAAACTCGGCGTAAATG‑TTTCCGTTGAACAACAATGCGCTCACGGCATTGACGGGCGTTTAA‑CGGCGG‑GTCAG‑CTC‑CACAGTACa                                                                                                                                                                                                                                    >  2:455882/1‑239 (MQ=255)            ttGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAAGGCAA‑CCGCCTGTA‑t                                                                                                                                                                                                                       >  1:66418/1‑239 (MQ=255)                 ctcGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCCGTA‑tttttt                                                                                                                                                                                                                  >  1:425947/1‑238 (MQ=255)                          ggATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTTCGGCTGAAg                                                                                                                                                                                                         >  1:722663/1‑239 (MQ=255)                          ggATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCCCCTGGA‑TTTTTGCGGCTGAAg                                                                                                                                                                                                         >  1:196573/1‑239 (MQ=255)                                     gAGGCACGCAACCTTCTGCATCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCCA‑CTC‑CATAATGCCA‑CCGCCTGTA‑TTTTTTCGGCTGGAGTCGCCATCaaa                                                                                                                                                                                              >  1:416112/1‑238 (MQ=255)                                         cACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AAAAGTGCAA‑CCGCCTGTA‑TTTTTTCGGCTGAAGTCGCCATCAAGgaga                                                                                                                                                                                          >  1:1150146/1‑239 (MQ=255)                                           cGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAATGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAAGAGA‑Gt                                                                                                                                                                                       >  1:1000855/1‑239 (MQ=255)                                              acccTTCTCCCTTTTTTCCCAACGCCCCAACTCTCGTCGGCGAAGCGGCCTGTCGGGGGGTTTTTGCTCACCCTTATAGGGCGGCGGAATATCGGAAAAAGGCCGGACGAAAACGGGGCGTAATTG‑TTTACTTAAAAAAAAAGGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGGGG‑ACCAA‑CTA‑AATGGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAg                                                                                                                                                                                    <  2:1150146/237‑1 (MQ=255)                                                          ctTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGCTTAA‑CCGCCTGTA‑TTTTTGCGGCTCAAGTCGCCATCTAGGAGA‑GAGAGAGTTTTAATgg                                                                                                                                                                         >  1:1414563/1‑238 (MQ=255)                                                                    cgcgCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑GTTAA‑ATCAAAAAAGGCAA‑CCGCC‑GTATTTTTTTGGGCTGAGGTCCCCA‑CAAAGAGAGGGGGGCGTTTTAAGGCCTTTc‑cccg                                                                                                                                                               >  2:1305034/1‑236 (MQ=255)                                                                                      cGTAACGGCCTGTCCGGCGGGGATTGCTCACACCTATATGGATGCGGAATATCTGGAAAAAGACGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑AGCAA‑ATC‑AAATATGAAA‑CCCGCCGGA‑TTTTTGGGGCGGGAGTCGGCATCAAGGGGA‑GGGGGGGGTTTAATGGGAATCACCGCGCTTTCatgccggtgtccg                                                                                                                                            >  2:87635/1‑227 (MQ=255)                                                                                                tGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑GGCAA‑CTC‑AATAGTGCAA‑CCCCCGGTA‑TTTTTTTGGCTGGAGGCGCCCTCAATGGGG‑GGGAGGGGGTTGATGGCATTA‑CCGCGCTTTTAGTTCCTGTTtcggagacgggtt                                                                                                                                  >  2:574530/1‑234 (MQ=255)                                                                                                tGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGCTTAA‑CGGCGG‑CTCCA‑ACC‑AATAGTGCCA‑CCGCCCGTA‑TTTTTTTTGCCGTAGTCGCCATCAATGAGG‑GGGGGCGCTTTAATTGCATTT‑CCCAGCTTTTCCGGTCTGTTGTCTAGCCCATTc                                                                                                                                  >  2:620504/1‑238 (MQ=255)                                                                                                ggccgggggggTATTGCCCACACTAAAGTGCCGGGGAAATATCGGAAAAAAGCGGGACGAAACCGCGGTGTAAATT‑TTTACGTAAAAAAACAGGGCGCTACCGGCA‑TGGAGGGGGTTTAAGGGGGGG‑AAAAA‑CAC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCg                                                                                                                                 <  1:721439/231‑1 (MQ=255)                                                                                                   cgggggggTATTGCCCACCCATATATGGCGGCAGAATATCAGAAAAAGGCCGGACGAAACCTGGGGAAAAATG‑TTTACGTTAAAAAACAGGCCGCAAGCGGTTTGGAGGGGTTTTGAA‑CGGGGA‑ACAAA‑CTC‑AATAGTGCAC‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGcc                                                                                                                               <  1:872955/233‑1 (MQ=255)                                                                                                         gggtttacaccaaacTATATGGCGGCGAAATATGTGAAAAAGGCCGGGCGCAACCCCGGGGTAAATG‑TGTTCGTTAAAGACCGAGGCACAACGGGCATGGAGGGGCGTTTAA‑GGGGGG‑TAAAC‑CTA‑AATGGCGCAC‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAg                                                                                                                         <  2:66418/224‑1 (MQ=255)                                                                                                                     acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGG‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATTGAGAAACGTTTATCca                                                                                                            >  1:1067208/1‑239 (MQ=255)                                                                                                                     acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCca                                                                                                            >  1:581798/1‑239 (MQ=255)                                                                                                                          tatGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCcaacaa                                                                                                        >  1:1092097/1‑238 (MQ=255)                                                                                                                          tatGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑CTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATGAA‑CTC‑AAAAGTTCAA‑CCCCCGGTA‑TTTTTTAGGCTGAAGTTGCCATCAAGGAGA‑GGAAGACTTTTAAATGGATTC‑CCGCGGTTTCAGAGCCGGTTGCCTCGGCGGTTCGCCTTTGAGAAGCGTTAATTcaaaaa                                                                                                        >  2:1436157/1‑238 (MQ=255)                                                                                                                               cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATTAA‑CCC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGGAGTCGCCATCAAGGAGG‑GGGGGCGTTTTAATTGCATTC‑CCGCGCCTTCAGTGGCTGTGGCCGGGCCGATTCGGCAATCAGAAGCGTTGCTcacaccagacct                                                                                                   >  2:184173/1‑227 (MQ=255)                                                                                                                               cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGGAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAACGCtt                                                                                                  >  1:1187434/1‑239 (MQ=255)                                                                                                                               cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGG‑AGGAGCGGTTTAATGGCATTT‑CCCCGCTTTCAGTGCCTGTTGCCGGGCCGAGTCGCCATGCAGAAGCGTTGGTCCCACAAACGCt                                                                                                   >  2:104489/1‑238 (MQ=255)                                                                                                                               cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGGGCCGATTCGCCCTGCAGAAGCGTTGATCCAACAAGCGCtt                                                                                                  >  1:1432887/1‑239 (MQ=255)                                                                                                                                tGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑ATC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATTGCATTC‑CCGCGCTTTCAGTGCATGTTGCCGAGCCGATTCGCCATGCAGAAACGTTGATCCAACAAGAGCTTa                                                                                                 >  1:580958/1‑239 (MQ=255)                                                                                                                                 gCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 >  1:418599/1‑238 (MQ=255)                                                                                                                                 gCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 >  1:1055168/1‑238 (MQ=255)                                                                                                                                 gCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GGGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 >  2:1444316/1‑238 (MQ=255)                                                                                                                                 gCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGGCGCCCTCAAGGAGA‑GGGAGGGGTTTAATGGCATTC‑CCGCGCTTTCAGTGGCTGTTGCCGAACCGATTCGCCCTTCAGAAGCGGTGATCCAACAAGCGCTTa                                                                                                 >  2:14681/1‑238 (MQ=255)                                                                                                                                  cGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑CATAGGGCAA‑CCGCCTGTA‑TTTTTTCGGCTTAAGGCGCCATCAAGGGGA‑GAGAGCGTTTTAATGGAACTT‑CCGCGCGTTTCGTGCCTGGGGCCGAGACGCTTCGCCCTGCAGACGCGTTGTTCCAACAAGCGCCTAc                                                                                                >  2:1187955/1‑238 (MQ=255)                                                                                                                                  cGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑CATAGAGCAA‑CCGCCTGTT‑TTTTTTCGGCTGAAGTCGCCCTCCAGGAGA‑GTGAGCGTTTTTATGGCATTT‑CCGCGCTTTCAGTGCCTTTTTCCGAGCCCATTCGCCCTGCAGAAGCCTTGGTCCAACAAACGCTTAc                                                                                                >  2:1380144/1‑238 (MQ=255)                                                                                                                                  cGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTTCGGCTGAAGGCGCCATCAAGGAGA‑GGGAGGGGTTTAATGGCCTTC‑CCGCGCTTTTAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAACGTTGTTCCAAAAAGGGCTTAc                                                                                                >  2:1245289/1‑238 (MQ=255)                                                                                                                                  cGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAAGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGCTTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTAc                                                                                                >  1:721160/1‑238 (MQ=255)                                                                                                                                            ggAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGGGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCCCAAGCGTa                                                                                     >  1:996199/1‑239 (MQ=255)                                                                                                                                                 aggCCGGACGAAAAC‑GGGGGTAAATGTTTTACGTTAAAAAAAAGGGCGCTAACGGAATTGAGGGGCG‑TTAA‑CGGGGGAAACAA‑CAA‑AATGTGGAAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGAt                                                                                <  1:567515/237‑1 (MQ=255)                                                                                                                                                   gCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGga                                                                                                                                                                                            >  1:596877/1‑131 (MQ=255)                                                                                                                                                   gCCGGACGCAAACTCGGCGTAAAAG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑ATC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GGGAGGGGTTTAATTGCATTC‑CCCCGCTTTCCGTGCCTTTTGCCCAGCCGATTCTCCATGCAGAAGCGGTGATCCAACAAGCGCTTACCCCCAAGCGGAGCGATga                                                                              >  2:910696/1‑239 (MQ=255)                                                                                                                                                        aCGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGAGCAA‑ACGCCTGTA‑TTTTTGCTGCTGAAAGCGGCGACAAGGAGA‑GGGGGCGGTTTAATGGGCTTC‑CCGCGCTTTGAGGGACATTTGCCCGGCCGATTCCCCATGCCGAACCGTTGATCCACCCAACGCTTACCCTCACGCGTAGCggttgcagt                                                                          >  2:732682/1‑230 (MQ=255)                                                                                                                                                                    cGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGc                                                             >  1:1044371/1‑239 (MQ=255)                                                                                                                                                                     gcgaaatgttttacgtaaaaaaaaaGGGGCAAAGGGGATGTAAGGGGGTTTAA‑AGGGGG‑AACAACCAA‑AAGAGTGCAT‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACa                                                               <  2:416112/221‑1 (MQ=255)                                                                                                                                                                      aaaaTG‑TTTACGTTAAAAAAAAGGGCCATAACGGAATTGAGGGGGTTTAAA‑GGGGGG‑ATCAA‑CTC‑AAAGGTGAAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCa                                                            <  1:104489/237‑1 (MQ=255)                                                                                                                                                                               tACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTT‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAACCGCTTACCCTCAAGCGTAGCGATGAGAAGCGGACCGTAAAGCAAGATACGCa                                                   >  2:712327/1‑239 (MQ=255)                                                                                                                                                                               tACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCa                                                   >  1:737198/1‑239 (MQ=255)                                                                                                                                                                                                    aacgtaaTTGAGGGCCTTTTAA‑CGGGGG‑AACAA‑CCA‑AATAGTCCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAAc                              <  1:87635/233‑1 (MQ=255)                                                                                                                                                                                                          cttGAAGGGGGTTTAA‑CGGGGGAAACAC‑CTA‑AAT‑GTGGAATCCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAg                          <  1:736115/237‑1 (MQ=255)                                                                                                                                                                                                           tgagggggttttaaa‑CGGGGA‑AACAA‑CTA‑AAGGGCGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGTAACCGGTGAAAAGTGTCAAAACGGAGCt                        <  2:1055168/225‑1 (MQ=255)                                                                                                                                                                                                               aGGGCGTTTAA‑CGGGGC‑ATCAA‑CTA‑AAAAGATCAACCCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGaa                     <  1:1170467/238‑1 (MQ=255)                                                                                                                                                                                                                     tttAA‑CGGCGG‑ATCAA‑CTC‑CATAGTGCAA‑ACGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GGGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCCACCAACGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATAACCAACCGGTGAAAAGTGTCAAAACGGAGCTGCAACAggaa              >  2:190861/1‑236 (MQ=255)                                                                                                                                                                                                                            ggcgg‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCCTCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCCGAAGCGGTGATCCAACAAGCGCTTACCCCCAAGCGTAGCGATGAGACGCGTACCGTACAGCCAGATACGCCAACGGGGAAAAGTGGCAAAACGGAGCTGAAACAGGccccggtga       >  2:483749/1‑231 (MQ=255)                                                                                                                                                                                                                                  aacAA‑CAC‑AAAAGTCAAC‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGa  <  1:802376/237‑1 (MQ=255)                                                                                                                                                                                                                                   atca‑CTA‑AAAAGTGCCACCCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGa  <  1:910696/236‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         CGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATG‑TTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAA‑CGGCGG‑ATCAA‑CTC‑AATAGTGCAA‑CCGCCTGTA‑TTTTTGCGGCTGAAGTCGCCATCAAGGAGA‑GTGAGCGTTTTAATGGCATTC‑CCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGA  >  NC_000913/766462‑766922

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