Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,670	T→A	29.6%	pseudogene (934/1617 nt)	glvC ←	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,670	0	T	A	29.6%	41.2 / 7.7	51	pseudogene (934/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (21/9);  new base A (10/3);  total (35/16)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCATATCCGCGTGGTTGCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862446‑3862905                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         catcatATCCGCGTGGTTGCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGTT‑TCGGG‑GAA‑C‑AGC‑AA‑ACCC‑G‑cc                                                                                                                                                                                                                                       >  2:442200/1‑238 (MQ=255)      tatCCGCGTGGTTGCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTAAGGGT‑TGGGG‑AAA‑C‑ACC‑AACGCCC‑G‑CAC‑Cttt                                                                                                                                                                                                                                 >  2:403333/1‑239 (MQ=255)         ccggttGGTTGCTG‑ACATCCGGGATCCAGTTTGGCGAAAAAACCAAGTCAATCAGCCCCCCGCCCATGTGCCCCACCACACCAAAGAAATACTTTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTCCCGCAAAC‑AGCAACGGTAAAATGAACAGGAAGGTAAATTCC‑GGCGGTTCGGTAATCCCCACCAGCATGGC‑GT‑TTAAGGTT‑GCGGG‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑cc                                                                                                                                                                                                                            <  2:923562/234‑1 (MQ=255)            cgTGGTTGCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGCGGG‑TTAAGG‑T‑GTGGG‑AAT‑C‑AGC‑CA‑CCCCAC‑CCC‑CTTTCA‑CCCg                                                                                                                                                                                                                          >  2:653615/1‑239 (MQ=255)              tGGTTGCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAAAGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGTT‑GCGGG‑GAC‑C‑AGC‑AA‑GCCC‑C‑CCC‑CCTTTT‑CCCGCtt                                                                                                                                                                                                                       >  1:1076523/1‑239 (MQ=255)                   gCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTTAGGTTGGCGGG‑AA‑‑C‑AGC‑‑A‑GCCC‑G‑CCACCTTTCC‑CCCGGTTTTCCg                                                                                                                                                                                                                  >  1:475769/1‑239 (MQ=255)                   gCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTAAGGTT‑GGGGG‑GAT‑T‑AGC‑CA‑ACCC‑G‑CCA‑CCTTTA‑CACGTttttgcc                                                                                                                                                                                                                  >  1:470756/1‑236 (MQ=255)                        aCAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTAAAGGT‑GCGGG‑GAT‑C‑AAC‑AA‑ACCC‑C‑CCA‑CTTTTA‑TCCCAtttttccgggcg                                                                                                                                                                                                             >  1:934176/1‑231 (MQ=255)                                    cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTTAAGGT‑GCGGG‑GAA‑T‑AGC‑AA‑ACCC‑GCCCA‑CTTTAA‑CCGGATTTCCGGGGGCgccgggaaaga                                                                                                                                                                                                   >  1:1233456/1‑229 (MQ=255)                                     aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGGT‑GTCGG‑AAA‑C‑AAC‑AA‑ACCC‑G‑CCA‑ACTTTA‑CCCCATTTTTCGGTGCgcgactggagta                                                                                                                                                                                                  >  2:672722/1‑227 (MQ=255)                                            cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTTAAGGT‑GCGGG‑AAA‑C‑AGC‑AA‑GCCC‑C‑CCC‑CTTTTA‑CCCGCTTTTCCGGTGCGGCAGTGAAGGCCATCGCTa                                                                                                                                                                                          >  1:817127/1‑239 (MQ=255)                                            cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGGT‑GCCGG‑GAT‑C‑AGC‑AA‑GACC‑G‑CCC‑CATTTA‑CCCCGATTTCCGGTGCGGCGGGGAAGGACATcccca                                                                                                                                                                                          >  1:1161620/1‑235 (MQ=255)                                                 aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑GAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTTCGGTGCGGCAGTGGAGTACATCGCTAAAGaa                                                                                                                                                                                     >  1:337780/1‑239 (MQ=255)                                                         tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGGT‑GGGGG‑GAA‑T‑AAC‑AA‑AGCC‑G‑CCC‑CCTTTA‑CCCGATTTTCCGGGGCGGCAGTGAAGTACAACGCTAAAAAAA‑AGCCCAc                                                                                                                                                                            >  2:654443/1‑239 (MQ=255)                                                             tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGTT‑GGGGG‑GAT‑C‑AGC‑AA‑ACCC‑C‑CCC‑CCTTTT‑CCCCATTTTCCGGTGCGGCGGGGAAGTACAACGCCAAAAAAA‑AGCCCACGGc                                                                                                                                                                         >  1:893952/1‑238 (MQ=255)                                                                                gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTTAAGGT‑GCGGG‑GAT‑A‑AGC‑AA‑ACCC‑G‑CCC‑CTTTTT‑CCCCCTTTTTCGGGGCGGCAGGGGAGGACATCGCCAAAGAAA‑TGACCACGGCACC‑AAAAATTTTTTAGTTa                                                                                                                                                     >  1:1183247/1‑238 (MQ=255)                                                                                       acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑GAT‑C‑AGC‑AA‑GCCC‑G‑CCC‑CCTTTATCCCCATTTT‑CTGTGCGGCCGGGGAGTACATCCCCAAAGAAA‑AACCCCCCGCCCC‑AAAGATTTTTTAGTTTCCCT‑Gccg                                                                                                                                            >  1:265052/1‑237 (MQ=255)                                                                                        caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTAAAGGT‑GCGGG‑AAT‑C‑AAC‑AA‑AACC‑C‑CCC‑CCTTTA‑CCCGCTTTTCCGCGGCGGCGGTGAGAGACAACGCCAAAAAAA‑AACCCACGGCCCC‑AAAAATTTTTTAGTTAcccggacag                                                                                                                                            >  2:996977/1‑232 (MQ=255)                                                                                                  aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGAT‑GCGGG‑AATGA‑AGC‑AG‑GCCC‑C‑CCC‑CTTTAC‑CCCGTTTTTCCGG‑GGGGCGGGGGAGAAAAAAGGCAAAAAAA‑AACCCAACGCACC‑AAAAATGTTTGAGATACCGC‑GC‑Agcggacaaacg                                                                                                                                 >  2:1237030/1‑228 (MQ=255)                                                                                                  aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTTTGGGT‑GCGGG‑GAC‑A‑AGC‑AA‑ACCC‑C‑CCC‑CTTTTT‑ACCGATTTTTCTGTGCGGCGAGGAAGGAACACGCTTAAAAAA‑TTGCCACCGCCAC‑ACAGATGTTTTTGTTAACGT‑GCAAGGGCAAAAcc                                                                                                                                  >  2:370681/1‑238 (MQ=255)                                                                                                   tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGGTT‑‑CGGG‑GAT‑C‑AGC‑CA‑ACCC‑C‑CCC‑CTTTTA‑CCCGGTTTTTCCGTGCGGCCGTGAAGTACAACGCTAAAAAAA‑AAACCACCGCCCC‑AAAAAATTTTTGGTTAACGT‑g                                                                                                                                               >  1:600534/1‑224 (MQ=255)                                                                                                   tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGGT‑GCGGG‑GAT‑C‑AGC‑AA‑ACCC‑C‑CCC‑CTTTTA‑CCCCGTTTTCCCGGGCGGCAGGGAAGGACAACCCCAAAAAAA‑ATGCCACGGCACCAAAAGATTTTTTAG‑TACCCT‑GC‑AAGGCAAAAcccc                                                                                                                                >  1:884846/1‑236 (MQ=255)                                                                                                     cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑GTAAAGGT‑GTGGG‑AAA‑A‑AGCAAA‑GCCC‑C‑CAC‑CTTTTA‑CCC‑AGTTTTCGGGGCGGGAGGGGAGGAACACGCCCAAGAAA‑AGACCACGCGCCC‑CAAGAATTTTTAGTTGACGG‑GC‑AGGGCAAAAcccgcc                                                                                                                              >  2:781763/1‑234 (MQ=255)                                                                                                          cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑GAT‑A‑AGG‑AA‑GCCC‑G‑CCC‑CCTTTA‑CCCCGTTTTTCGGTGCGGCCGGGAAGGAAATCGCTCAAAAAA‑AACCCCCCGGACC‑CAAGAGTTTTTAGGTACCCTGGC‑CGGGAAAAACGCCCTTCCgg                                                                                                                         >  2:223641/1‑239 (MQ=255)                                                                                                          cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGGT‑GCGGG‑GAT‑C‑AGG‑AA‑ACCC‑G‑CCC‑CCTTTT‑CCCGCTTTTTCGGTGCGGGAGTGGAGTACATCGCTAAAGAAA‑AGCCCCCGGCACC‑CAAAAATTTTTAGTTACCGT‑GC‑AGGGCAAAAACGCCCTCCggg                                                                                                                        >  1:118564/1‑239 (MQ=255)                                                                                                          cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGGT‑GCGGG‑GAT‑C‑AGC‑AA‑ACCC‑G‑CCA‑CTTTTA‑CCCCATTTTCCGGTGCGGGCGGGAAGTACATCCCTAAAGAAA‑TTCCCACCGCACC‑ACAGATTTTTTAGTTACCGT‑GC‑AGGGCAAAAACGCCTTCCgcg                                                                                                                        >  1:230200/1‑237 (MQ=255)                                                                                                             tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑GAT‑C‑AAC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTTAGTTACCGG‑GC‑AGGGCAAAACCGCCTTCCGGGgac                                                                                                                     >  1:379086/1‑239 (MQ=255)                                                                                                                aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGGT‑TCGGG‑AAT‑C‑AGA‑CA‑GACC‑G‑CCC‑CATTTT‑TCCGCTTTTTCCGTGGGGCAGGGGAGGAAACTCCCCAAGAAA‑ATCCCACCGCACC‑ACAGACTTTTTTGTTTCCCG‑GC‑CGCGCAACAACCCCCCCCCGGaacaac                                                                                                                  >  2:204889/1‑239 (MQ=255)                                                                                                                           gccAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑GAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTT‑CCCCCTTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑AGCCCACGGCACC‑CAAGATTTTTTAGGTACCCT‑TC‑AGGGCAAAACCCCCTTCCCGGAACAACGATTTctacc                                                                                                        >  1:67557/1‑234 (MQ=255)                                                                                                                                  cGGGGTGTACCCCAAACAAGCAACGGAAAAAGGAACAGGAAGGTAAATCCC‑AGCGTGTCGGTACCCCCCACCACCAGGGC‑GG‑TTAGGGTT‑GCGGA‑AAC‑C‑AGC‑AA‑CCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc                                                                                                    <  2:265052/236‑1 (MQ=255)                                                                                                                                     cttttACCGCAAAC‑AGCCACGGTAAAATGAACAGAAAGTTAAATTCC‑AGCTTTTCGTTAATCCCCCCCAGCATGGC‑GT‑TTAAGTTT‑CCGGG‑AAT‑C‑ACC‑AA‑CCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                              <  1:1179218/234‑1 (MQ=255)                                                                                                                                             gaaaaC‑AGAAACGGTAAAAGGAACAGGAAGGAAAATCCC‑AGGGGTTCGGTAATCCCCCCCAGCGTGGC‑GG‑TTAAGGTT‑GCGGG‑AAC‑C‑AGC‑AA‑CCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                     <  1:442200/237‑1 (MQ=255)                                                                                                                                              cAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGGT‑GCGGG‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CCTTTA‑CCCCCTTTTTCCGTGCGGCAGTGAAAGACACCGCTAAAAAAATTACCCAC‑GCCCC‑AAAGAATTTTTAGGTAGCCC‑GC‑AGGGCAAAACCCCCCTCCCGCGACAAAGATTTCAGCcgcgccgcgccccaaacgaa                                                                                     >  2:640172/1‑228 (MQ=255)                                                                                                                                                      cAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGGT‑TCGGG‑AAA‑C‑AGC‑AA‑GCCC‑C‑CCA‑CTTTTA‑CCCCATTTTCCGGTGCGGCAGTGAAATACATCCCTAAAGAAA‑TGCCCACGGCACC‑AAAAATTTTTTAGTTACCGT‑GC‑AGGGGAAAAACCCCCTCCGGGAACAACGATTTCAGCGGCCCCCCGCTCCAACTGAACCCTTTCa                                                                             >  2:1208360/1‑239 (MQ=255)                                                                                                                                                             aaagggAACGGGAGGGTAACTTCC‑AGTGGCGCGAGACTCACCACAGGCAGGGC‑GG‑TGAGGGTT‑GGGGT‑AAG‑A‑AGA‑AA‑GGCC‑G‑CCA‑CGTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGACCTCTTGCAGAtgct                                                                       <  2:1242715/233‑1 (MQ=255)                                                                                                                                                               aaagAACAGGAA‑GTAAATTCCGGGCGGCTCGGTAACTCCCCCCAGCATGGC‑GG‑TAAAGGTG‑GCGGG‑AAT‑C‑ACC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                    <  1:573248/236‑1 (MQ=255)                                                                                                                                                               aaTGAACAGGAAGGTAAATTCC‑AGGGGTTCGGTAATTCCCCCCAGCATGGC‑GG‑TTAAGGTT‑GGGGG‑AAA‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                    <  1:572913/239‑1 (MQ=255)                                                                                                                                                                 tGAACAGGAAGATAACTTCC‑AGGGTTCCGGTAATCCCCACCAGCTTGGC‑GG‑TTAGGTTT‑GCGGA‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                  <  1:654443/239‑1 (MQ=255)                                                                                                                                                                    aCAGGAAGATAAATCCC‑AGCGGTTCGGTAATCCCCCCAAGCATGGC‑GG‑TTAGGTTG‑GGGGG‑AAA‑C‑AGA‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                               <  1:314736/239‑1 (MQ=255)                                                                                                                                                                      aGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGAT‑TCGGG‑AAA‑T‑AGC‑AA‑GCCC‑G‑CCC‑CTTTTA‑CCCGATTTTTCGGTGCCGCAGGGaa                                                                                                                                                                                                     >  2:298351/1‑107 (MQ=255)                                                                                                                                                                       ggagagTAAATTCG‑GGCGGTTCGATAACCCCCCCCAGCAGGGG‑GG‑TAAGGGTG‑GGGGA‑AAA‑C‑AGA‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                             <  2:1234551/233‑1 (MQ=255)                                                                                                                                                                        gAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAAGTT‑GCGGG‑AAA‑C‑AGC‑AA‑GCCC‑C‑CCA‑CTTTTA‑CCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACCGCCCC‑AAAGATTTTTTAGTTACCGC‑TG‑AGGGCAAAACCGCCCTCCCGGAACAACGATTTCAGCGGCTCGGCGCCCCAACTTAACTCCTGCAGAAGCTGCCCCCAgcaa                                                            >  2:807529/1‑235 (MQ=255)                                                                                                                                                                                 atccgggCGGTCAGGTACTCCCCACCGGCATGGG‑GG‑TAAAGGTT‑GC‑GG‑AAT‑CAAAC‑AC‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATg                                                  <  1:177743/233‑1 (MQ=255)                                                                                                                                                                                      aGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTTAGGTT‑GCGGG‑AAG‑A‑AGG‑AA‑GACC‑C‑CCA‑CTTTTA‑CCCGACTTTCCGGGGGGGCCGGGAAGTAGCACGCTTAACAAA‑ATCCCCCGGGCCC‑AAAGAATTTTGTGTTTCCGT‑TC‑AGGGCCAAACGGCCCTCCGGGAACAACATTTTCAGCGGGTCGGGGTCAAAACTCAACTCTTGCCAGTGCTCCGCCCCGCA‑CCACTGACTgcacc                                              >  2:273552/1‑235 (MQ=255)                                                                                                                                                                                       gCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑TCGGG‑AAT‑C‑AGC‑AA‑ACCC‑G‑CCC‑CTTTTT‑CCCGATTTTTCGGTGCGGCCGTGAAGTACATCGCTCAAAAAA‑TGCCCCCGGCCCC‑CAAGATTTTTGAGGTACCCT‑GC‑AGGGCAAAAACGCCTTCCGGGGACAAAGATTTTCGAGGCTCCGCGCCCCAAACAACCTCCTTCCGATGCTGCGCCCCATC‑CCTCTGAAAgcccc                                              >  2:229027/1‑235 (MQ=255)                                                                                                                                                                                                       caccaGCAGGGC‑GG‑TTAAGGTT‑GGGGG‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  1:807529/239‑1 (MQ=255)                                                                                                                                                                                                         ccgtgcgggC‑GT‑TTAGGGTG‑GGGGG‑AAC‑C‑AGC‑AG‑CCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                          <  2:884846/232‑1 (MQ=255)                                                                                                                                                                                                          caGCATGGC‑GG‑TTTAGGTT‑GCGGG‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCC‑CTTTTA‑CCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCCCC‑CAAGATTTTTTAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCCGCGCTCCAACTGAACTCTTTCAGATGCCGCGCCCCGTACCACCTG‑ATGGCGCCTTCCACAACTGCCCGAccca                         >  2:1046543/1‑237 (MQ=255)                                                                                                                                                                                                                         aaGGGT‑GGGGGAAAT‑C‑AGA‑AA‑CCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc             <  1:781763/239‑1 (MQ=255)                                                                                                                                                                                                                               ggCGGG‑AAT‑G‑CGC‑AC‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAg        <  1:572222/237‑1 (MQ=255)                                                                                                                                                                                                                                      aaT‑C‑AGC‑AA‑CCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTCTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg  <  1:265365/238‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         CATCATATCCGCGTGGTTGCTGAACAT‑CGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAC‑AGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGGTAATTCCCACCAGCATGGC‑GG‑TTAAGGTT‑GCGGG‑AAT‑C‑AGC‑AA‑GCCC‑G‑CCA‑CTTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGATTTTTGAGTTACCGT‑GC‑AGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTA‑CATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862446‑3862905

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: ‑