Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,385,964 A→T 13.4% L56L (CTA→CTT metR → DNA‑binding transcriptional activator, homocysteine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,385,9640AT13.4% 158.2 / 14.1 59L56L (CTA→CTTmetRDNA‑binding transcriptional activator, homocysteine‑binding
Reads supporting (aligned to +/- strand):  ref base A (31/20);  new base T (4/4);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTCACCAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAACCAGGTACTGCCGCAAATTAGC  >  minE/2385895‑2386027
                                                                     |                                                               
gTCTCACCAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:4175068/1‑70 (MQ=255)
gTCTCACCAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTAc                                                                >  1:3608811/1‑71 (MQ=255)
     aCCAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCTGCTATTTGTGCGTAAGAGCCAGCCGCTACGCt                                                             <  1:4692691/69‑1 (MQ=255)
     aCCAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGc                                                              <  1:4421151/68‑1 (MQ=255)
      ccAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          <  1:242924/71‑1 (MQ=255)
      ccAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          <  1:720770/71‑1 (MQ=255)
       cAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCtt                                                            <  1:4709587/68‑1 (MQ=255)
        aGTTTAGCGTTCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:642158/1‑62 (MQ=255)
        aGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:2843290/1‑62 (MQ=255)
        aGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:2460218/1‑62 (MQ=255)
        aGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:1175499/1‑62 (MQ=255)
        aGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:3173875/1‑62 (MQ=255)
        aGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:4315913/1‑62 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:1700211/1‑66 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:4649305/1‑66 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:438283/1‑66 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:1346449/1‑66 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:4364045/1‑66 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:4889021/1‑66 (MQ=255)
           ttAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:1781117/1‑66 (MQ=255)
                    tGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGgaga                                             >  1:3104413/1‑70 (MQ=255)
                      gAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGga                                               <  1:1846477/66‑1 (MQ=255)
                      gAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGga                                               <  1:3422439/66‑1 (MQ=255)
                      gAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGga                                               <  1:2987568/66‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCtt                                                            <  1:4688231/49‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCtt                                                            <  1:783167/49‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCtt                                                            <  1:1649464/49‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCtt                                                            <  1:613820/49‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCtt                                                            <  1:178418/49‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAAAGCCAGCCGCTTCGGTTTACTCCGCAGGGAGAAATCCTg                                      <  1:4448841/71‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAAAGCCAGCCGCTTCGCTTTACACCGCAGGGAGAAATCCTg                                      <  1:780440/71‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAAAGCCAGCCGCTTCGCTTTACACCGCAGGGAGAAATCCTg                                      <  1:578227/71‑1 (MQ=255)
                          aacGCCTTGGCTTCCGGCTATTTGTGCGTAAAAGCCAGCCGCTTCGCTTTACACCGCAGGGAGAAATCCTg                                      <  1:1875468/71‑1 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:108138/1‑50 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:5110326/1‑50 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:3859143/1‑50 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:950770/1‑50 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:2046033/1‑50 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:1948650/1‑50 (MQ=255)
                           acGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:2030769/1‑50 (MQ=255)
                            cGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTa                                                                 >  1:4228431/1‑42 (MQ=255)
                            cGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:3357442/1‑49 (MQ=255)
                            cGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTa                                                          >  1:3102332/1‑49 (MQ=255)
                             gCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGc                                                    <  1:2733181/54‑1 (MQ=255)
                               cTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCa                                 <  1:2199478/71‑1 (MQ=255)
                                ttGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGgaga                                             >  1:1940155/1‑58 (MQ=255)
                                  ggCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGtt                                    >  1:3912613/1‑65 (MQ=255)
                                  ggCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGtt                                    >  1:2326966/1‑65 (MQ=255)
                                  ggCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGtt                                    >  1:4546835/1‑65 (MQ=255)
                                  ggCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               <  1:1313580/70‑1 (MQ=255)
                                  ggCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCAACt                              >  1:2409197/1‑71 (MQ=255)
                                  ggCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCAACt                              >  1:47720/1‑71 (MQ=255)
                                                               gccGCTTCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAATCAGGTGCTGCCGCAAATTAGc  >  1:4392339/1‑70 (MQ=255)
                                                               gccGCTTCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAATCAGGTGCTGCCGCAAATTAGc  >  1:3062036/1‑70 (MQ=255)
                                                               gccGCTTCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAATCAGGTGCTGCCGCAAATTAGc  >  1:830278/1‑70 (MQ=255)
                                                               gccGCTTCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAATCAGGTGCTGCCGCAAATTAGc  >  1:2606519/1‑70 (MQ=255)
                                                                  gCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCa                                 <  1:2788119/36‑1 (MQ=255)
                                                                  gCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCa                                 <  1:4238889/36‑1 (MQ=255)
                                                                  gCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCa                                 <  1:4178148/36‑1 (MQ=255)
                                                                  gCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCa                                 <  1:2668242/36‑1 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:1368081/1‑36 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:1248060/1‑36 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:3648176/1‑36 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:443357/1‑36 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:3152684/1‑36 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:1004273/1‑36 (MQ=255)
                                                                    tACGCTTTACACCGCAGGGAGAAATCCTGTTGCAAc                               >  1:995511/1‑36 (MQ=255)
                                                                     aCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAAc                       >  1:813540/1‑43 (MQ=255)
                                                                     aCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAAc                       >  1:459645/1‑43 (MQ=255)
                                                                     aCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAAc                       >  1:407619/1‑43 (MQ=255)
                                                                     aCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAAc                       >  1:4062615/1‑43 (MQ=255)
                                                                     aCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAAc                       >  1:3918116/1‑43 (MQ=255)
                                                                     aCGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAAc                       >  1:3487087/1‑43 (MQ=255)
                                                                     |                                                               
GTCTCACCAGTTTAGCGATCTGGAACAACGCCTTGGCTTCCGGCTATTTGTGCGTAAGAGCCAGCCGCTACGCTTTACACCGCAGGGAGAAATCCTGTTGCAACTGGCAAACCAGGTACTGCCGCAAATTAGC  >  minE/2385895‑2386027

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: