Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 913,128 G→A 31.8% S258S (AGC→AGT ycjW ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE913,1280GA31.8% 85.2 / 45.6 94S258S (AGC→AGTycjWpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (47/17);  new base A (23/7);  total (70/24)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CACCACCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTG  >  minE/913071‑913193
                                                         |                                                                 
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:1313914/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:1873007/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:3338431/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:1608838/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:526963/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:57940/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:1289037/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:299406/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:2634769/1‑67 (MQ=255)
caccacCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:258796/1‑67 (MQ=255)
caccacAGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAAt                                                          >  1:1571484/1‑67 (MQ=255)
    accGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCAt                                                              >  1:2083578/1‑59 (MQ=255)
    accGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCa                                                      >  1:1513530/1‑67 (MQ=255)
    accGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAgg                                                    >  1:2144948/1‑69 (MQ=255)
    accGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAgg                                                    >  1:3071148/1‑69 (MQ=255)
    accGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAgg                                                    >  1:994302/1‑69 (MQ=255)
    accGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAgg                                                    >  1:3243400/1‑69 (MQ=255)
       gCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCt                                                  >  1:982383/1‑68 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:2299572/69‑1 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:2697844/69‑1 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:105674/69‑1 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:313307/69‑1 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:3173068/69‑1 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:1151123/69‑1 (MQ=255)
        cTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                <  1:610724/69‑1 (MQ=255)
         tAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGt                                                >  1:980715/1‑68 (MQ=255)
         tAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTc                                               >  1:3307896/1‑69 (MQ=255)
         tAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTc                                               <  1:1868210/69‑1 (MQ=255)
           aCTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTc                                               <  1:2766441/67‑1 (MQ=255)
           aCTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTc                                               <  1:2888612/67‑1 (MQ=255)
               ctgcGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTc                                               <  1:1803285/63‑1 (MQ=255)
                     aaTAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAacac                                   <  1:1887910/69‑1 (MQ=255)
                     aaTAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAacac                                   <  1:1824873/69‑1 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCTCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:3223044/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:458528/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:3326042/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:667197/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:455441/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2872244/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2819403/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:363657/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:904660/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2735041/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:3349445/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:995923/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2564060/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:1160422/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:1794699/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:1846566/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:1971371/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:205288/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2090814/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2193029/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2283437/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2566741/1‑68 (MQ=255)
                         gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                >  1:2395207/1‑68 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTc                                               >  1:787630/1‑53 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                >  1:1444230/1‑68 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                >  1:1747874/1‑68 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                >  1:934858/1‑68 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                >  1:3228163/1‑68 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                >  1:1168473/1‑68 (MQ=255)
                         gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                >  1:2170423/1‑68 (MQ=255)
                               tttttCTCTGCTGCCGCAGGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:2059588/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGcc                               >  1:409359/1‑63 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:952814/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:993620/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1172037/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:811409/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         <  1:800806/69‑1 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1142495/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1250054/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1044638/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:129831/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:2210100/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1809245/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:2853894/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:2890064/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:2387927/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1887843/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         <  1:26251/69‑1 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:332196/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:1489955/1‑69 (MQ=255)
                               tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                         >  1:2605905/1‑69 (MQ=255)
                                      ctgctgCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCGGGCAACGGt                   <  1:434800/68‑1 (MQ=255)
                                      ctgctgCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGt                   <  1:521871/68‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:759275/69‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:703026/69‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:1178557/69‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:1772127/69‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:1961725/69‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:2549699/69‑1 (MQ=255)
                                              gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt          <  1:2744521/69‑1 (MQ=255)
                                                      aaTGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTg  >  1:1545288/1‑69 (MQ=255)
                                                         |                                                                 
CACCACCGCTAACTGCTGCGGAATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTG  >  minE/913071‑913193

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: