Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,634,919 G→C 53.1% G19A (GGG→GCG) ‡ ugpB → glycerol‑3‑phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,634,9190GC53.1% 0.2 / 35.0 32G19A (GGG→GCG) ‡ugpBglycerol‑3‑phosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (4/11);  new base C (11/6);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

AACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTA  >  minE/2634860‑2634978
                                                           |                                                           
aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                                    >  1:3316586/1‑69 (MQ=255)
aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                                    >  1:3386718/1‑69 (MQ=255)
aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                                    >  1:945423/1‑69 (MQ=255)
aaCGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGcaca                                                    >  1:1390509/1‑69 (MQ=255)
  caaTGAAACCGTTACGTTATACAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGcgca                                                    >  1:36396/3‑67 (MQ=255)
  caaTGAAACCGTTACGTTATACAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGcgca                                                    >  1:62210/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:1661524/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)
                   tataCAGCTTCAGCACTGGCGCTCGGACTGGCTTTAATGGCAAATGCGCAAGCAGTGACGACCATTCCg                                 >  1:1303927/1‑69 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:636555/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:593123/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:3365996/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:2772621/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:135184/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:1508770/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:1531690/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTc                      <  1:1047574/68‑1 (MQ=255)
                               gCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTCACGACCATTCCGTTCTGGCATTc                      <  1:2989662/68‑1 (MQ=255)
                                             aCTGGCTTTAATGGCAATTGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc        <  1:1545621/68‑1 (MQ=255)
                                             aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc        <  1:965446/68‑1 (MQ=255)
                                             aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc        <  1:2962540/68‑1 (MQ=255)
                                             aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc        <  1:20607/68‑1 (MQ=255)
                                             aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc        <  1:1954225/68‑1 (MQ=255)
                                             aCTGGCTTTAATGGCAAATGCGCAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAAc        <  1:514084/68‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:1666810/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:884900/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:1575609/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:1650626/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:3263388/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:2650858/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:3017730/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:2699041/61‑1 (MQ=255)
                                                          gggAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:2785442/61‑1 (MQ=255)
                                                          gggAATGCAAAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTa  <  1:734330/61‑1 (MQ=255)
                                                           |                                                           
AACGATGAAACCGTTACATTATACAGCTTCAGCACTGGCGCTCGGACTGGCGTTAATGGGGAATGCACAGGCAGTGACGACCATTCCGTTCTGGCATTCTATGGAAGGGGAACTGGGTA  >  minE/2634860‑2634978

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: