Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,053,813 T→C 32.7% L6L (TTA→CTA) ‡ ydhM → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,053,8130TC32.7% 41.8 / 29.2 58L6L (TTA→CTA) ‡ydhMpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (15/4);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 5.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGC  >  minE/1053747‑1053881
                                                                  |                                                                    
cATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAAtt                                                                     >  1:444637/1‑68 (MQ=255)
cATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAAtt                                                                     >  1:2062890/1‑68 (MQ=255)
cATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAAtt                                                                     >  1:1375225/1‑68 (MQ=255)
   ctcCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACt                                                                  >  1:1212965/1‑68 (MQ=255)
           gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:1008473/68‑1 (MQ=255)
           gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:2462003/68‑1 (MQ=255)
           gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:3424511/68‑1 (MQ=255)
           gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCg                                                          <  1:2370844/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:1103375/1‑69 (MQ=255)
                           tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:219108/1‑69 (MQ=255)
                           tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:888620/1‑69 (MQ=255)
                           tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAggg                                         >  1:469612/1‑69 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:1874319/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:1199054/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:1311416/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:2454001/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:256294/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:2717614/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:2833049/68‑1 (MQ=255)
                             ccTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGt                                        <  1:3320436/68‑1 (MQ=255)
                              cTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAATGGt                                        <  1:3407114/67‑1 (MQ=255)
                                 cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:610911/1‑69 (MQ=255)
                                 cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:48351/1‑69 (MQ=255)
                                 cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:906531/1‑69 (MQ=255)
                                 cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:3419975/1‑69 (MQ=255)
                                 cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:580387/1‑69 (MQ=255)
                                 cAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTc                                   >  1:1771595/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2900470/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1011245/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:901675/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1267001/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:12696/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1296514/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:534519/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:1955296/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2345846/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3431263/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2619517/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:2835458/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3238095/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3042078/1‑69 (MQ=255)
                                    cGTGGATTCACCGGGATGGGGCTAAGCGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAt                                >  1:3077225/1‑69 (MQ=255)
                                           tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:2631105/1‑69 (MQ=255)
                                           tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:2946531/1‑69 (MQ=255)
                                           tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:3008817/1‑69 (MQ=255)
                                           tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:672780/1‑69 (MQ=255)
                                           tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:3165146/1‑69 (MQ=255)
                                           tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCactact                         >  1:3384086/1‑69 (MQ=255)
                                                              cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:1391079/67‑1 (MQ=255)
                                                              cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:3170225/67‑1 (MQ=255)
                                                              cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:319838/67‑1 (MQ=255)
                                                              cGAACTCCTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGTTCTAAAGAAGCg        <  1:2867968/67‑1 (MQ=255)
                                                                 aTTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGGTCTAAAGAAGCGTTTgg   <  1:1623495/69‑1 (MQ=255)
                                                                  ttACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc  <  1:2421808/69‑1 (MQ=255)
                                                                  ttACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc  <  1:1194615/69‑1 (MQ=255)
                                                                  ttACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc  <  1:2715279/69‑1 (MQ=255)
                                                                  |                                                                    
CATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGC  >  minE/1053747‑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: