Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,664	0	C	G	29.2%	48.9 / 12.8	49	A314P (GCA→CCA)	glvC	putative PTS enzyme II component GlvC
Reads supporting (aligned to +/- strand):  ref base C (20/14);  new base G (4/10);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGT  >  NC_000913/3862414‑3862906                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          gTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑c                                                                                                                                                                                                                                                           >  1:165802/1‑251 (MQ=255)           aCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AACGGTTCGG‑TAATGCCCACCAGCATGGCGG‑GTAAGG‑TT‑G‑CGGG‑GA‑TCaac                                                                                                                                                                                                                                               >  2:84171/1‑249 (MQ=255)             ccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACCTTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGGGTGAAGG‑TT‑G‑ctgaagaagca                                                                                                                                                                                                                                                 >  2:772071/1‑240 (MQ=255)             ccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCaa                                                                                                                                                                                                                                             <  2:596886/251‑1 (MQ=255)                    cGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGcc                                                                                                                                                                                                                                     <  2:54662/251‑1 (MQ=255)                                 catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑GGGG‑GA‑ACAGCAA‑Accccgcccctttaccccgt                                                                                                                                                                                                                        >  1:228715/1‑235 (MQ=255)                                  atcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑ATAAAG‑TT‑T‑CGGG‑GA‑TCCGGCA‑GCCCCCCCCTTTT‑ACCCgctt                                                                                                                                                                                                                      >  1:67876/1‑248 (MQ=255)                                    catATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑GA‑TCAGCAA‑GCCCCCCACTTTT‑Accccaattt                                                                                                                                                                                                                    >  1:159507/1‑245 (MQ=255)                                       atCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGGGGTAAAACCAAGTCAATCAGACCTCCCCCCGTGTGGCCCACCACACCAAAGAGATACATGCCGGTCGACATGGAGGCC‑CCCAGCACGGGG‑TGTACCGCAAACAGCAACGGTAAAAGGAACAGGAAGGGAAATCCC‑GGGGGTTCGG‑TAATTCCCACAAGCAGGGGGG‑TTAAGG‑TG‑G‑GGGG‑AA‑CCAGAAA‑GCCCGCCACTTTT‑ACCCGATTTT‑cc                                                                                                                                                                                                                 <  1:601103/250‑1 (MQ=255)                                                      aaacttCGGGACCCATTTTTGCGATAAACCCTAGTAAATCAGCCCCCCCCCAATTTTCCCCACCACACAAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTCCCCCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCC‑AGCGGTTCGG‑TAATCCCCACCAGCATGGGGG‑TTA‑GG‑TT‑GGCGGG‑AA‑TCAGAAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGt                                                                                                                                                                                                 <  1:244301/246‑1 (MQ=255)                                                      aaaCATCGGGTTCCATTTTTGCGATAAAACATATTAACCCAGACCCCCGCCCGTGTTCCCCACCACACAAAAGAACTACATTCCGGTCGACTTTGAGGCC‑GCCAGCACGGGG‑TTTACCGCAAACAGCAACGGTGAAAGAAAAAGGAGGGTAAATTCC‑AGCGTTTCGG‑TATTCCCCCCCAGCATGGCGT‑TTAAGG‑TG‑G‑CGGG‑AA‑CCAAAAG‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGaa                                                                                                                                                                                                   <  1:287311/248‑1 (MQ=255)                                                       aaCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑GGGG‑GA‑ACAGAAA‑ACCCGCCACTTTT‑ACCCGATTTT‑TCGGGGCGGCAGGGAAgg                                                                                                                                                                                                 >  1:729070/1‑249 (MQ=255)                                                             gggATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑AG‑ACAGCAAACCCCCCCCCTTTT‑ACC‑GATTTT‑TCGGCGCGGCAGCGaacaaaaatgc                                                                                                                                                                                          >  2:253046/1‑242 (MQ=255)                                                                  ccAGTTTTGGGGAAAAACCTAGTAAACCAGCCCCCCGCCCATGTTCCCCACCCCCCAAAAGAAATAAATTCCGGTCGACATGGAGGCC‑GCCAGCACGGCG‑TGTCCCACAAACAGCAAGGGTGAAATAAACAGGAAGGTAAATTCC‑GGCGGTTCGG‑TAATCCCCCCCAGCAGGGCGG‑TTAGGG‑TT‑G‑GGGA‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAg                                                                                                                                                                                     <  1:380798/251‑1 (MQ=255)                                                                               aaaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑AA‑TCAGCCA‑GCCCGCCCCTTTT‑TACCGCTTTT‑CCGGTGCGGCAGGGGAGTACATACCTAAAAAAAAGCCC‑ACGGc                                                                                                                                                                       >  1:585227/1‑251 (MQ=255)                                                                                aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑GTAAAG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGGGGCAGTGAAGTACATCGCTAAAGAAAAGCCC‑ACGGcc                                                                                                                                                                      >  1:759752/1‑250 (MQ=255)                                                                                              aGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGCGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑ATAAAG‑TT‑T‑GGGG‑GA‑GCAGCCA‑ACCCGCCAATTTT‑TCCCCATTTT‑CCGGAGCGGGAGGGGAGTACCACGCTAAAGaaaaaaccacggggcccaagaattttg                                                                                                                                                          >  1:362858/1‑226 (MQ=255)                                                                                                 ccTCCGCCCATGTTCCCCCCCACACCAAAGAGATACATTACGGCCGACATGGAGCCC‑GCCAGCACGGCG‑TGTCCCGCAAACAGCAAGGGTAAAATGAACAGAAAGGTAAATCCC‑AGCGGTTCGT‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGtt                                                                                                                                                    <  1:517911/251‑1 (MQ=255)                                                                                                  ccccgccattgtgccccACCACACAAAAGAGATACATTCGGGTCGACATTGGGGCC‑GCCAGCCCGGGG‑GGTCCCGAAAACAGCAACGGTAAAATGAAAAGGAAGGTAAATCCC‑AGCGGTTCGG‑TAATCCCCACCAG‑AAGGGGG‑TTTAGG‑TTTG‑GGGA‑AA‑CCAGCAG‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTa                                                                                                                                                   <  1:38413/237‑1 (MQ=255)                                                                                                       ccctttttGCCCACCCCACCAAAGAGATCAATTGCGGCCGAC‑TTGAGGCCCGCAAGCA‑GGCGGTCTCCCCCAAACAGAAACGGAAAAAAGAAAAGGAAGTTAATTTCG‑GGGGGTTGGG‑AAATTCCCACCCGTAGGGGGG‑TTGGGG‑TG‑G‑CGGA‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTg                                                                                                                                              <  2:165802/245‑1 (MQ=255)                                                                                                        ccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑GA‑ACAGCAA‑GCCCCCCCCTTTT‑ACCCGGTTTT‑TCGGGGGGGCCGGGAAGGACATCGCCAAAAAAATGCCC‑ACCGCACCCAAAAATTTTT‑AGTTACCGTGc                                                                                                                                             >  1:613149/1‑251 (MQ=255)                                                                                                              tcccccccacaccaAAGAGATACATTACGGTCGACATTGGGGCG‑GCCAGCTGGGGG‑TGTACCGCAAAAAACAAGGGAAAAATGAACAG‑AAGGGAAAATCCGAGTGTTCCGG‑TAATCCCCCCCAGCGGGGCGT‑TTAGGT‑TT‑G‑CGGG‑AA‑TCAGAAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGc                                                                                                                                        <  1:679661/244‑1 (MQ=255)                                                                                                              tGCCCCCCACACAAAAGATATACATTACGTCCGACATTGAGGCC‑GCCAGCACGGCG‑TGTCCCGAAAACAGCAACGGTGAAATGAACAGAAAGGTAAATCCC‑AGCGGTTCGG‑TAATCCCCACCAGCATGGCGG‑TTAGGG‑TG‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGc                                                                                                                                        <  1:297924/250‑1 (MQ=255)                                                                                                                                     attccgctggACATGGAGCCC‑GCCAGCACGGCG‑TGTCCCGCAAACAGAAAGGGGAAAAGGAACAGGAAGGAAAATCCC‑AGCGTTTCGG‑TAATCCCCACCAGCATGGCGT‑TGAAGT‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGa                                                                                                                <  1:253046/242‑1 (MQ=255)                                                                                                                                        gcGGTGGACGTTGGGGCC‑GCCGGCACGGCG‑TGTACCGCAAACATCAACGGAAAAATGAAGAGAAAGATAAATCCC‑GGCGGTTGGG‑TAATCCCCCCCAGCTGGGGGG‑TTAAGG‑TG‑G‑CGGA‑AA‑TCAGCAG‑CCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAttt                                                                                                             <  2:159507/250‑1 (MQ=255)                                                                                                                                              gACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑GA‑ACAGCAA‑GCCCCCCCCTTTT‑TCCCCATTTT‑TCgg                                                                                                                                                                                                               >  1:300295/1‑149 (MQ=255)                                                                                                                                                 aTTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑GA‑TCAGCCA‑ACCCCCCACTTTTTCCCCGTTTTT‑CCG‑TGGGGGAGTGAAGTACATCGCTAAAAAAATGACC‑ACGGCACC‑AAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCTTCCCGGAACAACGAGTTCAGCGGCTc                                                                                                    >  1:738030/1‑251 (MQ=255)                                                                                                                                                             cAGCACGGGG‑TGTACCGAAACCAGCAACGGTAAAATAAACGGAAAGATAAATCCC‑GGCGGTTCGT‑TAATCCCCACCAGCATGGGGG‑TTAAGG‑TG‑G‑GGGA‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAAc                                                                                         <  1:84171/251‑1 (MQ=255)                                                                                                                                                                           aCCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑AT‑G‑CGGG‑AA‑ACAGCAA‑ACCCCCCCCCTTT‑ACCCGATTTTTCCGG‑GCGGCAGTGGAGTACATCGCTAAAGAAAAGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCCAGGCCAAACCTCCCTCCGGGAACAACGAATTTAGCGGCCCGgcgccccaacatcaaccttgc                                                                              >  1:250300/1‑232 (MQ=255)                                                                                                                                                                              gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCCCTTTT‑ACCCGAATTT‑TCGGGGCGGCAGTGGAGTACAACCCCCAAAAAATTCCC‑CCCGCCCC‑ACAGAATTTTTAGTTACCGTGCAGGGGCAAACCGCCCTCCGGGGACAACGATTTCCTCGGCCCCGCGCCCACAATGAACTCTTGCAGa                                                                           >  2:440643/1‑249 (MQ=255)                                                                                                                                                                               aaaaCAGCAACGGTGAAATGAACGGGAAGTTAAATTCC‑‑GGGGTTTGGATAATCCCCCCCAGCATGGCG‑‑TTAAGGGTG‑G‑GGGG‑AA‑CCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgc                                                                        <  2:613149/250‑1 (MQ=255)                                                                                                                                                                                 aacagcaagggaaaaataacaaGGAAGGTAAATTCC‑AGGGGTTCGG‑TAATCCCCCCCAGCAGGGCGG‑TAAAGG‑TT‑G‑GGGA‑AA‑CCAGCAG‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgctg                                                                      <  1:120866/230‑1 (MQ=255)                                                                                                                                                                                  acagaaacggaaaaaTAAACAGGAAGATAAATCCC‑AGCGGTTCGG‑TAATCCCCCCCAGCATGGCGG‑TTAAGT‑TT‑G‑CGGA‑AA‑TCAGCAC‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgctgc                                                                     <  1:109937/239‑1 (MQ=255)                                                                                                                                                                                    aGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCCCTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAAACGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACCGAACTCTTGCAGAAGCTgcgc                                                                   >  1:286274/1‑251 (MQ=255)                                                                                                                                                                                      caacggaaaaaagaaaaGGAAGGTAAATCCC‑AGCGGTTCGG‑TAATCCCCACCAGCAGGGCGG‑TTAAGG‑TG‑G‑GGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGccc                                                                 <  1:215102/235‑1 (MQ=255)                                                                                                                                                                                       agggggaaaaagaacgagAAGGTAAACTCC‑CGCGGTGCGG‑AAATCCCCCCGAGCATGGCGG‑TTG‑GG‑TT‑T‑GGGG‑AACGCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCa                                                                <  2:241756/234‑1 (MQ=255)                                                                                                                                                                                        aCGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑AA‑TCAGCCA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGGAGGGAAGGACCTCGCTAAAGAAATGCCCCACGGC‑CC‑AAAAATTTTTTAGTTAACGTGCAGGGCAAAAACGCCTTCCGGGAAAAAAAATTTCAGCGGCTCGGCGCGCAAAATGGACCCTTGCAGAggctggcccccc                                                               >  1:444900/1‑244 (MQ=255)                                                                                                                                                                                              aaaTGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑ACAGCAA‑GCCCGCCACTTTT‑ACCCCATTTT‑TCGGTGCGGCAGTGAAGGACATCGCTAAAGAAAAGCCC‑ACGGGACC‑AAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAATCTTGCAAAGTCCGCGCCCAGTACATc                                                         >  1:405325/1‑251 (MQ=255)                                                                                                                                                                                                  gAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑TGGG‑AA‑TCAGGAA‑GCCCGCCACTTTT‑AACCGATTTT‑CCGGTGGGGCAGTGGAGGACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTTAGTTACCGGTCAGGGCAAAAACGCCTTCCGGGAACAACGATTTTAGCGGCCCGGCGCTCAAAATGAACTCTTGGAGATGCTGCGGCCCGTACATCTGa                                                      >  1:483928/1‑250 (MQ=255)                                                                                                                                                                                                    aCAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑AA‑TCAGCAA‑ACCCGCCACTTTT‑ACCCGATTTT‑CCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATTCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGGGCAGGGCAAAAACGCCTTCCGGGGACAAAGATTTCAGCGGCTCGGGGCTCAAAATGAACTCTTGCAGAAGGTGCGCCCAGGACCACTGAATg                                                   >  1:527382/1‑251 (MQ=255)                                                                                                                                                                                                        aaaGGTAAATCCC‑GGCGGTTCGG‑TAATCCCCACCAGCAGGGCGG‑TTAAGG‑TG‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgc                                               <  1:316632/250‑1 (MQ=255)                                                                                                                                                                                                             ggtgtgctcttcctatctg‑TAATCCCCCCCAGCAGGGCGG‑TTAAGG‑TG‑G‑GGGA‑AA‑TCAGAAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTc                                           <  1:307559/233‑1 (MQ=255)                                                                                                                                                                                                                    tccatgagtacgtaaacccccACCAGCATGGCGG‑T‑AAGG‑GT‑G‑CGGGCAA‑TCAACAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACan                                      <  2:362858/232‑2 (MQ=255)                                                                                                                                                                                                                        gcgagccgttaatcccccccggcgtGGCGG‑TTAGGG‑TG‑G‑CGAG‑AC‑ACAGCAG‑CCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGcc                                 <  2:431302/227‑1 (MQ=255)                                                                                                                                                                                                                              gg‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTTAGTTACCCGGCAGGGCAAAACCGCCTTCCGGGGACAACGATTTTAGCGGGTCGGCGCTCAAAATGAAATCTTGCAGAtgct                                                                       >  1:513669/1‑206 (MQ=255)                                                                                                                                                                                                                                     tCCCACCAGCATGGCGG‑TTAAGG‑GT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑AACCGATTTT‑TCGGTGCGGCAGGGAAGTACATCGCTAAAAAAATACCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCCTCCGGGAACAACGATTTCAGCGGCTCGGCGCCCCAAATGAACTCCTGCAGATGCTGCGCCCCGGACATCTGAATTCCGCCTTTCACAGCTGCCGGGCCAAAGATa                     >  1:526975/1‑250 (MQ=255)                                                                                                                                                                                                                                       ccaccaGCATGGCGG‑TTTAGG‑GT‑G‑GGGG‑AA‑TCAGCAA‑GCCCCCCCCTTTT‑ACCCGGTTTT‑CCGGTGCGGCAGGGAAGTACATCGCCAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAAACGCCTTCCGGGAACAACGATTTCAGCGGCTCAGCGCTCCAACCGAAACCTTTCAGAGGCTGCGCCCAGGACATCTGAATGCCGCCCTCCACCGATGCCGGGCCAAAAATaaa                   >  2:134255/1‑250 (MQ=255)                                                                                                                                                                                                                                             gCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAA‑CTGTc              >  1:572282/1‑248 (MQ=255)                                                                                                                                                                                                                                                tGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt           >  1:165823/1‑249 (MQ=255)                                                                                                                                                                                                                                                    gg‑TTGAGG‑TT‑G‑GGGG‑AA‑CCGGCAA‑TCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGa     <  1:720485/251‑1 (MQ=255)                                                                                                                                                                                                                                                        tAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgg  >  1:268660/1‑250 (MQ=255)                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          GTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCC‑GCCAGCACGGCG‑TGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCC‑AGCGGTTCGG‑TAATTCCCACCAGCATGGCGG‑TTAAGG‑TT‑G‑CGGG‑AA‑TCAGCAA‑GCCCGCCACTTTT‑ACCCGATTTT‑CCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCC‑ACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGT  >  NC_000913/3862414‑3862906

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