Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	766,690	T→A	16.0%	L236H (CTC→CAC)	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,690	0	T	A	16.0%	88.4 / 2.5	52	L236H (CTC→CAC)	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 (31/10);  new base A (3/5);  total (36/16)
Fisher's exact test for biased strand distribution p-value = 4.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGACGACGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGA  >  NC_000913/766453‑766922                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             cGGCGACGCCGCCTGATGATTTGCTCTCGGCGACGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑t                                                                                                                                                                                                                                              <  2:516280/238‑1 (MQ=255)    cgacgaCGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑Caa                                                                                                                                                                                                                                          <  2:1042807/238‑1 (MQ=255)                          ctcGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCCCACACCTATCTGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGAGTTTAAC‑GGCGG‑ATACAC‑A‑CAATATAGCAACCCCCTGTAtttttt                                                                                                                                                                                                                   >  1:530517/1‑238 (MQ=255)                                     aggACAGGGGAGGCCCCCACCCTTCTGCCTCTTTTCCCACCGCGCCAACTATCGTCGGCGCACCGGCCTGTCCGGGGGGTATTGCACCCCCTTATAGGGCTGCGAAATATATAAAAAAGGCCGGACGAAACCGC‑GGCGAAAT‑TTTTTACGTAAAAAAACAAGG‑CGCAAACGGCATTGAAGGTCGTTAAAC‑GGGGG‑ATCACC‑A‑AAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGt                                                                                                                                                                                                         <  2:583480/236‑1 (MQ=255)                                        aaaggggAGGCACCCACCCTCCTGCCTTTTTTCCCACCGCGCCAACTATGGTCGGCTTAACGGCGGGTCGGGGGGGTATTGCACACCCCTATATGGGTGCGAAATATGTGAAAAACGCCGGACGAAACCCC‑GGCGTAAT‑TGTTTAGGTTAAAAAACAAGG‑CTCTAACGGCTTTGAGGGGCTTTAAAC‑GGCGG‑AACAAC‑A‑CAATGGGGCAACCGCCTGTATTTTTGCGGCTGAAGTCGcc                                                                                                                                                                                                     <  2:1073399/235‑1 (MQ=255)                                              gAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGGAGTCGCCATCAAg                                                                                                                                                                                               >  1:583480/1‑239 (MQ=255)                                                       acccTTCTGCCTCTTTCACCAACGCCCCAATTTGCGTGTCCATACCGGCCTCGCGGGGGGGTTTGCCCCCCACATATGTGGGGGCGAAATATCAGAAAAAGGCCGGACGCAAACGG‑GGGATAAA‑TTTTTAGGTTAAAAAACAAGG‑CGCAAACGGCATGGAGGGGTTTTTAAC‑GGGGG‑TCCAAC‑T‑AAATGGGGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAg                                                                                                                                                                                      <  2:1383485/237‑1 (MQ=255)                                                            tCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGGAACCCCCTGTATTTTTGCGGCTGAAGAAGCCATCAAGGAGAGTGAGC‑Gttt                                                                                                                                                                                >  1:506746/1‑239 (MQ=255)                                                            tCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAAAAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAAACCCCTGTATTTTTTCGGCTGAAGTCGCCCTCCAGGAGAGGGAGC‑Gttt                                                                                                                                                                                >  1:782982/1‑239 (MQ=255)                                                               gCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑GTCCAC‑C‑CACTAAGGCAACCCCCCGGATTTTTGTGGGTGTAGACCGCAACAAAGAGAGGGGGG‑GCTTTaa                                                                                                                                                                             >  2:1154028/1‑239 (MQ=255)                                                                   tttttcCCACCGCTCCAATTTCCGTCTGGGTACCGGCCTGTCCGGTGGGTTTTGCTCACCCTTATGTGGCTGCAAAATATCTAAAAAAACCGGGACGAAACCTC‑GGAGAAAA‑TTTTTACGTTAAAAAAAAGGG‑CGCTACCGGCGTTGAGGGGCGTTTAGC‑GGGGA‑ATCAAC‑T‑CAAAGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGc                                                                                                                                                                         <  2:914158/238‑1 (MQ=255)                                                                    tttttCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑AACAAA‑T‑CAATAGTTCAACCCCCTGTATTTTTGGGGCTGAAGTAGCCATCAAGGCGAGTGAGC‑GTTTTGATTGc                                                                                                                                                                         >  1:384328/1‑238 (MQ=255)                                                                              gcgcCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTGCAACCGCCTGGATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑TCCGCGCt                                                                                                                                                             >  1:1511530/1‑239 (MQ=255)                                                                                    aCTATGGTCTGCGAAGCCGCTTGTCGGGCGGGTATTCCGCACACCTAGATGGCGGCGGAATAGCTGAAAAAAGCCGGACCAAAAGTG‑GGCGTAAA‑TGTTAACGTTAAAAAACAAGC‑GGCAAGCGGCATTGAAGGGGGTTAAAG‑GGGGG‑ACCACC‑T‑AAAGAGTGCCGCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGt                                                                                                                                                       <  2:1395316/239‑1 (MQ=255)                                                                                               cGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGGAACCCCCTGTATTTTTGTGGCCGCAGTCGGCCTCCAGGAGGGTGAGC‑GGTTTTATGGCATT‑CCCGCGCTTTCCAGGCCT‑GTTtccg                                                                                                                                           >  1:543662/1‑239 (MQ=255)                                                                                               cGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAAATC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAAACGCCTGTATTTTTTGGGCTGAAGGCGCCATCAAGGGGGGTGAGC‑GGTTTTATGGGCTT‑CCCGGGCTTTCAGTGCCT‑GTTgccg                                                                                                                                           >  1:47887/1‑239 (MQ=255)                                                                                                          gcccGGGGGATATTCCCCACCCCTATAGGGCTGCAGAATATCTGAAAAAAGCCGGACCAAACCGC‑GGGGTAAA‑TTTTTACTTAAAAAAACAAGG‑CGCAAACGGCTTTGAGGGGTTTTAAAC‑GGGGG‑ATCAAC‑T‑AAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCg                                                                                                                                 <  1:1321666/236‑1 (MQ=255)                                                                                                           tCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑GTAAAA‑T‑CAATAAGGCAAACCCCTGTATTTTTTCGGCTGAAGGGGCCCCCAAGGAGAGGGAGC‑GGTTTTATGGCATTCCCCGCCCTTT‑AGTGCCC‑GTTGCCGGGGCGAGTCgg                                                                                                                                >  2:1055061/1‑237 (MQ=255)                                                                                                                             cacCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑AGCAAA‑C‑CAAAAGAGCAACCGCCCGTAGTTTTGCGGCTGAAGCAGGCAACCAGGGGGGGGAGC‑GGTTTCATTGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGTGACGGTTCGCCATGCAGAAGCGTTGATcc                                                                                                             >  2:1194709/1‑239 (MQ=255)                                                                                                                              acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGGGCCGATTTGCCATGccg                                                                                                                         >  1:1147672/1‑224 (MQ=255)                                                                                                                              acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCca                                                                                                            >  1:1373290/1‑239 (MQ=255)                                                                                                                              acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAAAAGTGCCACCCCCTGTATTTTTTCGGCTGAAGGCGCCATCAAGGAGAGTGAGC‑GGTTTTATGTCATT‑CTCGCGCTTTCAGTGGCT‑GTTGGCACGCCGCTTCGCCCTGCCGAAGCGTTGTTctc                                                                                                            >  2:1026576/1‑237 (MQ=255)                                                                                                                                       ggtgggaattatCGGAAAAAGGCCGGACAAAAACCC‑GGAGTAAA‑TGTTTACGTAAAAAAACAAGG‑CGCTAAGGGCATTGAGGGGCGTTAAAC‑GGCGG‑ATAAAC‑A‑AAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAgcgc                                                                                                    <  1:441141/229‑1 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGGAGTCGCCCTCCAGGAGAGTGAGC‑GTTTTAATTGCATT‑TCCGCGCTTTCAGTGCCC‑GTTGCCGAGGCGATTCGCCATGCCGGAGCGGTGTTCCAACCAGCCCtt                                                                                                  >  2:441141/1‑239 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGGCGGCGCCATCCATGAGAGTGGGC‑GTTTTAATTGCATT‑CCCGCGCTTTCAGTTCCTGGTTCCCGTCCCG‑TTTGCCATGCCGCAGCGTTGTTCCAACAAGCGCtt                                                                                                  >  2:43329/1‑239 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCtt                                                                                                  >  1:485056/1‑239 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAAAAGTGCAACCGCCTGTATTTTTTTGGCTGAAGTCGCCATCGAGGAGAGTGAGC‑GGGTTAATGGCATT‑ACCGCGCTTTCTGTGCCT‑GTCGCCGCGCCGGTTTTGCACGCAGGAGCGTTGATCCAAAAAgcgcct                                                                                                  >  2:957575/1‑237 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAAAAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCa                                                                                                                          >  1:1202592/1‑215 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATAAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCCGAAGc                                                                                                                     >  1:225466/1‑220 (MQ=255)                                                                                                                                         ggggaaaTATCGGAAAAAGGCCGGACGAAAACTC‑GGCGTAAA‑TTTTTACTAAAAAAAAAAAGG‑CGCTACCGTCATGGAGGGGCGTTAACC‑GGCGA‑ATAAAC‑T‑AAATAGTAAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 <  1:1314106/234‑1 (MQ=255)                                                                                                                                          gCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGGCGCCATCAAGGGGAGGGAGC‑GGTTTTATGGCATT‑CCCGCGCTTTCCGTCCCT‑GTTGCCGGGGCGGTTTGCCACGCCGAACGGTTGATTCAACAAACGCTTa                                                                                                 >  2:1186174/1‑238 (MQ=255)                                                                                                                                          gCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTACAACCGGCTGTATTTTTTTGGCTGTAGTCGCCATCAAGGGGAGTGAGC‑GGTTTATTGGCAGT‑CCCGCGCGTTCAGTGCCT‑GTTGCCGAGCCGGTTCGCCATGCAGAAGCGGTGGTCCAAGAAGCGCGTa                                                                                                 >  2:266912/1‑238 (MQ=255)                                                                                                                                          gCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑A‑CAAACGGGCAAACGCCTGTATTTTTTCGGCTGAAGTAGCCACCTAAAGGCGGGAACGGTTTTAAAGGCTTC‑CCCTCCCTTCCAGTCGCG‑GTTCCCG‑GCCGCATCGCCATGTCCAAGCACTGCTCCAACAAGCGCTTa                                                                                                 >  2:1134261/1‑238 (MQ=255)                                                                                                                                           gcgaaagacTGAAAAAAGGCGGACACAAACCCGGGAGAAAT‑TGTTTAGGAT‑AAAAACAGGG‑GCGTAAAGGGATTGGAGGGGTGTTAAG‑GGGGG‑AACCAC‑TAAAATGGGGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 <  2:485056/230‑1 (MQ=255)                                                                                                                                           gcaaaattcgagAAAAAGGCGGACCAAAACCC‑GGAGTTAA‑TTTTTTCGTTTAAAAAAAAGGCCGCTAAGGGTTTTGAGGGTCGTTAAAA‑GG‑GG‑TACAACTT‑AAATTGTACAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 <  1:803907/227‑1 (MQ=255)                                                                                                                                           cGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑GTTAAC‑T‑CAAACGTGCAACCGCCTGTATTTTTGCGGAGGAAATCGCCCTCAAGGCACGGGAAC‑GGTTTAAAGGCAGT‑CCCCCGCCTTTTATACCC‑GTTGTCCAGCCGCTTTGCCATGCAGAAGAGTGGATCctaccaaccgctac                                                                                                >  2:1336829/1‑228 (MQ=255)                                                                                                                                              attagtgagAAAAGGCGGGACGAAACGGC‑GGCAAAAT‑TGTTGAGGTAAAAAAAAAAGG‑CTAAAGCGGTATGGAGGGGCGTTAAGC‑GGCGG‑AAAAAC‑A‑AAATAGCGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTc                                                                                            <  2:182843/231‑1 (MQ=255)                                                                                                                                              aaTATCTGAAAAACGCGGGACACAAA‑GA‑GGGGTAAAGTGTTTACTTTGAAAAACGAGG‑CGCTAACGGCTTTGAGGGGCTTTTAGC‑GGCGA‑ATCAAC‑T‑CAATGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTc                                                                                            <  1:290165/239‑1 (MQ=255)                                                                                                                                                        aaaaGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑A‑CAATAGTGCAACCGCCTGTATTTTTGCGGCCGAAGGCGCCCTCAAGGGGGGGGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGGGCCGATTCGCCATGCAGGAGCGTTGATCCAACCAACGCTTACCCTCAAACGTAGCg                                                                                  >  2:1356320/1‑239 (MQ=255)                                                                                                                                                           aGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATAAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGAGAGC‑GTTTTAAAGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGGCGGGCCGATTAGCCATGCAGAAGCGTTGATCCAACAAACGCTTACCCTCAAGCGTAGCGATg                                                                               >  2:1450707/1‑239 (MQ=255)                                                                                                                                                                                         tACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GGTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGGTTCGCCATGCAGAAACGGTGGTCCCACAAGCGCCTACCCCCCAGCGGAGCGGTGGGACGCGGACCGCACCGCAAGATACGCa                                                   >  2:606473/1‑239 (MQ=255)                                                                                                                                                                                         tACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GGTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAg                                                              >  1:922045/1‑228 (MQ=255)                                                                                                                                                                                         tACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCGGAAGGCGCCATCAAGGAGAGTGAGC‑GTTTTAAAGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATgag                                                                             >  1:13481/1‑213 (MQ=255)                                                                                                                                                                                              taaaaaaCAAGG‑CCGCAACCGGATTGAAGGGGCTTTAACGGGGGG‑ACAAAC‑T‑AAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAAc                                                 <  1:1055061/235‑1 (MQ=255)                                                                                                                                                                                                aaaaaCAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGGCGCCATCCAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGGTCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGCACAGCAAGATAAGCAACCGGTg                                            >  1:1045199/1‑239 (MQ=255)                                                                                                                                                                                                  aaaCACGG‑CTCTAGCGGGTTTGAAGGGCGTTAAAG‑GGTGT‑CACATC‑A‑AAAGAGAGCAGCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGGCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGaa                                          <  1:381756/239‑1 (MQ=255)                                                                                                                                                                                                        gg‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CCATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑TCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGAc                                                                           >  1:838417/1‑200 (MQ=255)                                                                                                                                                                                                        gg‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCCAGCGTAGCGATGAGAAGCGTACCGTAACGCAAGATACGCAACCGGGGAAAAgag                                     >  1:1288685/1‑236 (MQ=255)                                                                                                                                                                                                                     cattgaggggttttAAC‑GGGGGAAACAAC‑G‑AAATAGCAAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAg                          <  2:782982/227‑1 (MQ=255)                                                                                                                                                                                                                                         cgg‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCCAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGGTTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAACGTAGCGATGGGACGCGTACCGCACAAGAAGATACGCAACCGGGGAAAAGTGACAAAACGGAGCTGGAAAAGGACCTGGTGAga     >  2:309483/1‑238 (MQ=255)                                                                                                                                                                                                                                         cgg‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGc     >  1:1092822/1‑239 (MQ=255)                                                                                                                                                                                                                                           g‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCgg   >  1:148991/1‑239 (MQ=255)                                                                                                                                                                                                                                             aTCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGGCGCCATCAAGGAGAGGGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGGTTTGCCATGCAGAAGGGGTGTTTCAACCAGCGCTTTCCCTCAAACGGAGGGCGGGGACGCGTACCGGCCAACACGAAACGCCACCGGGGAAAAATGTCCAAAAGGGGCTGAAACAAGCACAGTTGAGCGGa  >  2:282214/1‑239 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             CGGCGACGACGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGCTTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGA  >  NC_000913/766453‑766922

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