Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 800,666 C→T 36.9% G221G (GGG→GGA dhaH ← fused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE800,6660CT36.9% 51.1 / 47.9 81G221G (GGG→GGAdhaHfused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base T (14/16);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TCCTCCGTTTCACCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATA  >  minE/800598‑800729
                                                                    |                                                               
tcctccGTTTCACCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGc                                                                 <  1:1488058/69‑1 (MQ=255)
        ttCACCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGc                                                         <  1:2748415/69‑1 (MQ=255)
           aCCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCaa                                                       <  1:1180429/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGcc                                           <  1:2863186/68‑2 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:20176/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2030303/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2049702/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:172245/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2377722/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2443062/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2592635/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2646498/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:1345051/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:1337067/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:1331440/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:2950304/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:3303494/68‑1 (MQ=255)
                       ttCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:690354/68‑1 (MQ=255)
                        tCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:1302530/68‑1 (MQ=255)
                           gccagcTGACGGAAAGCGATCAGTGCTTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTAt                                        <  1:3131887/67‑1 (MQ=255)
                           gccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:657805/65‑1 (MQ=255)
                           gccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:216088/65‑1 (MQ=255)
                           gccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:796020/65‑1 (MQ=255)
                           gccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:992531/65‑1 (MQ=255)
                           gccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:1603421/65‑1 (MQ=255)
                            ccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                          <  1:1872696/64‑1 (MQ=255)
                            ccagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGt                                     <  1:1684866/69‑1 (MQ=255)
                            ccagcTGACGGAAAGCGATCAGTGCATCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGt                                     <  1:2906203/69‑1 (MQ=255)
                                         aGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        <  1:3205199/69‑1 (MQ=255)
                                         aGCGATCAGTGCATCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        <  1:951448/69‑1 (MQ=255)
                                          gCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        <  1:2562522/68‑1 (MQ=255)
                                             aTCAGTGCCTCTTCAGCTTCTGGTCCTTTCGCAATCAGGCGCAGCg                                           <  1:3041843/46‑1 (MQ=255)
                                               cAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCg                              <  1:242628/57‑1 (MQ=255)
                                                aGTGCCTCTTTAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:3004697/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:2732034/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:27070/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:2865226/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:1261778/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:2967217/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:2994995/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:3087548/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:2202601/1‑66 (MQ=255)
                                                aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtagt                    >  1:3418887/1‑66 (MQ=255)
                                                   gCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCa              <  1:2988502/69‑1 (MQ=255)
                                                   gCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCa              <  1:2266185/69‑1 (MQ=255)
                                                   gCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCa              <  1:2882896/69‑1 (MQ=255)
                                                    ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCaa             >  1:3156922/1‑69 (MQ=255)
                                                    ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCaa             >  1:3387166/1‑69 (MQ=255)
                                                    ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCaa             >  1:1494252/1‑69 (MQ=255)
                                                    ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCaa             >  1:2044782/1‑69 (MQ=255)
                                                    ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCaa             >  1:2258429/1‑69 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:2721440/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:2724056/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:2626703/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:2097217/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:1134057/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:1123125/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTgg        >  1:2369143/1‑68 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTg         >  1:2363585/1‑67 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTg         >  1:3423869/1‑67 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTg         >  1:2076108/1‑67 (MQ=255)
                                                          cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTg         >  1:1393026/1‑67 (MQ=255)
                                                            gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGtt      <  1:2131304/68‑1 (MQ=255)
                                                            gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGtt      <  1:3343861/68‑1 (MQ=255)
                                                            gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGt       <  1:115688/67‑1 (MQ=255)
                                                            gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGt       <  1:1778835/67‑1 (MQ=255)
                                                            gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGt       <  1:2643515/67‑1 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:1057596/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:762522/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:1107719/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:1091136/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:973072/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:621116/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:392418/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:3247342/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:1766277/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:2387674/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:2424795/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:2433031/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:1567584/1‑47 (MQ=255)
                                                               tCTGGTCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTg                        >  1:1559646/1‑47 (MQ=255)
                                                               tCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATa  <  1:428910/69‑1 (MQ=255)
                                                                    |                                                               
TCCTCCGTTTCACCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATA  >  minE/800598‑800729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: