Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,965	T→A	I423N (ATC→AAC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,965	0	T	A	100.0%	62.5 / NA	18	I423N (ATC→AAC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/11);  total (7/11)

GGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTG  >  NC_000913/1756888‑1757036                                                                              |                                                                        ggCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAacac                                                                      >  2:284360/1‑81 (MQ=255)  gCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACt                                                                     >  2:284369/1‑81 (MQ=255)   cACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACtt                                                                    <  1:284338/81‑1 (MQ=255)    aCTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTc                                                                   <  2:284339/81‑1 (MQ=255)        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACt                                                               >  2:284361/1‑81 (MQ=255)          caccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACtga                                                             <  2:284340/81‑1 (MQ=255)                cGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGTCGCAGCTTGTTAAAGAGAACACTTCTACTGATGAttt                                                       <  1:284341/81‑1 (MQ=255)                 gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTc                                                      >  1:284366/1‑81 (MQ=255)                    aaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTAc                                                   >  2:284365/1‑81 (MQ=255)                     aaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTAcc                                                  >  1:284376/1‑81 (MQ=255)                         ggCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCt                                              <  2:284342/81‑1 (MQ=255)                              atcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGt                                          <  1:284343/80‑1 (MQ=255)                                caGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTa                                         >  2:284368/1‑79 (MQ=255)                                       tACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTgg                                <  1:284344/81‑1 (MQ=255)                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc                          <  1:284346/81‑1 (MQ=255)                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc                          <  1:284345/81‑1 (MQ=255)                                                        gtgCCGCAGCTTGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCt                <  1:284347/80‑1 (MQ=255)                                                                  ttGTTAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAg     >  2:284377/1‑81 (MQ=255)                                                                     ttAAAGAGAACACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACGGGCTCTGCAGAGCGGTCTGGCACACAAAGGTg  <  2:284348/81‑1 (MQ=255)                                                                              |                                                                        GGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTG  >  NC_000913/1756888‑1757036

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