Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 913,152 G→A 41.6% A250A (GCC→GCT ycjW ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE913,1520GA41.6% 26.9 / 50.3 72A250A (GCC→GCTycjWpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (36/6);  new base A (23/7);  total (59/13)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTC  >  minE/913092‑913217
                                                            |                                                                 
aaTAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAacac                                                           <  1:1824873/69‑1 (MQ=255)
aaTAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAacac                                                           <  1:1887910/69‑1 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCTCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:3223044/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2283437/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:3349445/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2193029/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:3326042/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:363657/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2395207/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2564060/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2566741/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2735041/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2872244/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:2819403/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:1846566/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:455441/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:995923/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:1794699/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:458528/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:667197/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:1160422/1‑68 (MQ=255)
    gCAATGTTTTTCTCTGCCGCCGCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGc                                                        >  1:904660/1‑68 (MQ=255)
    gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                                        >  1:1747874/1‑68 (MQ=255)
    gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                                        >  1:2170423/1‑68 (MQ=255)
    gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                                        >  1:1444230/1‑68 (MQ=255)
    gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                                        >  1:1168473/1‑68 (MQ=255)
    gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                                        >  1:934858/1‑68 (MQ=255)
    gCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGc                                                        >  1:3228163/1‑68 (MQ=255)
          tttttCTCTGCTGCCGCAGGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:2059588/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGcc                                                       >  1:409359/1‑63 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:332196/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:952814/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:993620/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:2890064/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:811409/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 <  1:800806/69‑1 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:2853894/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 <  1:26251/69‑1 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1809245/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1142495/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1172037/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1250054/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:129831/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1489955/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1887843/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:2210100/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:2387927/1‑69 (MQ=255)
          tttttCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGc                                                 >  1:1044638/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:2549699/69‑1 (MQ=255)
                         gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt                                  <  1:2744521/69‑1 (MQ=255)
                         gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt                                  <  1:1178557/69‑1 (MQ=255)
                         gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt                                  <  1:759275/69‑1 (MQ=255)
                         gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt                                  <  1:703026/69‑1 (MQ=255)
                         gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt                                  <  1:1772127/69‑1 (MQ=255)
                         gCACGATAAATACTCATTAATTTCAGGCTGTCAGTAGCAAACACTGCCTCAGGCAACGGTGACTGGCtt                                  <  1:1961725/69‑1 (MQ=255)
                                 aaTGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTg                          >  1:1545288/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:66212/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:2783117/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:870538/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:1290186/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:2982/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:940310/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:400632/1‑69 (MQ=255)
                                                     tGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAg      >  1:2167684/1‑69 (MQ=255)
                                                         aGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTc  >  1:3145461/1‑69 (MQ=255)
                                                         aGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTc  >  1:3098196/1‑69 (MQ=255)
                                                            |                                                                 
AATAGCAATATTTTTCTCTGCTGCCGCACGATAAATGCTCATTAATTTCAGGCTGTCAGTGGCAAACACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTC  >  minE/913092‑913217

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: