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%	172.5 / NA	54	R34C (CGC→TGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/31);  total (23/31)

CCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGT  >  NC_000913/4116987‑4117237                                                                                                                                         |                                                                                                                   ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACt                                                                                                                  <  1:620656/139‑1 (MQ=255) ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACt                                                                                                                  >  1:386140/1‑139 (MQ=255)  cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTg                                                                                                                 >  1:1118286/1‑139 (MQ=255)    aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               >  2:987315/1‑139 (MQ=255)    aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               >  1:1180060/1‑139 (MQ=255)     tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              <  1:287621/139‑1 (MQ=255)     tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              <  2:219258/139‑1 (MQ=255)         tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                          <  1:739026/139‑1 (MQ=255)            aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                       <  2:97368/139‑1 (MQ=255)                         cTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCAt                                                                                          >  1:679357/1‑139 (MQ=255)                                 ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              <  2:475407/111‑1 (MQ=255)                                 ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                              >  1:475407/1‑111 (MQ=255)                                      tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                       <  1:957107/113‑1 (MQ=255)                                      tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                             >  1:493997/1‑139 (MQ=255)                                      tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                             >  1:605575/1‑139 (MQ=255)                                      tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                             >  1:1065306/1‑139 (MQ=255)                                      tCGCCAGCACTTCTACCAGCGAGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                       >  2:957107/1‑113 (MQ=255)                                        gCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGa                                                                           <  2:679357/139‑1 (MQ=255)                                                ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                                                   <  1:302729/139‑1 (MQ=255)                                                    aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                               >  1:9738/1‑139 (MQ=255)                                                    aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                               <  2:605575/139‑1 (MQ=255)                                                       aGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                            <  2:528607/91‑1 (MQ=255)                                                       aGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                            >  1:528607/1‑91 (MQ=255)                                                        gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     <  2:9285/113‑1 (MQ=255)                                                        gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                     >  1:9285/1‑113 (MQ=255)                                                         cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                          <  1:423832/139‑1 (MQ=255)                                                                  ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTgg                                                 <  2:493997/139‑1 (MQ=255)                                                                  ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTgg                                                 <  2:1065306/139‑1 (MQ=255)                                                                   tGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTg                                                                                              >  1:1251727/1‑93 (MQ=255)                                                                   tGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTg                                                                                              <  2:1251727/93‑1 (MQ=255)                                                                   tGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGt                                                >  1:53609/1‑139 (MQ=255)                                                                       tGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAg                                            <  2:9738/139‑1 (MQ=255)                                                                       tGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTAGTCGAg                                            <  2:386140/139‑1 (MQ=255)                                                                        ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                          >  1:21288/1‑76 (MQ=255)                                                                        ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                          <  2:21288/76‑1 (MQ=255)                                                                        ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc                                           >  1:911905/1‑139 (MQ=255)                                                                               aTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGa                                    >  2:442107/1‑139 (MQ=255)                                                                                         ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttt                          >  2:96483/1‑139 (MQ=255)                                                                                              tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTTGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                     <  1:689402/139‑1 (MQ=255)                                                                                              tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                     <  2:1118286/139‑1 (MQ=255)                                                                                               gtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGt                    <  2:911905/139‑1 (MQ=255)                                                                                                ttCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTc                   <  1:96483/139‑1 (MQ=255)                                                                                                          ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               <  1:567064/37‑1 (MQ=255)                                                                                                          ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                               >  2:567064/1‑37 (MQ=255)                                                                                                          ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTgg                                                                                             <  2:1264984/55‑1 (MQ=255)                                                                                                          ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTgg                                                                                             >  1:1264984/1‑55 (MQ=255)                                                                                                             ggTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTc      <  1:442107/139‑1 (MQ=255)                                                                                                               ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa           <  2:858188/132‑1 (MQ=255)                                                                                                               ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa           >  1:858188/1‑132 (MQ=255)                                                                                                                 ttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt  <  1:36587/139‑1 (MQ=255)                                                                                                                 ttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt  <  1:265653/139‑1 (MQ=255)                                                                                                                 ttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt  <  2:280228/139‑1 (MQ=255)                                                                                                                        gATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                          >  2:742428/1‑76 (MQ=255)                                                                                                                        gATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                          <  1:742428/76‑1 (MQ=255)                                                                                                                                         |                                                                                                                   CCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGT  >  NC_000913/4116987‑4117237

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