Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,369,370 G→A 52.8% D577N (GAT→AAT)  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,3700GA52.8% ‑2.6 / 28.7 36D577N (GAT→AAT) fadBfused 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 (7/10);  new base A (15/4);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 3.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAAGAAGACAGCAAA  >  minE/2369309‑2369432
                                                             |                                                              
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:2665202/1‑69 (MQ=255)
                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt                           >  1:1218301/1‑69 (MQ=255)
                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt                           >  1:424478/1‑69 (MQ=255)
                              aGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt                           >  1:2937383/1‑69 (MQ=255)
                              aGATTAACGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGttt                           >  1:442899/1‑69 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:3334407/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:26785/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:2847426/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:982027/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:416959/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:905400/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:957317/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:2478603/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:2213119/1‑60 (MQ=255)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGCCt                                 >  1:1100751/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)
                                 ttATCGCGATGCCATCGACGCGCTGTTTAATGCCAACCGCTTTGGTCAGAAGAACGGACt                                 >  1:1025709/1‑60 (MQ=255)
                                            ccATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAaag             <  1:2492583/69‑1 (MQ=255)
                                            ccATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAaag             <  1:2736363/69‑1 (MQ=255)
                                            ccATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAaag             <  1:1832883/69‑1 (MQ=255)
                                            ccATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAaag             <  1:1283554/69‑1 (MQ=255)
                                            ccATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAaag             <  1:1157161/69‑1 (MQ=255)
                                                       cTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAAGAAGACAGCaaa  <  1:1020320/69‑1 (MQ=255)
                                                       cTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAAGAAGACAGCaaa  <  1:2493037/69‑1 (MQ=255)
                                                       cTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAAGAAGACAGCaaa  <  1:2030980/69‑1 (MQ=255)
                                                       cTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAAGAAGACAGCaaa  <  1:2729988/69‑1 (MQ=255)
                                                       cTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAACAAGACAGCaaa  <  1:2579586/69‑1 (MQ=255)
                                                             |                                                              
GGCAGCAGGCTTCCCGCAGCGGATGCAGAAAGATTACCGCGATGCCATCGACGCGCTGTTTGATGCCAACCGCTTTGGTCAGAAGAACGGCCTCGGTTTCTGGCGTTATAAAGAAGACAGCAAA  >  minE/2369309‑2369432

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: