Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,248,380 A→C 100% D295E (GAT→GAG dhaM ← fused predicted dihydroxyacetone‑specific PTS enzymes: HPr component/EI component

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,248,3800A→C100.0% 112.9 32885dhaMfused predicted dihydroxyacetone‑specific PTS enzymes: HPr component/EI component
Reads supporting (aligned to +/- strand):  new base (19/13):  ref base (0/0):  total (19/13)

CTGGCTTCTGCTTTCGCTGTTAACGTCATCAGATCTAATAACGTGAAGTCAATAGCCTGGCGTA  >  REL606/1248348‑1248411
                                |                               
cTGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTCt                              <  1:4779847/36‑1 (MQ=255)
 tGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTCta                             >  1:4840322/1‑36 (MQ=255)
  ggCTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaa                            >  1:4108887/1‑36 (MQ=255)
   gCTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaat                           <  1:4123774/36‑1 (MQ=255)
   gCTTCTGCTTTCGCTATTAACGTCATCAGCTCtaat                           <  1:3066381/36‑1 (MQ=37)
    cTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaata                          <  1:1473983/36‑1 (MQ=255)
    cTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaata                          >  1:5197882/1‑36 (MQ=255)
     ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                         <  1:1354107/36‑1 (MQ=255)
     ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                         >  1:4741444/1‑36 (MQ=255)
         gCTTTCGCTGTTAACGTCATCAGCTCTAATAACGTg                     >  1:1927389/1‑36 (MQ=255)
           tttCGCTGTTAACGTCATCAGCTCTAATAACGTGaa                   <  1:1350780/36‑1 (MQ=255)
           tttCGCTGTTAACGTCATCAGCTCTAATAACGTGaa                   >  1:1657580/1‑36 (MQ=255)
            ttCGCTGTTAACGTCATCAGCTCTAATAACGTGAAg                  >  1:873013/1‑36 (MQ=255)
            ttCGCTGTTAACGTCATCAGCTCTAATAACGTGAAg                  <  1:2874057/36‑1 (MQ=255)
            ttCGCTGTTAACGTCATCAGCTCTAATAACGTGAAg                  >  1:2860951/1‑36 (MQ=255)
             tCGCTGTTAACGTCATCAGCTCTAATAACGTGAAGt                 >  1:2207196/1‑36 (MQ=255)
             tCGCTGTTAACGTCATCAGCTCTAATAACGTGAAGt                 >  1:2092315/1‑36 (MQ=255)
              cGCTGTTAACGTCATCAGCTCTAATAACGTGAAGTc                >  1:2674855/1‑36 (MQ=255)
               gCTGTTAACGTCATCAGCTCTAATAACGTGAAGTCa               <  1:2668382/36‑1 (MQ=255)
               gCTGTTAACGTCATCAGCTCTAATAACGTGAAGTCa               <  1:4129153/36‑1 (MQ=255)
                 tGTTAACGTCATCAGCTCTAATAACGTGAAGTCAAt             >  1:5826995/1‑36 (MQ=255)
                 tGTTAACGTCATCAGCTCTAATAACGTGAAGTCAAt             <  1:987076/36‑1 (MQ=255)
                  gTTAACGTCATCATCTCTAATAACGTGAAGTCAATa            <  1:873807/36‑1 (MQ=37)
                  gTTAACGTCATCAGCTCTAATAACGTGAAGTCAATa            >  1:2819112/1‑36 (MQ=255)
                     aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc         >  1:4299216/1‑36 (MQ=255)
                     aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc         >  1:3512028/1‑36 (MQ=255)
                      aCGTCATCAGCTCTAATAACGTGAAGTCAATAGCCt        <  1:5709303/36‑1 (MQ=255)
                      aCGTCATCAGCTCTAATAACGTGAAGTCAATAGCCt        >  1:5752041/1‑36 (MQ=255)
                        gTCATCAGCTCTAATAACGTGAAGTCAATAGCCTgg      >  1:5804324/1‑36 (MQ=255)
                        gTCATCAGCTCTAATAACGTGAAGTCAATAGCCTgg      >  1:2110137/1‑36 (MQ=255)
                          catcaGCTCTAATAACGTGAAGTCAATAGCCTGGCg    >  1:3475016/1‑36 (MQ=255)
                            tcaGCTCTAATAACGTGAAGTCAATAGCCTGGCGTa  <  1:1171271/36‑1 (MQ=255)
                                |                               
CTGGCTTCTGCTTTCGCTGTTAACGTCATCAGATCTAATAACGTGAAGTCAATAGCCTGGCGTA  >  REL606/1248348‑1248411

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG