Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,634,911 G→T 31.4% A16A (GCG→GCT ugpB → glycerol‑3‑phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,634,9110GT31.4% 64.7 / 30.4 54A16A (GCG→GCTugpBglycerol‑3‑phosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base T (11/6);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACAAGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAC  >  minE/2634845‑2634972
                                                                  |                                                             
acaacaAGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGtt                                                             <  1:349104/69‑1 (MQ=255)
acaacaAGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGtt                                                             <  1:2964495/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:3038010/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:2042462/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:2054401/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:2057340/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:2130741/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:218670/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:2560513/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:2890521/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:1109508/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:3141206/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:478535/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:531355/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:533871/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:691408/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:982335/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:99874/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            >  1:2019872/1‑69 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:135429/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:2019942/69‑1 (MQ=255)
 caacaaGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTa                                                            <  1:162796/69‑1 (MQ=255)
               aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                              >  1:3316586/1‑69 (MQ=255)
               aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                              >  1:945423/1‑69 (MQ=255)
               aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                              >  1:3386718/1‑69 (MQ=255)
               aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                              >  1:1390509/1‑69 (MQ=255)
                 caaTGAAACCGTTACGTTATACAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGcgca                                              >  1:62210/3‑67 (MQ=255)
                 caaTGAAACCGTTACGTTATACAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGcgca                                              >  1:36396/3‑67 (MQ=255)
                           gTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTgacgac                                  <  1:1941179/69‑1 (MQ=255)
                             tACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCa                                <  1:1149695/69‑1 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTcc                            >  1:3008322/1‑68 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTcc                            >  1:1229075/1‑68 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCg                           >  1:1619304/1‑69 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCg                           >  1:706968/1‑69 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCg                           >  1:1227032/1‑69 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCg                           >  1:2176066/1‑69 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCg                           >  1:3153349/1‑69 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCg                           >  1:1661524/1‑69 (MQ=255)
                                  tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAAGCAGTGACGACCATTCCg                           >  1:1303927/1‑69 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:1047574/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:636555/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:593123/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:135184/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:3365996/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:1508770/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:1531690/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                <  1:2772621/68‑1 (MQ=255)
                                              gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTCACGACCATTCCGTTCTGGCATTc                <  1:2989662/68‑1 (MQ=255)
                                                            aCTGGCTTTAATGGCAATTGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc  <  1:1545621/68‑1 (MQ=255)
                                                            aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc  <  1:514084/68‑1 (MQ=255)
                                                            aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc  <  1:20607/68‑1 (MQ=255)
                                                            aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc  <  1:2962540/68‑1 (MQ=255)
                                                            aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc  <  1:1954225/68‑1 (MQ=255)
                                                            aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc  <  1:965446/68‑1 (MQ=255)
                                                                  |                                                             
ACAACAAGAGAGATAAACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAC  >  minE/2634845‑2634972

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: