Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	766,690	T→A	16.3%	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.3%	95.6 / 3.1	57	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 (33/12);  new base A (4/5);  total (39/18)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGGCGACGACGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑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‑TCCGCGC‑t                                                                                                                                                             >  1:1511530/1‑239 (MQ=255)                                                                                    aCTATGGTCTGCGAAGCCGCTTGTCGGGCGGGTATTCCGCACACCTAGATGGCGGCGGAATAGCTGAAAAAAGCCGGACCAAAAGTG‑GGCGTAAA‑TGTTAACGTTAAAAAACAAGC‑GGCAAGCGGCATTGAAGGGGGTTAAAG‑GGGGG‑ACCACC‑T‑AAAGAGTGCCGCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGt                                                                                                                                                       <  2:1395316/239‑1 (MQ=255)                                                                                               cGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGGAACCCCCTGTATTTTTGTGGCCGCAGTCGGCCTCCAGGAGGGTGAGC‑GGTTTTATGGCATT‑CCCGCGC‑TTTCCAGGCCT‑GTTtccg                                                                                                                                           >  1:543662/1‑239 (MQ=255)                                                                                               cGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGTGGGTGAAGGCGCCCTCAAGGAGGGGGGGCGGTTTTAAT‑GCCTT‑CCCCCGC‑ttttctggcccggtggcg                                                                                                                                            >  2:1456012/1‑224 (MQ=255)                                                                                               cGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAAATC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAAACGCCTGTATTTTTTGGGCTGAAGGCGCCATCAAGGGGGGTGAGC‑GGTTTTATGGGCTT‑CCCGGGC‑TTTCAGTGCCT‑GTTgccg                                                                                                                                           >  1:47887/1‑239 (MQ=255)                                                                                                          gcccGGGGGATATTCCCCACCCCTATAGGGCTGCAGAATATCTGAAAAAAGCCGGACCAAACCGC‑GGGGTAAA‑TTTTTACTTAAAAAAACAAGG‑CGCAAACGGCTTTGAGGGGTTTTAAAC‑GGGGG‑ATCAAC‑T‑AAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCg                                                                                                                                 <  1:1321666/236‑1 (MQ=255)                                                                                                           tCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑GTAAAA‑T‑CAATAAGGCAAACCCCTGTATTTTTTCGGCTGAAGGGGCCCCCAAGGAGAGGGAGC‑GGTTTTATGGCATTCCCCGCCC‑TTT‑AGTGCCC‑GTTGCCGGGGCGAGTCgg                                                                                                                                >  2:1055061/1‑237 (MQ=255)                                                                                                                       tGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGGATTTTTTCGGCTGGAGTCGGCATCAAGGGGAGGGGAC‑GGTTTTATGGCATT‑TCCGCGCTTTTCAGGGCC‑‑GTTGCCCAGCCGGTTTGcccttgcgaaaggg                                                                                                                   >  2:1364771/1‑227 (MQ=255)                                                                                                                             cacCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑AGCAAA‑C‑CAAAAGAGCAACCGCCCGTAGTTTTGCGGCTGAAGCAGGCAACCAGGGGGGGGAGC‑GGTTTCATTGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGTGACGGTTCGCCATGCAGAAGCGTTGATcc                                                                                                             >  2:1194709/1‑239 (MQ=255)                                                                                                                              acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGGGCCGATTTGCCATGccg                                                                                                                         >  1:1147672/1‑224 (MQ=255)                                                                                                                              acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCca                                                                                                            >  1:1373290/1‑239 (MQ=255)                                                                                                                              acCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAAAAGTGCCACCCCCTGTATTTTTTCGGCTGAAGGCGCCATCAAGGAGAGTGAGC‑GGTTTTATGTCATT‑CTCGCGC‑TTTCAGTGGCT‑GTTGGCACGCCGCTTCGCCCTGCCGAAGCGTTGTTctc                                                                                                            >  2:1026576/1‑237 (MQ=255)                                                                                                                                       ggtgggaattatCGGAAAAAGGCCGGACAAAAACCC‑GGAGTAAA‑TGTTTACGTAAAAAAACAAGG‑CGCTAAGGGCATTGAGGGGCGTTAAAC‑GGCGG‑ATAAAC‑A‑AAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAgcgc                                                                                                    <  1:441141/229‑1 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGGAGTCGCCCTCCAGGAGAGTGAGC‑GTTTTAATTGCATT‑TCCGCGC‑TTTCAGTGCCC‑GTTGCCGAGGCGATTCGCCATGCCGGAGCGGTGTTCCAACCAGCCCtt                                                                                                  >  2:441141/1‑239 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGGCGGCGCCATCCATGAGAGTGGGC‑GTTTTAATTGCATT‑CCCGCGC‑TTTCAGTTCCTGGTTCCCGTCCCG‑TTTGCCATGCCGCAGCGTTGTTCCAACAAGCGCtt                                                                                                  >  2:43329/1‑239 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCtt                                                                                                  >  1:485056/1‑239 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAAAAGTGCAACCGCCTGTATTTTTTTGGCTGAAGTCGCCATCGAGGAGAGTGAGC‑GGGTTAATGGCATT‑ACCGCGC‑TTTCTGTGCCT‑GTCGCCGCGCCGGTTTTGCACGCAGGAGCGTTGATCCAAAAAgcgcct                                                                                                  >  2:957575/1‑237 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAAAAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCa                                                                                                                          >  1:1202592/1‑215 (MQ=255)                                                                                                                                        cTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATAAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCCGAAGc                                                                                                                     >  1:225466/1‑220 (MQ=255)                                                                                                                                         ggggaaaTATCGGAAAAAGGCCGGACGAAAACTC‑GGCGTAAA‑TTTTTACTAAAAAAAAAAAGG‑CGCTACCGTCATGGAGGGGCGTTAACC‑GGCGA‑ATAAAC‑T‑AAATAGTAAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 <  1:1314106/234‑1 (MQ=255)                                                                                                                                          gCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGGCGCCATCAAGGGGAGGGAGC‑GGTTTTATGGCATT‑CCCGCGC‑TTTCCGTCCCT‑GTTGCCGGGGCGGTTTGCCACGCCGAACGGTTGATTCAACAAACGCTTa                                                                                                 >  2:1186174/1‑238 (MQ=255)                                                                                                                                          gCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑T‑CAATAGTACAACCGGCTGTATTTTTTTGGCTGTAGTCGCCATCAAGGGGAGTGAGC‑GGTTTATTGGCAGT‑CCCGCGC‑GTTCAGTGCCT‑GTTGCCGAGCCGGTTCGCCATGCAGAAGCGGTGGTCCAAGAAGCGCGTa                                                                                                 >  2:266912/1‑238 (MQ=255)                                                                                                                                          gCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑A‑CAAACGGGCAAACGCCTGTATTTTTTCGGCTGAAGTAGCCACCTAAAGGCGGGAACGGTTTTAAAGGCTTC‑CCCTCCC‑TTCCAGTCGCG‑GTTCCCG‑GCCGCATCGCCATGTCCAAGCACTGCTCCAACAAGCGCTTa                                                                                                 >  2:1134261/1‑238 (MQ=255)                                                                                                                                           gcgaaagacTGAAAAAAGGCGGACACAAACCCGGGAGAAAT‑TGTTTAGGAT‑AAAAACAGGG‑GCGTAAAGGGATTGGAGGGGTGTTAAG‑GGGGG‑AACCAC‑TAAAATGGGGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 <  2:485056/230‑1 (MQ=255)                                                                                                                                           gcaaaattcgagAAAAAGGCGGACCAAAACCC‑GGAGTTAA‑TTTTTTCGTTTAAAAAAAAGGCCGCTAAGGGTTTTGAGGGTCGTTAAAA‑GG‑GG‑TACAACTT‑AAATTGTACAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTa                                                                                                 <  1:803907/227‑1 (MQ=255)                                                                                                                                           cGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑GTTAAC‑T‑CAAACGTGCAACCGCCTGTATTTTTGCGGAGGAAATCGCCCTCAAGGCACGGGAAC‑GGTTTAAAGGCAGT‑CCCCCGC‑CTTTTATACCC‑GTTGTCCAGCCGCTTTGCCATGCAGAAGAGTGGATCctaccaaccgctac                                                                                                >  2:1336829/1‑228 (MQ=255)                                                                                                                                              attagtgagAAAAGGCGGGACGAAACGGC‑GGCAAAAT‑TGTTGAGGTAAAAAAAAAAGG‑CTAAAGCGGTATGGAGGGGCGTTAAGC‑GGCGG‑AAAAAC‑A‑AAATAGCGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTc                                                                                            <  2:182843/231‑1 (MQ=255)                                                                                                                                              aaTATCTGAAAAACGCGGGACACAAA‑GA‑GGGGTAAAGTGTTTACTTTGAAAAACGAGG‑CGCTAACGGCTTTGAGGGGCTTTTAGC‑GGCGA‑ATCAAC‑T‑CAATGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTc                                                                                            <  1:290165/239‑1 (MQ=255)                                                                                                                                                     ggAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGA‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAA‑A‑CTACCAGGATATCTCCAGTATTTTTGGGGCTGAACAACCCATTAAAGAAAGGGAGA‑GTCTTTATGGTATT‑TCCCAGT‑GTTTAGGGCCT‑GTTGGCTAGGCGGCTAACACTGCAGAAGGGTTCATCCAAAAAggtcctaacccatggagg                                                                                      >  2:1133174/1‑221 (MQ=255)                                                                                                                                                        aaaaGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑A‑CAATAGTGCAACCGCCTGTATTTTTGCGGCCGAAGGCGCCCTCAAGGGGGGGGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGGGCCGATTCGCCATGCAGGAGCGTTGATCCAACCAACGCTTACCCTCAAACGTAGCg                                                                                  >  2:1356320/1‑239 (MQ=255)                                                                                                                                                           aGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATAAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGAGAGC‑GTTTTAAAGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGGCGGGCCGATTAGCCATGCAGAAGCGTTGATCCAACAAACGCTTACCCTCAAGCGTAGCGATg                                                                               >  2:1450707/1‑239 (MQ=255)                                                                                                                                                                                         tACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GGTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGGTTCGCCATGCAGAAACGGTGGTCCCACAAGCGCCTACCCCCCAGCGGAGCGGTGGGACGCGGACCGCACCGCAAGATACGCa                                                   >  2:606473/1‑239 (MQ=255)                                                                                                                                                                                         tACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GGTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAg                                                              >  1:922045/1‑228 (MQ=255)                                                                                                                                                                                         tACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCGGAAGGCGCCATCAAGGAGAGTGAGC‑GTTTTAAAGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATgag                                                                             >  1:13481/1‑213 (MQ=255)                                                                                                                                                                                              taaaaaaCAAGG‑CCGCAACCGGATTGAAGGGGCTTTAACGGGGGG‑ACAAAC‑T‑AAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAAc                                                 <  1:1055061/235‑1 (MQ=255)                                                                                                                                                                                                aaaaaCAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGGCGCCATCCAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGGTCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGCACAGCAAGATAAGCAACCGGTg                                            >  1:1045199/1‑239 (MQ=255)                                                                                                                                                                                                  aaaCACGG‑CTCTAGCGGGTTTGAAGGGCGTTAAAG‑GGTGT‑CACATC‑A‑AAAGAGAGCAGCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGGCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGaa                                          <  1:381756/239‑1 (MQ=255)                                                                                                                                                                                                    acaagggcgtgaagggcatgaaggggCGTTAACC‑GGGGA‑ACAAAC‑T‑AAAAAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGaaa                                         <  1:1356320/216‑1 (MQ=255)                                                                                                                                                                                                        gg‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CCATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑TCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGAc                                                                           >  1:838417/1‑200 (MQ=255)                                                                                                                                                                                                        gg‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCCAGCGTAGCGATGAGAAGCGTACCGTAACGCAAGATACGCAACCGGGGAAAAgag                                     >  1:1288685/1‑236 (MQ=255)                                                                                                                                                                                                                     cattgaggggttttAAC‑GGGGGAAACAAC‑G‑AAATAGCAAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAg                          <  2:782982/227‑1 (MQ=255)                                                                                                                                                                                                                                tgtttcagggtgaaTCAAC‑T‑CAAGAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATG‑CCCGCGC‑TTTCAGTGCCT‑GTTGCGGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAgg                <  1:838585/225‑1 (MQ=255)                                                                                                                                                                                                                                         cgg‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCCAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGGTTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAACGTAGCGATGGGACGCGTACCGCACAAGAAGATACGCAACCGGGGAAAAGTGACAAAACGGAGCTGGAAAAGGACCTGGTGAga     >  2:309483/1‑238 (MQ=255)                                                                                                                                                                                                                                         cgg‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGc     >  1:1092822/1‑239 (MQ=255)                                                                                                                                                                                                                                           g‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCgg   >  1:148991/1‑239 (MQ=255)                                                                                                                                                                                                                                             aTCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGGCGCCATCAAGGAGAGGGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑GTTGCCGAGCCGGTTTGCCATGCAGAAGGGGTGTTTCAACCAGCGCTTTCCCTCAAACGGAGGGCGGGGACGCGTACCGGCCAACACGAAACGCCACCGGGGAAAAATGTCCAAAAGGGGCTGAAACAAGCACAGTTGAGCGGa  >  2:282214/1‑239 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              CGGCGACGACGCCTGATGAGTTGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTC‑GGCGTAAA‑TGTTTACGTTGAAAAACAAGG‑CGCTAACGGCATTGAAGGGCGTTTAAC‑GGCGG‑ATCAAC‑T‑CAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGC‑GTTTTAATGGCATT‑CCCGCGC‑TTTCAGTGCCT‑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: ‑