Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,445	G→T	E250* (GAA→TAA)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,445	0	G	T	100.0%	106.2 / NA	29	E250* (GAA→TAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/9);  total (20/9)

ATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTG  >  NC_000913/1756298‑1756592                                                                                                                                                    |                                                                                                                                                    aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATaaa                                                                                                                                                  >  2:160670/1‑151 (MQ=255)     aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATcc                                                                                                                                               <  2:347694/150‑1 (MQ=255)        tCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGt                                                                                                                                           >  1:229156/1‑151 (MQ=255)        tCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGt                                                                                                                                           <  2:543798/151‑1 (MQ=255)        tCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGt                                                                                                                                           >  2:493524/1‑151 (MQ=255)        accGTGAGCACCTGAAAGCGAACGGCGTCGAAAACATCCACATCATCACCAAAATCGAAAACCAGGAAGGCGTCAACAACTTCGACGAAATCCTCGAACCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGt                                                                                                                                           <  2:80106/150‑1 (MQ=255)              aGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTagaag                                                                                                                                      <  2:121693/150‑1 (MQ=255)                  ccTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGtt                                                                                                                                 >  2:272735/1‑151 (MQ=255)                  ccTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGtt                                                                                                                                 >  2:188536/1‑151 (MQ=255)                                  ggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCagaaga                                                                                                                 >  2:381310/1‑151 (MQ=255)                                     gAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatga                                                                                                              >  1:464051/1‑151 (MQ=255)                                              cacaTCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaa                                                                                                     <  2:32967/151‑1 (MQ=255)                                                          aaaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTg                                                                                         <  2:9981/151‑1 (MQ=255)                                                                             gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATc                                                                        <  1:448795/149‑1 (MQ=255)                                                                              ccTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCAc                                                                      <  1:93469/150‑1 (MQ=255)                                                                               cTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCa                                                                       >  2:207704/1‑148 (MQ=255)                                                                                                aaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGAt                                                   <  1:381315/151‑1 (MQ=255)                                                                                                 aTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGAtt                                                  >  1:61726/1‑151 (MQ=255)                                                                                                 aTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGAtt                                                  >  1:577952/1‑151 (MQ=255)                                                                                                  tCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGAtt                                                  >  1:284388/1‑150 (MQ=255)                                                                                                        aaGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGAt                                           >  1:124808/1‑151 (MQ=255)                                                                                                           ccTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaa                                        >  2:292302/1‑151 (MQ=255)                                                                                                                    gCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACg                               >  1:368934/1‑151 (MQ=255)                                                                                                                       tcatGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGcc                             >  2:602681/1‑150 (MQ=255)                                                                                                                                              ggTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg     >  1:286448/1‑151 (MQ=255)                                                                                                                                                 gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  1:616815/1‑151 (MQ=255)                                                                                                                                                 gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  2:313631/1‑151 (MQ=255)                                                                                                                                                 gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  2:406917/1‑151 (MQ=255)                                                                                                                                                 gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATCCTGGATTCAATGATAAAAAACCCACGCCCGACTCGCGCAGAAGCGGGTGACGTTg  >  1:192023/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTG  >  NC_000913/1756298‑1756592

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