Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 7,498 A→G 31.5% Y154Y (TAT→TAC yaaJ ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE7,4980AG31.5% 40.2 / 23.1 52Y154Y (TAT→TACyaaJpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (26/9);  new base G (9/7);  total (36/16)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7445‑7550
                                                     |                                                    
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)
                                                     |                                                    
AAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7445‑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: