Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,795,002 A→G 37.2% R113R (CGT→CGC plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,795,0020AG37.2% 21.0 / 24.5 43R113R (CGT→CGCplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (17/10);  new base G (7/9);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG  >  minE/2794938‑2795054
                                                                |                                                    
gCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACgg                                                   <  1:1439501/68‑1 (MQ=255)
         cgcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCga                                         <  1:2881599/69‑1 (MQ=255)
          gcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCga                                         <  1:1329716/68‑1 (MQ=255)
          gcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCga                                         <  1:2185980/68‑1 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:1026679/1‑69 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:723105/1‑69 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:662286/1‑69 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:647822/1‑69 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:606411/1‑69 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:583444/1‑69 (MQ=255)
                aaaCATCACCGATACAGGCACCATCTGCACATCCAGATTTGGGTTGCTGCGGTGCAAATCGAGATAGTc                                  >  1:1023091/1‑69 (MQ=255)
                             acTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTa                     <  1:2883423/69‑1 (MQ=255)
                                ggCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATa                  >  1:1524510/1‑69 (MQ=255)
                                   aCCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGa                <  1:914362/68‑1 (MQ=255)
                                     cATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGActc             <  1:954128/69‑1 (MQ=255)
                                      aTCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGActc             <  1:1190734/68‑1 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGa                             >  1:906074/1‑49 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:2166989/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:73809/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:2880894/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:2093799/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:1763832/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:2926163/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:1861612/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:567851/1‑69 (MQ=255)
                                         tGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttct         >  1:407555/1‑69 (MQ=255)
                                           cacaTCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAg                         <  1:2228424/51‑1 (MQ=255)
                                               tCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCttt       >  1:2436568/1‑65 (MQ=255)
                                               tCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCgg    >  1:514240/1‑68 (MQ=255)
                                               tCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGc   >  1:89222/1‑69 (MQ=255)
                                                 cAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCg  >  1:422816/1‑68 (MQ=255)
                                                 cAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCg  >  1:469493/1‑68 (MQ=255)
                                                 cAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCg  >  1:278334/1‑68 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:564720/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:3403136/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:2034429/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:2906176/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:2062212/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:689842/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:2696633/43‑1 (MQ=255)
                                                        gggTTGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:2308448/43‑1 (MQ=255)
                                                        gggTCGCTGCGGTGCAAATCGAGATAGTCGTGGAACAGCTTaa                    <  1:3164934/43‑1 (MQ=255)
                                                             gCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGa                <  1:2317326/42‑1 (MQ=255)
                                                                |                                                    
GCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCG  >  minE/2794938‑2795054

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