Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,248,380 A→C 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% 96.4 34885dhaMfused predicted dihydroxyacetone‑specific PTS enzymes: HPr component/EI component
Reads supporting (aligned to +/- strand):  new base (16/16):  ref base (0/0):  total (18/16)

CTGGCTTCTGCTTTCGCTGTTAACGTCATCAGATCTAATAACGTGAAGTCAATAGCCTGGCGT  >  REL606/1248348‑1248410
                                |                              
cTGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTCt                             >  1:490537/1‑36 (MQ=255)
cTGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTCt                             <  2:1826723/36‑1 (MQ=255)
 tGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTCta                            >  2:1567754/1‑36 (MQ=255)
  ggCTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaa                           >  1:2731878/1‑36 (MQ=255)
   gCTTCTGCTTTCGCTGTTAACGTCATCAGTTCtaat                          >  1:3730187/1‑36 (MQ=255)
    cTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaata                         <  2:2486006/36‑1 (MQ=255)
    cTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaata                         >  2:3424376/1‑36 (MQ=255)
     ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                        >  1:1717733/1‑36 (MQ=255)
     ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                        <  1:1930320/36‑1 (MQ=255)
     ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                        >  1:3427656/1‑36 (MQ=255)
      tCTGCTTTCGCTGTTAACGTCATCAGCTCTAATAAc                       <  2:1364637/36‑1 (MQ=255)
       cTGCTTTCGCTGTTAACGTCATCAGTTCtaataatg                      >  1:590347/1‑34 (MQ=37)
       cTGCTTTCGCTGTTAACGTCATCAGCTCTAACAACg                      >  1:2767545/1‑36 (MQ=37)
       cTGCTTTAGATGTTAACGTCATCAGCTCTAATAACg                      <  1:237367/36‑1 (MQ=25)
        ttcTTTCGCTGTTAACGTCATCAGCTCTAATAACGt                     <  1:390440/34‑1 (MQ=37)
         gCTTTCGCTGTTAACGTCATCAGCTCTAATAACGTg                    <  1:1874264/36‑1 (MQ=255)
          cTTTCGCTGTTAACGTCATCAGCTCTAATAACGTGa                   >  1:2426042/1‑36 (MQ=255)
            ttCGCTGTTAACGTCATCAGCTCTAATAACGTGAAg                 >  1:225033/1‑36 (MQ=255)
            ttCGCTGTTAACGTCATCAGCTCTAATAACGTGAAg                 <  1:749685/36‑1 (MQ=255)
               gCTGTTAACGTCATCAGCTCTAATAACGCGAAGTCa              >  2:887531/1‑36 (MQ=37)
                 tGTTAACGTCATCAGCTGTAATAACGTGAAGTCAAt            <  2:2895710/36‑1 (MQ=37)
                  gTTAACGTCATCAGCTCTAATAACGTGAAGTCAATa           <  2:974576/36‑1 (MQ=255)
                   ttAACGTCATCAGCTCTAATAACGTGAAGTCAATAg          <  1:1991602/36‑1 (MQ=255)
                     aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc        >  2:2041645/1‑36 (MQ=255)
                     aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc        <  1:3325816/36‑1 (MQ=255)
                      aCGTCATCAGCTCTAATAACGTGAAGTCAATAGCCt       <  1:386732/36‑1 (MQ=255)
                      aCGTCATCAGCTCTAATAACGTGAAGTCAATAGCCt       >  1:1073422/1‑36 (MQ=255)
                       cGTCATCAGCTCTAATAACGTGAAGTCAATAGACTg      >  1:3632350/1‑36 (MQ=255)
                        gTCATCAGCTCTAATAACGTCAAGTCAACACCCTgg     >  2:2424303/1‑36 (MQ=19)
                         tcatcaGCTCTAATAACGTGAAGTCAATAGCCTGGc    >  2:1258534/1‑36 (MQ=255)
                         tcatcaGCTCTAATAACGTGAAGTCAATAGCCTGGc    <  2:2421195/36‑1 (MQ=255)
                         tcatcaGCTCTAATAACGGTAAGTCAATAGCCTGGc    <  2:3458009/36‑1 (MQ=25)
                           atcaGCTCTAATAACGTGAAGTCAATAGCCTGGCGt  >  1:757064/1‑36 (MQ=255)
                           atcaGCTCTAATAACGTGAAGTCAATAGCCTGGCGt  <  1:2648148/36‑1 (MQ=255)
                                |                              
CTGGCTTCTGCTTTCGCTGTTAACGTCATCAGATCTAATAACGTGAAGTCAATAGCCTGGCGT  >  REL606/1248348‑1248410

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG