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_0009134,117,1230GA100.0% 154.5 / NA 48R34C (CGC→TGC) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/24);  total (24/24)

TAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAT  >  NC_000913/4116991‑4117245
                                                                                                                                    |                                                                                                                          
tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                      <  2:1639244/139‑1 (MQ=255)
   cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGacac                                                                                                                   <  1:1231617/139‑1 (MQ=255)
                            ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                    >  2:373404/1‑113 (MQ=255)
                            ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                    <  1:373404/113‑1 (MQ=255)
                              cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCAt                                                                                        >  1:1273778/1‑139 (MQ=255)
                                 tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                      <  2:397208/106‑1 (MQ=255)
                                 tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                      >  1:397208/1‑106 (MQ=255)
                                 tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                                     >  2:284780/1‑139 (MQ=255)
                                             cTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                    <  1:1349321/96‑1 (MQ=255)
                                             cTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                    >  2:1349321/1‑96 (MQ=255)
                                                   gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATc                                                                                                 >  1:87656/1‑109 (MQ=255)
                                                   gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATc                                                                                                 <  2:87656/109‑1 (MQ=255)
                                                     gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                                     >  1:882946/1‑87 (MQ=255)
                                                     gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                                     <  2:882946/87‑1 (MQ=255)
                                                             ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                               <  1:1375428/85‑1 (MQ=255)
                                                             ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                               >  2:1375428/1‑85 (MQ=255)
                                                                     cccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcgc                                                 >  1:826351/1‑139 (MQ=255)
                                                                       cAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCaa                                               <  2:1273778/139‑1 (MQ=255)
                                                                                      tCGTTTTCTACCCAACCTGGTCTTGGGTAGATTTGCTCAAATTCGCACTTCGACACTCTAATAATATTGTCATCGTCATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttttt                                <  2:781204/139‑1 (MQ=255)
                                                                                         tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                             <  2:826351/139‑1 (MQ=255)
                                                                                                cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATg                        >  1:953369/1‑137 (MQ=255)
                                                                                                cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATg                        <  2:953369/137‑1 (MQ=255)
                                                                                                cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGtt                      >  1:1001875/1‑139 (MQ=255)
                                                                                                   aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                    >  1:1615875/1‑42 (MQ=255)
                                                                                                   aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                    <  2:1615875/42‑1 (MQ=255)
                                                                                                   aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                             >  1:1626070/1‑65 (MQ=255)
                                                                                                   aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                             <  2:1626070/65‑1 (MQ=255)
                                                                                                    aCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGTTAATGATATTGGCa                                                                                                   <  1:1027198/58‑1 (MQ=255)
                                                                                                    aCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                   >  2:1027198/1‑58 (MQ=255)
                                                                                                     ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                                     >  1:700666/1‑39 (MQ=255)
                                                                                                     ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGac                                                                                                                     <  2:700666/39‑1 (MQ=255)
                                                                                                     ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                   >  2:975145/1‑57 (MQ=255)
                                                                                                     ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                   <  1:975145/57‑1 (MQ=255)
                                                                                                        ggTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                                                       <  2:218303/98‑1 (MQ=255)
                                                                                                        ggTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                                                       >  1:218303/1‑98 (MQ=255)
                                                                                                          ttttGGGTAGATTTGCTCAAATTCGCACTGCGacac                                                                                                                   >  1:147576/1‑36 (MQ=255)
                                                                                                          ttttGGGTAGATTTGCTCAAATTCGCACTGCGacac                                                                                                                   <  2:147576/36‑1 (MQ=255)
                                                                                                          ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                               >  2:405603/1‑40 (MQ=255)
                                                                                                          ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                               <  1:405603/40‑1 (MQ=255)
                                                                                                               tgTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                       <  2:765984/74‑1 (MQ=255)
                                                                                                               ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                       >  1:765984/1‑75 (MQ=255)
                                                                                                               ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttttt                                <  1:529592/114‑1 (MQ=255)
                                                                                                               ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttttt                                >  2:529592/1‑114 (MQ=255)
                                                                                                               ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTccc            >  2:484823/1‑134 (MQ=255)
                                                                                                               ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTccc            <  1:484823/134‑1 (MQ=255)
                                                                                                                    aTTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCatat  >  1:471023/1‑139 (MQ=255)
                                                                                                                         cTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgcg                                                                            >  1:1369922/1‑60 (MQ=255)
                                                                                                                         cTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgcg                                                                            <  2:1369922/60‑1 (MQ=255)
                                                                                                                                    |                                                                                                                          
TAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAT  >  NC_000913/4116991‑4117245

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