Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,053,815 A→C 31.1% L6F (TTA→TTC) ‡ ydhM → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,053,8150AC31.1% 56.6 / 40.1 61L6F (TTA→TTC) ‡ydhMpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base C (15/4);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 2.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGC  >  minE/1053750‑1053881
                                                                 |                                                                  
ctcCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACt                                                                  >  1:1212965/1‑68 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:1008473/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:2462003/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:2370844/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:3424511/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:1845240/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:1838884/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:163864/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:1582858/68‑1 (MQ=255)
        gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:89007/68‑1 (MQ=255)
           gggCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGa                                                      >  1:3402628/1‑69 (MQ=255)
                       ttGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAgg                                          >  1:627333/1‑69 (MQ=255)
                        tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:518159/1‑69 (MQ=255)
                        tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:888620/1‑69 (MQ=255)
                        tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:469612/1‑69 (MQ=255)
                        tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:219108/1‑69 (MQ=255)
                        tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:1103375/1‑69 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:2833049/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:1311416/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:2717614/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:3320436/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:256294/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:2454001/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:1874319/68‑1 (MQ=255)
                          ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:1199054/68‑1 (MQ=255)
                           cTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAATGGt                                        <  1:3407114/67‑1 (MQ=255)
                              cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:3419975/1‑69 (MQ=255)
                              cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:906531/1‑69 (MQ=255)
                              cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:48351/1‑69 (MQ=255)
                              cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:610911/1‑69 (MQ=255)
                              cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:580387/1‑69 (MQ=255)
                              cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:1771595/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1267001/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1296514/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:534519/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:12696/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1955296/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3431263/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2345846/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2619517/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3238095/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:901675/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1011245/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3077225/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3042078/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2900470/1‑69 (MQ=255)
                                 cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2835458/1‑69 (MQ=255)
                                        tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:3384086/1‑69 (MQ=255)
                                        tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:2631105/1‑69 (MQ=255)
                                        tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:672780/1‑69 (MQ=255)
                                        tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:3165146/1‑69 (MQ=255)
                                        tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:3008817/1‑69 (MQ=255)
                                        tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:2946531/1‑69 (MQ=255)
                                                           cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:2867968/67‑1 (MQ=255)
                                                           cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:1391079/67‑1 (MQ=255)
                                                           cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:319838/67‑1 (MQ=255)
                                                           cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:3170225/67‑1 (MQ=255)
                                                              aTTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGGTCTAAAGAAGCGTTTgg   <  1:1623495/69‑1 (MQ=255)
                                                               ttACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc  <  1:2421808/69‑1 (MQ=255)
                                                               ttACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc  <  1:2715279/69‑1 (MQ=255)
                                                               ttACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc  <  1:1194615/69‑1 (MQ=255)
                                                                 |                                                                  
CTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGC  >  minE/1053750‑1053881

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: