Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,123	G→A	R34C (CGC→TGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,123	0	G	A	100.0%	123.2 / NA	39	R34C (CGC→TGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/24);  total (15/24)

CCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGT  >  NC_000913/4116987‑4117237                                                                                                                                         |                                                                                                                   ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCCACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACt                                                                                                                  >  2:399993/1‑139 (MQ=255) ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACt                                                                                                                  >  2:551144/1‑139 (MQ=255)    aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               >  1:110933/1‑139 (MQ=255)    aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               >  2:215424/1‑139 (MQ=255)    aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               >  1:517858/1‑139 (MQ=255)     tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              <  1:215424/139‑1 (MQ=255)     tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              <  2:517858/139‑1 (MQ=255)      aGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                             <  1:100419/139‑1 (MQ=255)       gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                            <  1:444449/139‑1 (MQ=255)         tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                          <  2:294323/139‑1 (MQ=255)          gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGc                                                                                                         <  2:781028/139‑1 (MQ=255)            aaTGTGATCGGGACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                       <  1:399993/139‑1 (MQ=255)                   tCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGatat                                                                                                <  2:648153/139‑1 (MQ=255)                              atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     <  2:723295/139‑1 (MQ=255)                                 ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                  >  1:575025/1‑139 (MQ=255)                                         ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                           >  1:756565/1‑122 (MQ=255)                                         ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                           <  2:756565/122‑1 (MQ=255)                                         ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGAc                                                                          <  1:795337/139‑1 (MQ=255)                                                        gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATAttc                                                                                              >  1:560916/1‑103 (MQ=255)                                                        gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTg                                                                                              <  2:560916/104‑1 (MQ=255)                                                        gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggt                                                           <  2:110933/139‑1 (MQ=255)                                                         cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     >  1:686458/1‑112 (MQ=255)                                                         cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     <  2:686458/112‑1 (MQ=255)                                                                  ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTgg                                                 <  1:402992/139‑1 (MQ=255)                                                                       tGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAg                                            <  2:575025/139‑1 (MQ=255)                                                                              gATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACg                                     <  2:632612/139‑1 (MQ=255)                                                                                        gggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                               >  1:767787/1‑103 (MQ=255)                                                                                        gggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                               <  2:767787/103‑1 (MQ=255)                                                                                            cGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              <  2:852735/52‑1 (MQ=255)                                                                                            cGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              >  1:852735/1‑52 (MQ=255)                                                                                                ttCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGAGACAGGCTAATGATATTGGCATCGTGATCCATTACGGCCGCGCGGGAGCTGGTGTTGCCCTGGTCGAGCGGAACGATATATTTTTTTTCAGTc                   <  2:837434/139‑1 (MQ=255)                                                                                                             ggTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                       >  2:596384/1‑42 (MQ=255)                                                                                                             ggTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                       <  1:596384/42‑1 (MQ=255)                                                                                                                         aTTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt  >  2:183221/1‑131 (MQ=255)                                                                                                                         aTTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt  <  1:183221/131‑1 (MQ=255)                                                                                                                             gCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     <  2:247358/44‑1 (MQ=255)                                                                                                                             gCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     >  1:247358/1‑44 (MQ=255)                                                                                                                                   aTTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                               >  1:768773/1‑60 (MQ=255)                                                                                                                                   aTTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                               <  2:768773/60‑1 (MQ=255)                                                                                                                                         |                                                                                                                   CCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGT  >  NC_000913/4116987‑4117237

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