Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,369,345 C→T 41.0% Y568Y (TAC→TAT fadB → fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,369,3450CT41.0% 15.9 / 30.0 47Y568Y (TAC→TATfadBfused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (11/16);  new base T (15/4);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

ATTGATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTT  >  minE/2369277‑2369407
                                                                    |                                                              
aTTGATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTAc                                                                >  1:3107345/1‑69 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             >  1:3169200/1‑69 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             <  1:2947032/69‑1 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             <  1:512203/69‑1 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             <  1:2611284/69‑1 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             <  1:2597881/69‑1 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             >  1:540492/1‑69 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             >  1:628611/1‑69 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             <  1:1308427/69‑1 (MQ=255)
   gATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACcgc                                                             >  1:1380228/1‑69 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGcc                                                       <  1:1017520/69‑1 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGcc                                                       <  1:1720250/69‑1 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGcc                                                       <  1:95026/69‑1 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGcc                                                       <  1:3115188/69‑1 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGcc                                                       <  1:1346062/69‑1 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGcc                                                       <  1:2418064/69‑1 (MQ=255)
         gcgcATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGc                                                        <  1:3059809/68‑1 (MQ=255)
                     cAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTg                                           >  1:316325/1‑69 (MQ=255)
                        gCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGttt                                        <  1:569179/69‑1 (MQ=255)
                        gCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGttt                                        <  1:3128914/69‑1 (MQ=255)
                        gCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGttt                                        <  1:1493067/69‑1 (MQ=255)
                             cATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGttt                                        <  1:1263899/64‑1 (MQ=255)
                                ggCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCaa                                >  1:2549185/1‑69 (MQ=255)
                                 gcagcaGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAAc                               >  1:1570903/1‑69 (MQ=255)
                                                   cGGATGCAGAAAGATTATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAAc             <  1:545761/69‑1 (MQ=255)
                                                   cGGATGCAGAAAGATTATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAAc             <  1:2019229/69‑1 (MQ=255)
                                                   cGGATGCAGAAAGATTATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAAc             <  1:3158792/69‑1 (MQ=255)
                                                   cGGATGCAGAAAGATTATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAAc             <  1:2966025/69‑1 (MQ=255)
                                                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt  >  1:424478/1‑69 (MQ=255)
                                                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt  >  1:2937383/1‑69 (MQ=255)
                                                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt  >  1:2665202/1‑69 (MQ=255)
                                                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt  >  1:1218301/1‑69 (MQ=255)
                                                              aGATTAACGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt  >  1:442899/1‑69 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:905400/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:957317/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:982027/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:2847426/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:416959/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:3334407/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:26785/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:2478603/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:2213119/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:2025818/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:2024998/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:2004314/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:1950839/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt        >  1:1100751/1‑60 (MQ=255)
                                                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGACt        >  1:1025709/1‑60 (MQ=255)
                                                                    |                                                              
ATTGATACCGCGCATCACGCTCAGGCTGTCATGGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTT  >  minE/2369277‑2369407

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: