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

ATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAG  >  NC_000913/4116990‑4117255                                                                                                                                      |                                                                                                                                     aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                                                 >  2:10411/1‑139 (MQ=255) aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                                                 >  2:363667/1‑139 (MQ=255)   aGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                                               >  1:5288/1‑139 (MQ=255)   aGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                                               <  1:748708/139‑1 (MQ=255)    gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                              <  2:723372/139‑1 (MQ=255)    gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                              <  1:222286/139‑1 (MQ=255)    gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                              <  2:582054/139‑1 (MQ=255)     cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGacac                                                                                                                             <  2:5288/139‑1 (MQ=255)     cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAG‑‑CTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGc                                                                                                                           <  1:588413/139‑1 (MQ=255)      tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                                            >  1:490175/1‑139 (MQ=255)                   gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTgg                                                                                                               <  2:490175/139‑1 (MQ=255)                         atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                                         <  2:407920/123‑1 (MQ=255)                         atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                                         >  1:407920/1‑123 (MQ=255)                         atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGt                                                                                                         >  1:35773/1‑139 (MQ=255)                         atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGt                                                                                                         >  1:555690/1‑139 (MQ=255)                           atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                       <  1:635944/139‑1 (MQ=255)                              ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                                    >  1:435437/1‑139 (MQ=255)                                cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCAt                                                                                                  >  2:473776/1‑139 (MQ=255)                                   tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTg                                                                                                        <  1:729757/130‑1 (MQ=255)                                   tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTg                                                                                                        >  2:729757/1‑130 (MQ=255)                                                 aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                                 <  1:108204/139‑1 (MQ=255)                                                    aGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtgg                                                                              >  1:41247/1‑139 (MQ=255)                                                      cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCt                                                                                                                          <  1:613876/93‑1 (MQ=255)                                                      cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCt                                                                                                                          >  2:613876/1‑93 (MQ=255)                                                      cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                             >  1:386629/1‑106 (MQ=255)                                                      cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                             <  2:386629/106‑1 (MQ=255)                                                       gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                                                                           <  1:363667/139‑1 (MQ=255)                                                         ggAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGgtggtg                                                                         <  2:435437/139‑1 (MQ=255)                                                           aGCTTTGGGTGGCCCAGAATTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCAGGTGGTGcc                                                                       <  1:518019/139‑1 (MQ=255)                                                                 gggTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                                                                 <  2:593678/139‑1 (MQ=255)                                                                      gCCCAGATTTCCATTGGGTTGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcg                                                            <  2:59497/139‑1 (MQ=255)                                                                       cccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                       <  1:611265/95‑1 (MQ=255)                                                                       cccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                       >  2:611265/1‑95 (MQ=255)                                                                                     gggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                              <  1:371115/58‑1 (MQ=255)                                                                                     gggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                              >  2:371115/1‑58 (MQ=255)                                                                                           tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                                       >  1:580228/1‑139 (MQ=255)                                                                                           tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                                       <  2:46073/139‑1 (MQ=255)                                                                                           tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                                       <  1:140188/139‑1 (MQ=255)                                                                                            gtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGt                                      <  1:306272/139‑1 (MQ=255)                                                                                             ttCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTc                                     <  1:269213/139‑1 (MQ=255)                                                                                                  cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGtt                                >  1:312080/1‑139 (MQ=255)                                                                                                       ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                           <  1:473776/139‑1 (MQ=255)                                                                                                       ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                           >  1:657568/1‑139 (MQ=255)                                                                                                            ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                       >  2:451222/1‑138 (MQ=255)                                                                                                            ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                       <  1:451222/138‑1 (MQ=255)                                                                                                              ttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt                    <  2:35773/139‑1 (MQ=255)                                                                                                              ttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt                    <  1:646449/139‑1 (MQ=255)                                                                                                                gggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAg                  <  2:510009/139‑1 (MQ=255)                                                                                                                gggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAg                  <  1:10411/139‑1 (MQ=255)                                                                                                                                aTTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAg  >  1:691337/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     ATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAG  >  NC_000913/4116990‑4117255

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