Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,542,969 A→G 42.3% G481G (GGT→GGC yfcX ← fused enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, dodecenoyl‑CoA D‑isomerase, and 3‑hydroxyacyl‑CoA dehydrogenase, subunit of anaerobic fatty acid oxidation complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,542,9690AG42.3% 15.0 / 42.9 59G481G (GGT→GGCyfcXfused enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, dodecenoyl‑CoA D‑isomerase, and 3‑hydroxyacyl‑CoA dehydrogenase, subunit of anaerobic fatty acid oxidation complex
Reads supporting (aligned to +/- strand):  ref base A (26/8);  new base G (19/6);  total (45/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AATGTAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGAATA  >  minE/1542903‑1543035
                                                                  |                                                                  
aaTGTAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTAcc                                                                  >  1:2603046/1‑69 (MQ=255)
aaTGTAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTAcc                                                                  >  1:1876937/1‑69 (MQ=255)
    tAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGt                                                              >  1:2542564/1‑69 (MQ=255)
    tAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGt                                                              >  1:935674/1‑69 (MQ=255)
    tAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGt                                                              >  1:669733/1‑69 (MQ=255)
    tAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGt                                                              >  1:1472720/1‑69 (MQ=255)
      aGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGttt                                                            >  1:985332/1‑69 (MQ=255)
      aGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGttt                                                            >  1:64543/1‑69 (MQ=255)
             aaaaTGCGATTGACGTAAAACCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGcc                                                     <  1:1594046/66‑1 (MQ=255)
             aaaaTGCGATTGACGTAAAACCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGcc                                                     <  1:2235984/66‑1 (MQ=255)
             aaaaTGCGATTGACGTAAAACCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGcc                                                     <  1:1807219/66‑1 (MQ=255)
             aaaaTGCGATTGACGTAAAACCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGcc                                                     <  1:419576/66‑1 (MQ=255)
             aaaaTGCGATTGACGTAAAACCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGcc                                                     <  1:850202/66‑1 (MQ=255)
             aaaaTGCGATTGACGTAAAACCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGcc                                                     <  1:127989/66‑1 (MQ=255)
                        gACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGtttt                                               >  1:294478/1‑64 (MQ=255)
                        gACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACtgt                                          >  1:3081777/1‑69 (MQ=255)
                        gACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACtgt                                          >  1:100630/1‑69 (MQ=255)
                        gACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACtgt                                          >  1:2436774/1‑69 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCTGTTTTACTGTGGTa                                      >  1:2302/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGt                                       >  1:1591839/1‑67 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:3179390/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:2683530/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:3124138/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:3401350/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:2765773/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:2667818/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:987926/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:1242514/1‑68 (MQ=255)
                             aaaacCCGGCTTTGTCACGCACGACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTa                                      >  1:1535816/1‑68 (MQ=255)
                              aaaaCCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGc                                    >  1:2199981/1‑69 (MQ=255)
                              aaaaCCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGc                                    >  1:762101/1‑69 (MQ=255)
                              aaaaCCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGc                                    >  1:1851330/1‑69 (MQ=255)
                              aaaaCCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGc                                    >  1:1709498/1‑69 (MQ=255)
                              aaaaCCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGc                                    >  1:153256/1‑69 (MQ=255)
                                       tttGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTg                           <  1:371909/69‑1 (MQ=255)
                                            cacgcacgACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGcc                       <  1:7240/68‑1 (MQ=255)
                                            cacgcacgACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGcc                       <  1:1473321/68‑1 (MQ=255)
                                            cacgcacgACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGcc                       <  1:519301/68‑1 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTgcg                         >  1:2053032/1‑63 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATg                   >  1:2325530/1‑69 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATg                   >  1:957165/1‑69 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATg                   >  1:454040/1‑69 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATg                   >  1:3258680/1‑69 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATg                   >  1:2857139/1‑69 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATg                   >  1:2989547/1‑69 (MQ=255)
                                               gcacgACAATGGGCGTTTTGCCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGAGCCGATg                   >  1:1388735/1‑69 (MQ=255)
                                                   gACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTccc               <  1:1564015/69‑1 (MQ=255)
                                                   gACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTccc               <  1:3062167/69‑1 (MQ=255)
                                                   gACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTccc               <  1:303512/69‑1 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAgg      <  1:2937067/65‑1 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:3349762/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:3161421/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:76299/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:828603/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:1369425/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:1318986/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaata  >  1:1067842/1‑69 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaat   >  1:2277139/1‑68 (MQ=255)
                                                                ttACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGaat   >  1:2568356/1‑68 (MQ=255)
                                                                  |                                                                  
AATGTAAGGCGCTAAGATGCGATTGACGTAAAAACCGGCTTTGTCACGCACGACAATTGGCGTTTTACCCTGTTTTTTCGCCAGTTTTACTGTGGTAGCGATGGTTTGCGCCGATGTCCCCGCATGAGGAATA  >  minE/1542903‑1543035

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