Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,033	T→C	G112G (GGT→GGC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,033	0	T	C	100.0%	109.5 / NA	32	G112G (GGT→GGC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/18);  total (14/18)

AAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAA  >  NC_000913/1755900‑1756164                                                                                                                                      |                                                                                                                                    aaGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCact                                                                                                                               >  1:551042/1‑140 (MQ=255) aaGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCact                                                                                                                               >  1:266146/1‑140 (MQ=255)             tCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACctgtct                                                                                                                   <  1:107186/140‑1 (MQ=255)             tCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACctgtct                                                                                                                   <  2:1094286/140‑1 (MQ=255)             tCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACAGAcaggtca                                                                                                                   <  2:800538/140‑7 (MQ=255)                   ccATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGc                                                                                                             <  2:767945/140‑1 (MQ=255)                      tGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAac                                                                                                          >  2:386663/1‑140 (MQ=255)                                      tAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACg                                                                                          <  1:386663/140‑1 (MQ=255)                                      tAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACg                                                                                          <  2:551042/140‑1 (MQ=255)                                      tAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACg                                                                                          >  2:550218/1‑140 (MQ=255)                                        acgacgTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGAt                                                                                        >  2:191944/1‑140 (MQ=255)                                                tctcTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGAt                                                                                <  1:191944/140‑1 (MQ=255)                                                tctcTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGAt                                                                                <  1:550218/140‑1 (MQ=255)                                                tctcTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGAt                                                                                <  1:1114952/140‑1 (MQ=255)                                                                 gACCTTTACTTTCACCACTGATAAATCTTTTATCGGCACCAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGTCGATGGTCTGATCGTTATGGAAGTTACCg                                                               <  2:172040/140‑1 (MQ=255)                                                                       tACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTg                                                         >  2:646637/1‑140 (MQ=255)                                                                         cTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGaa                                                       <  1:555457/140‑1 (MQ=255)                                                                              accacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTaa                                                  >  2:401279/1‑140 (MQ=255)                                                                                   tGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAg                                             >  1:1058302/1‑140 (MQ=255)                                                                                    gATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGt                                            <  1:401279/140‑1 (MQ=255)                                                                                          tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTAt                                         <  2:1135431/137‑1 (MQ=255)                                                                                          tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTAt                                         >  1:1135431/1‑137 (MQ=255)                                                                                               tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAAc                                                 <  2:185627/124‑1 (MQ=255)                                                                                               tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAAc                                                 >  1:185627/1‑124 (MQ=255)                                                                                               tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAg                                 >  2:539303/1‑140 (MQ=255)                                                                                                    gCAACAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTg                            >  2:158106/1‑140 (MQ=255)                                                                                                      aaCAGCGAAATGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGaa                          <  1:158106/140‑1 (MQ=255)                                                                                                             aaaTGGTTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGt                   >  2:395217/1‑140 (MQ=255)                                                                                                                      gggTAACGTATGAAGGCTTCACCACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGc          <  2:1058302/139‑1 (MQ=255)                                                                                                                        gTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGa        >  1:598412/1‑140 (MQ=255)                                                                                                                        gTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGa        <  1:539303/140‑1 (MQ=255)                                                                                                                             gTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGT‑AAGTGCTGAACAACGGTGACCTGGGCGAAAACaa  <  1:512439/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    AAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAA  >  NC_000913/1755900‑1756164

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