Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,359	T→C	I221T (ATC→ACC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,359	0	T	C	100.0%	76.4 / NA	21	I221T (ATC→ACC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/11);  total (10/11)

TCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATG  >  NC_000913/1756271‑1756486                                                                                         |                                                                                                                                tCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAAcc                                                                                                                         <  2:573578/97‑1 (MQ=255) tCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAAcc                                                                                                                         >  1:573578/1‑97 (MQ=255)                cGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAATAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTg                                                                >  1:377391/1‑139 (MQ=255)                     tGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatc                                                                       >  1:160238/1‑127 (MQ=255)                     tGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatc                                                                       <  2:160238/127‑1 (MQ=255)                                       tGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCg                                                                                        <  1:399081/92‑1 (MQ=255)                                       tGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCg                                                                                        >  2:399081/1‑92 (MQ=255)                                       tGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaa                                         >  1:476915/1‑139 (MQ=255)                                           cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATcc                                     >  2:420954/1‑139 (MQ=255)                                            aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATccc                                    <  1:572280/139‑1 (MQ=255)                                              cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGAc                                                     <  1:2744/120‑1 (MQ=255)                                              cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGAc                                                     >  2:2744/1‑120 (MQ=255)                                                 aaaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTag                               <  2:476915/139‑1 (MQ=255)                                                   aGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTag                               >  2:92876/1‑137 (MQ=255)                                                   aGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTag                               <  1:92876/137‑1 (MQ=255)                                                            cggcgAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATcc                                     <  1:262406/122‑1 (MQ=255)                                                            cggcgAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATcc                                     >  2:262406/1‑122 (MQ=255)                                                               cgAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcat                                                                        <  1:579099/84‑1 (MQ=255)                                                               cgAAAACATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcat                                                                        >  2:579099/1‑84 (MQ=255)                                                                  aaaCATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGcc              <  1:458455/139‑1 (MQ=255)                                                                  aaaCATCCACATCATCTCCAAAACCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGcc              <  1:420954/139‑1 (MQ=255)                                                                              catcGCCAAAACCGATAACCAGGACGGCGTCAACAACTTCGGTGAAATCCTCGAAGTCTCGGGCGTGTTCATGGTTGCGCGTGGCGAGCTTCGTGGATATATCCCGGTAGTAGTAGTTACCTTTGCCTAGAAgatgatg  <  1:323574/139‑1 (MQ=255)                                                                                         |                                                                                                                                TCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATG  >  NC_000913/1756271‑1756486

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