Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,661,032 G→A 70.6% L181L (CTG→CTA glpR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,661,0320GA70.6% 117.9 / 56.7 122L181L (CTG→CTAglpRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (34/2);  new base A (70/16);  total (104/18)
Fisher's exact test for biased strand distribution p-value = 9.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

AGTTCCGCCTTGATTTCGGCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAC  >  minE/2660974‑2661083
                                                          |                                                   
aGTTCCGCCTTGATTTCGGCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCg                                           >  1:2849721/1‑69 (MQ=255)
   tCCGCCTTGATTTCGGCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGAtt                                        <  1:2335543/69‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:2060887/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:1930140/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:2109767/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:572931/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:28441/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:368265/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:2909474/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:3053076/68‑1 (MQ=255)
                 ggCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTACGATTACCACGAAGTTCg                           <  1:1732709/68‑1 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCAGCaaa                    >  1:1046216/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:751829/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:2905837/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:2855160/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1433201/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:3193268/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1870372/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1554667/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1067985/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:156746/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1586976/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1901301/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:263716/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1663747/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaa                     >  1:2959956/1‑68 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaa                     >  1:2415748/1‑68 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCaa                     >  1:2820800/1‑68 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTAACACGAAGTTCGCACCaaa                    >  1:1611605/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:857992/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:951201/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:2135769/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:2230075/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:581497/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:2526875/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:2768301/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:413827/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:306104/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1415591/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1406191/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1633914/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1225514/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1348151/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaaa                    >  1:1300978/1‑69 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaa                     >  1:1502799/1‑68 (MQ=255)
                       cTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCaa                     >  1:1739018/1‑68 (MQ=255)
                        tGGGGATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:1826351/61‑1 (MQ=255)
                          gggATAAGCGGCATCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:2508182/59‑1 (MQ=255)
                               aaGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCat            <  1:3155744/69‑1 (MQ=255)
                                      aTCGATAGCGCCGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCCCCAAACGCGCCATCATTgag     >  1:1335504/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCCCGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:2445263/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCCCGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:271144/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:3059059/1‑62 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:3247445/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:1029183/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:3400016/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:3025125/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:548211/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:2949189/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:603330/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:716761/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:721953/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:2842751/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:1050736/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:192886/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:2765860/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:1845897/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:2165278/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:1794231/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:2608541/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTgag     >  1:1667004/1‑69 (MQ=255)
                                      aTCGATAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATCATTga      >  1:200873/1‑68 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGaa   >  1:1444915/1‑68 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:1555583/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:929723/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:2605157/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:238851/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:236827/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:123679/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:1439528/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAc  >  1:1237874/1‑69 (MQ=255)
                                         gATAGCGACGGCTCGCTGCTGGAGTTCGATGACCACGAAGTTCGCACCAAACGCGCCATTATTAAGAAc  >  1:719037/1‑69 (MQ=255)
                                           tAGCGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCg                           <  1:2094058/42‑1 (MQ=255)
                                              cGACGGCTCGCTACTTGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:2206359/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGTGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:686685/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:1075072/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:793714/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:912938/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:2856019/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:227647/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:2680681/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:3007769/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:1353924/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:3387607/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:3264861/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:2497128/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:2501163/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:1344350/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCatca          >  1:2558441/1‑54 (MQ=255)
                                              cGACGGCTCGCTACTGGAGTTCGATTACCACGAAGTTCGACCAAACGCGCCatca          >  1:1187579/1‑53 (MQ=255)
                                                          |                                                   
AGTTCCGCCTTGATTTCGGCATTCTGGGGATAAGCGGCATCGATAGCGACGGCTCGCTGCTGGAGTTCGATTACCACGAAGTTCGCACCAAACGCGCCATTATTGAGAAC  >  minE/2660974‑2661083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: