Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 7,496 C→G 32.7% G155A (GGC→GCC)  yaaJ ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE7,4960CG32.7% 54.0 / 35.9 52G155A (GGC→GCC) yaaJpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (26/9);  new base G (10/7);  total (36/16)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7428‑7550
                                                                    |                                                      
ggggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGc                                                        <  1:2000316/69‑1 (MQ=255)
                 aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:2144471/69‑1 (MQ=255)
                 aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:2231443/69‑1 (MQ=255)
                 aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:3074767/69‑1 (MQ=255)
                 aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:679474/69‑1 (MQ=255)
                  aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:901386/68‑1 (MQ=255)
                  aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:2216594/68‑1 (MQ=255)
                  aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:1511347/68‑1 (MQ=255)
                    cTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAg                                    <  1:2900651/69‑1 (MQ=255)
                      cAGGGCGCTGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:740452/69‑1 (MQ=255)
                      cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:44685/69‑1 (MQ=255)
                      cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:164911/69‑1 (MQ=255)
                      cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:2372362/69‑1 (MQ=255)
                               ggCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAg                         <  1:1409568/69‑1 (MQ=255)
                               ggCAACGGCGTTCGCTTGAACTCCCCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAg                         <  1:2131063/69‑1 (MQ=255)
                                    cGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGc                    <  1:3112990/69‑1 (MQ=255)
                                     ggCGTTCTCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:3166875/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:305304/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:845844/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1041534/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1018455/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1287037/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1520955/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:172562/1‑69 (MQ=255)
                                     ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGc                    >  1:2252914/1‑68 (MQ=255)
                                           cGCTTGCACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:439998/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:360526/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2010953/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:293591/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2994938/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:868008/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2027817/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2586750/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:1252826/1‑68 (MQ=255)
                                           cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:817177/1‑68 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:2032330/1‑58 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:638351/1‑58 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:754247/1‑58 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:3028450/1‑58 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCa        >  1:20065/1‑69 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCa        >  1:2753905/1‑69 (MQ=255)
                                                gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCa        >  1:2722561/1‑69 (MQ=255)
                                                 aaCTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGc                    <  1:2667871/56‑1 (MQ=255)
                                                      cGCTGAAAATTATGGCTTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:87435/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATTAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:3320544/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCa                 >  1:1653504/1‑54 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:622299/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:1534329/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:239230/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:3117022/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:302538/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:2610482/1‑69 (MQ=255)
                                                      cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGAGCATCCCCa  >  1:1580875/1‑69 (MQ=255)
                                                                    |                                                      
GGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7428‑7550

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: