Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,114,060 C→A 45.7% L248L (CTG→CTT pps ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,114,0600CA45.7% 6.5 / 33.5 35L248L (CTG→CTTppsphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (17/2);  new base A (10/6);  total (27/8)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCA  >  minE/1114003‑1114120
                                                         |                                                            
cATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGttt                                                   >  1:1133358/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:2620023/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:453238/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:3377156/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:2989575/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:2751304/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:2718366/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:2447483/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:1556679/1‑69 (MQ=255)
 aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                                  >  1:1667491/1‑69 (MQ=255)
 aTGCGGATTTTTTTAGACCCCATGGTGCGGCGCAAAATAGCCGGGCGATTCGCAGCAAGTGTCGGttt                                                   >  1:2762440/1‑68 (MQ=255)
                        ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt                           >  1:685895/1‑69 (MQ=255)
                        ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt                           >  1:1340944/1‑69 (MQ=255)
                        ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt                           >  1:1404951/1‑69 (MQ=255)
                        ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt                           >  1:3060073/1‑69 (MQ=255)
                        ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCTGGtt                           <  1:2503330/69‑1 (MQ=255)
                        ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt                           <  1:2171284/69‑1 (MQ=255)
                        ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt                           <  1:3202119/69‑1 (MQ=255)
                        ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt                           <  1:2886849/69‑1 (MQ=255)
                        ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt                           <  1:275483/69‑1 (MQ=255)
                        ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt                           <  1:2499302/69‑1 (MQ=255)
                                 cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc                  >  1:1120909/1‑69 (MQ=255)
                                 cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc                  >  1:3033651/1‑69 (MQ=255)
                                 cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc                  >  1:3121984/1‑69 (MQ=255)
                                 cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc                  >  1:2008341/1‑69 (MQ=255)
                                        gCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACg           >  1:729751/1‑69 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCa                    <  1:2296728/50‑1 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCa                    <  1:2228181/50‑1 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCa  >  1:1637574/1‑68 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCa  >  1:245147/1‑68 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCa  >  1:1715707/1‑68 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCa  >  1:2301383/1‑68 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCa  >  1:530686/1‑68 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCa  >  1:901100/1‑68 (MQ=255)
                                                  tCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTAAACCGCACCCTGCACGACCATCTCa  >  1:3107696/1‑68 (MQ=255)
                                                         |                                                            
CATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCA  >  minE/1114003‑1114120

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: