Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,369,021 G→A 76.8% L460L (TTG→TTA 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,0210GA76.8% 72.4 / 17.7 60L460L (TTG→TTAfadBfused 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 G (10/4);  new base A (20/26);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCTGGAACGCCCGGAAAACTTCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGAGCAA  >  minE/2368958‑2369089
                                                               |                                                                    
gcTGGAACGCCCGGAAAACTTCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGa                                                                 <  1:2004131/69‑1 (MQ=255)
 cTGGAACGCCCGGAAAACTTCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGa                                                                 <  1:2078232/68‑1 (MQ=255)
                aaCTTCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCgag                                                 <  1:2659177/69‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:876410/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1019684/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:2521034/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:913723/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:2359527/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:2931269/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:2976820/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1830213/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:446945/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1509488/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1362401/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:52559/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:825130/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:135934/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1268161/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1113003/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:630827/54‑1 (MQ=255)
                    tCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:1115754/54‑1 (MQ=255)
                     cTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGAAatta                                                            <  1:2812948/53‑1 (MQ=255)
                               cACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGaa                                  >  1:932399/1‑69 (MQ=255)
                               cACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGaa                                  >  1:195740/1‑69 (MQ=255)
                               cACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGaa                                  >  1:1222829/1‑69 (MQ=255)
                               cACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGaa                                  >  1:738849/1‑69 (MQ=255)
                                aCTTCTTTAACCCGGTCCACCGAATGCCGTTAGTAGATatta                                                            <  1:1715002/42‑1 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:254709/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:2834154/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:284983/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:3087430/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:596048/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:3125453/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:3311788/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:3392686/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:926369/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:1069276/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:1090254/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:959251/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:1501210/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:1665160/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:919733/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:2323715/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:2622132/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:2382058/1‑51 (MQ=255)
                                    cttTAACCCGGTCCACCGAATGCCGTTAGTAGAAATAATTCGCGGCGAGaa                                               >  1:1305552/1‑51 (MQ=255)
                                    cttTAACCCGGTCAACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGaa                                               >  1:1587532/1‑51 (MQ=255)
                                           ccGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGaaa                      <  1:3345192/69‑1 (MQ=255)
                                               tCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGtcg                  <  1:2951936/69‑3 (MQ=255)
                                               tCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGtcg                  <  1:2781882/69‑3 (MQ=255)
                                               tCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGtcg                  <  1:906488/69‑3 (MQ=255)
                                               tCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGtcg                  <  1:190262/69‑3 (MQ=255)
                                               tCCACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGtcg                  <  1:306670/69‑3 (MQ=255)
                                                ccACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATcgc                          >  1:1114317/1‑60 (MQ=255)
                                                ccACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGt                 >  1:3365763/1‑69 (MQ=255)
                                                ccACCGAATGCCATTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGt                 >  1:3087126/1‑69 (MQ=255)
                                                 cACCGAATGCCGTTAGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGtcg                  <  1:2777907/67‑3 (MQ=255)
                                                        tGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTgg          >  1:2469461/1‑68 (MQ=255)
                                                        tGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTgg          >  1:2234608/1‑68 (MQ=255)
                                                        tGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTgg          >  1:2952478/1‑68 (MQ=255)
                                                               ggTTGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGa      >  1:1805167/1‑65 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGa      >  1:2688705/1‑65 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGa      >  1:2852820/1‑65 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGa      >  1:2740188/1‑65 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGa      >  1:2130501/1‑65 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGa      >  1:1364220/1‑65 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGAGCaa  <  1:3291249/69‑1 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGAGCaa  <  1:1153365/69‑1 (MQ=255)
                                                               ggTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGAGCaa  <  1:987194/69‑1 (MQ=255)
                                                               |                                                                    
GCTGGAACGCCCGGAAAACTTCTGCGGGATGCACTTCTTTAACCCGGTCCACCGAATGCCGTTGGTAGAAATTATTCGCGGCGAGAAAAGCTCCGACGAAACCATCGCGAAAGTTGTCGCCTGGGCGAGCAA  >  minE/2368958‑2369089

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: