Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,051	T→C	S118S (TCT→TCC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,051	0	T	C	100.0%	102.9 / NA	30	S118S (TCT→TCC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/19);  total (11/19)

CTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAA  >  NC_000913/1755972‑1756120                                                                                |                                                                      cTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTT‑CGGTAACGTATGAAGGTTTCACTACTGACctgtccgt                                                                     <  2:318226/81‑1 (MQ=255)         acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAac                                                              <  2:318227/81‑1 (MQ=255)                      ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTgg                                                    <  2:318228/78‑1 (MQ=255)                      ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGtt                                                  <  1:318232/80‑1 (MQ=255)                      ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTg                                                 <  2:318231/81‑1 (MQ=255)                      ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTg                                                 <  2:318230/81‑1 (MQ=255)                      ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTg                                                 <  2:318229/81‑1 (MQ=255)                        aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGAc                                               >  2:318275/1‑81 (MQ=255)                           ggCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGAt                                            <  1:318233/81‑1 (MQ=255)                            gCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATg                                           >  1:318264/1‑81 (MQ=255)                                cAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCt                                       >  1:318277/1‑81 (MQ=255)                                     aaaTGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCg                                  <  1:318234/81‑1 (MQ=255)                                     aaaTGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCg                                  >  2:318267/1‑81 (MQ=255)                                       aTGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGt                                <  1:318235/81‑1 (MQ=255)                                        tGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTa                               >  1:318280/1‑81 (MQ=255)                                        tGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTa                               <  2:318236/81‑1 (MQ=255)                                          gTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTAt                              <  2:318237/80‑1 (MQ=255)                                             gCGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGa                           <  2:318238/80‑1 (MQ=255)                                              cGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGaa                          <  2:318239/80‑1 (MQ=255)                                              cGGTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAg                         >  1:318282/1‑81 (MQ=255)                                               ggTAACGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGt                        >  1:318265/1‑81 (MQ=255)                                                  aaCGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAc                     >  1:318279/1‑81 (MQ=255)                                                   aCGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAcc                    <  1:318240/81‑1 (MQ=255)                                                   aCGTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAcc                    <  1:318241/81‑1 (MQ=255)                                                     gTATGAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGc                  >  1:318276/1‑81 (MQ=255)                                                         gAAGGTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCAtt              >  1:318284/1‑81 (MQ=255)                                                             gTTTCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAg          <  1:318242/81‑1 (MQ=255)                                                              tttCACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAgg         <  1:318243/81‑1 (MQ=255)                                                                 cACTACTGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTa       <  1:318244/80‑1 (MQ=255)                                                                   ctactGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAAc     >  2:318273/1‑80 (MQ=255)                                                                      ctGACCTGTCCGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACaaa  <  1:318245/80‑1 (MQ=255)                                                                                |                                                                      CTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAA  >  NC_000913/1755972‑1756120

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