Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,248,380 A→C 87.0% 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→C87.0% 9.0 ‑0.4 ‑0.353885dhaMfused predicted dihydroxyacetone‑specific PTS enzymes: HPr component/EI component
Reads supporting (aligned to +/- strand):  new base (26/20):  ref base (5/2):  total (31/22)
Fisher's exact test strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.08e-01

CGCTGGCTTCTGCTTTCGCTGTTAACGTCATCAGATCTAATAACGTGAAGTCAATAGCCTGGCGTAATC  >  REL606/1248346‑1248414
                                  |                                  
cGCTGGCTTCTGCTTTCGCTGTTAACGTCATCAGAt                                   >  1:5909163/1‑36 (MQ=255)
 gCTGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTc                                  >  1:8101979/1‑36 (MQ=255)
  cTGGCTTCTGCTTTCGCTGTTAACGTCATCAGCTCt                                 <  1:3896446/36‑1 (MQ=255)
    ggCTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaa                               <  1:2672892/36‑1 (MQ=255)
    ggCTTCTGCTTTCGCTGTTAACGTCATCAGATCtaa                               >  1:8010840/1‑36 (MQ=255)
      cTTCTGCTTTCGCTGTTAACGTCATCAGCTCtaata                             >  1:7865718/1‑36 (MQ=255)
       ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                            >  1:6338502/1‑36 (MQ=255)
       ttCTGCTTTCGCTGTTAACGTCATCAGCTCtaataa                            <  1:2153469/36‑1 (MQ=255)
        tCTGCTTTCGCTGTTAACGTCATCAGCTCTAATAAc                           <  1:5974633/36‑1 (MQ=255)
        tCTGCTTTCGCTGTTAACGTCATCAGCTCTAATAAc                           <  1:5773722/36‑1 (MQ=255)
        tCTGCTTTCGCTGTTAACGTCATCAGATCTAATAAc                           <  1:1004442/36‑1 (MQ=255)
         cTGCTTTCGCTGTTAACGTCATCAGCTCTAATAACg                          <  1:2244480/36‑1 (MQ=255)
         cTGCTTTCGCTGTTAACGTCATCAGCTCTAATAACg                          >  1:7205757/1‑36 (MQ=255)
          tGCTTTCGCTGTTAACGTCATCAGCTCTAATAACGt                         <  1:6395997/36‑1 (MQ=255)
          tGCTTTCGCTGTTAACGTCATCAGCTCTAATAACGt                         <  1:1329848/36‑1 (MQ=255)
           gCTTTCGCTGTTAACGTCCTCAGCTCTAATAACGTg                        >  1:4858712/1‑36 (MQ=37)
           gCTTTCGCTGTTAACGTCATCAGCTCTAATAACGTg                        <  1:1469519/36‑1 (MQ=255)
           gCTTTCGCTGTTAACGTCATCAGCTCTAATAACGTg                        >  1:3645526/1‑36 (MQ=255)
            cTTTCGCTGTTAACGTCATCAGCTCTAATAACGTGa                       <  1:2660428/36‑1 (MQ=255)
             tttCGCTGTTAACGTCATCAGCTCTAATAACGTGaa                      >  1:6873888/1‑36 (MQ=255)
             tttCGCTGTTAACGTCATCAGCTCTAATAACGTGaa                      <  1:5551966/36‑1 (MQ=255)
              ttCGCTGTTAACGTCATCAGCTCTAATAACGTGAAg                     >  1:46095/1‑36 (MQ=255)
               tCGCTGTTAACGTCATCAGCTCTAATAACGTGAAGt                    >  1:1198397/1‑36 (MQ=255)
                cGCTGTTAACGTCATCAGCTCTAATAACGTGAAGTc                   >  1:3005469/1‑36 (MQ=255)
                cGCTGTTAACGTCATCAGCCCTAATAACGTGAAGTc                   >  1:5674847/1‑36 (MQ=37)
                 gCTGTTAACGTCATCAGCTCTAATAACGTTAAGTCa                  >  1:2549971/1‑36 (MQ=37)
                 gCTGTTAACGTCATCAGCTCTAATAACGTGAAGTCa                  >  1:5602141/1‑36 (MQ=255)
                 gCTGTTAACGTCATCAGCTCTAATAACGTGAAGTCa                  >  1:7264511/1‑36 (MQ=255)
                 gCTGTTAACGTCATCAGCTCTAATAACGTGAAGTCa                  >  1:447254/1‑36 (MQ=255)
                 gCTGTTAACGTCATCAGCTCTAATAACGTGAAGTCa                  <  1:8030299/36‑1 (MQ=255)
                   tGTTAACGTCATCAGCTCTAATAACGTGAAGTCAAt                <  1:659771/36‑1 (MQ=255)
                   tGTTAACGTCATCAGCTCTAATAACGTGAAGTCAAt                >  1:6867890/1‑36 (MQ=255)
                    gTTAACGTCATCAGCTCTAATAACGTGAAGTCAATa               <  1:4864430/36‑1 (MQ=255)
                    gTTAACGTCATCAGCTCTAATAACGTGAAGTCAATa               <  1:3405607/36‑1 (MQ=255)
                    gTTAACGTCATCAGATCTAATAACGTGAAGTCAATa               <  1:3143545/36‑1 (MQ=255)
                     ttAACGTCATCAGATCTAATTATGTGAAGTCAATAg              >  1:722057/1‑36 (MQ=38)
                      tAACGTCATCAGCTCTAATAACGTGAAGTCAATAGc             >  1:4170041/1‑36 (MQ=255)
                      tAACGTCATCAGCTCTAATAACGTGAAGTCAATAGc             >  1:6268168/1‑36 (MQ=255)
                       aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc            >  1:930115/1‑36 (MQ=255)
                       aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc            <  1:861375/36‑1 (MQ=255)
                       aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc            >  1:4327329/1‑36 (MQ=255)
                       aaCGTCATCAGCTCTAATAACGTGAAGTCAATAGcc            >  1:3550086/1‑36 (MQ=255)
                         cGTCATCAGCTCTAATAACGTGAAGTCAATAGCCTg          >  1:4211057/1‑36 (MQ=255)
                         cGTCATCAGCTCTAATAACGTGAAGTCAATAGCCTg          >  1:3190524/1‑36 (MQ=255)
                          gTCATCAGCTCTAATAACGTGAAGTCAATAGCCTgg         <  1:6610612/36‑1 (MQ=255)
                           tcatcaGCTCTAATAACGTGAAGTCAATAGCCTGGc        <  1:5131620/36‑1 (MQ=255)
                            catcaGCTCTAATAACGTGAAGTCAATAGCCTGGCg       >  1:725315/1‑36 (MQ=255)
                            catcaGCTCTAATAACGTGAAGTCAATAGCCTGGCg       >  1:2582161/1‑36 (MQ=255)
                            catcaGCTCTAATAACGTGAAGTCAATAGCCTGGCg       <  1:2248354/36‑1 (MQ=255)
                            catcaGATCTAATAACGTGAAGTCAATAGCCTGGCg       >  1:5165334/1‑36 (MQ=255)
                               caGCTCTAATAACGTGAAGTCAATAGCCTGGCGTaa    <  1:7112003/36‑1 (MQ=255)
                               caGCTCTAATAACGTGAAGTCAATAGCCTGGCGTaa    >  1:5924067/1‑36 (MQ=255)
                                 gATCTAATAACGTGAAGTCAATAGCCTGGCGTAATc  >  1:7588945/1‑36 (MQ=255)
                                  |                                  
CGCTGGCTTCTGCTTTCGCTGTTAACGTCATCAGATCTAATAACGTGAAGTCAATAGCCTGGCGTAATC  >  REL606/1248346‑1248414

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG