Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 800,303 T→C 65.5% A342A (GCA→GCG 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,3030TC65.5% 35.2 / 27.5 52A342A (GCA→GCGdhaHfused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component
Reads supporting (aligned to +/- strand):  ref base T (10/8);  new base C (20/14);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTG  >  minE/800242‑800370
                                                             |                                                                   
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:3264195/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:125158/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:2688270/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:1995181/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:799847/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:2109934/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:1729271/1‑69 (MQ=255)
catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                              >  1:1748912/1‑69 (MQ=255)
      ccAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGtt                                                        <  1:3308499/69‑1 (MQ=255)
      ccAGTTGCTGGTATTGATGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGtt                                                        <  1:2600959/69‑1 (MQ=255)
               ggTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGa                                                >  1:3068163/1‑68 (MQ=255)
               ggTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGa                                                >  1:249151/1‑68 (MQ=255)
               ggTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGa                                                >  1:2422024/1‑68 (MQ=255)
               ggTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATACCGCCGTGCAATGTTCATGCTGa                                                >  1:2211207/1‑68 (MQ=255)
                   ttGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAg                                           <  1:3053374/69‑1 (MQ=255)
                   ttGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAg                                           <  1:2050571/69‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAGGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:2255733/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:2207841/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:3168753/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:3098198/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:1835241/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:3257582/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:1901913/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:502002/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:1693485/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:69850/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:753576/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:1418188/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:917321/68‑1 (MQ=255)
                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAGGAGt                                          <  1:1193100/68‑1 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGc                                                   >  1:1384118/1‑53 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:946438/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:566980/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:588039/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:1017667/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:2577683/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:1408886/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:1412945/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:1464373/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:151089/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:1697669/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:1761904/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:2017035/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:2023911/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:2212591/1‑59 (MQ=255)
                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCCGCCGTGCAATGTTCATGCTGAAgg                                             >  1:2793489/1‑59 (MQ=255)
                             aGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgcc                                 >  1:3227443/1‑69 (MQ=255)
                             aGTTCTTTAAGAACTTGCTGCAAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgcc                                 >  1:451122/1‑69 (MQ=255)
                                                   aGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatca            <  1:3277839/68‑1 (MQ=255)
                                                   aGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatca            <  1:2539358/68‑1 (MQ=255)
                                                            cTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtg  <  1:1721774/69‑1 (MQ=255)
                                                            cTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtg  <  1:1489694/69‑1 (MQ=255)
                                                             tGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtg  <  1:1778380/68‑1 (MQ=255)
                                                             |                                                                   
CATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTG  >  minE/800242‑800370

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