Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,440,792 T→A 43.7% intergenic (‑149/‑1) ilvC ← / → ilvY ketol‑acid reductoisomerase, NAD(P)‑binding/DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,440,7920TA43.7% 21.1 / 67.5 64intergenic (‑149/‑1)ilvC/ilvYketol‑acid reductoisomerase, NAD(P)‑binding/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base A (19/9);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 5.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

ATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGCCGCAGCGCG  >  minE/2440738‑2440858
                                                      |                                                                  
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGATGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:1067610/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcta                                                      <  1:2848497/69‑3 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:701379/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:2249599/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:2305615/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:2052505/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:3051439/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:2777371/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:2554182/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:268196/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:527617/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:1414745/68‑1 (MQ=255)
atatTCACAACGTCACATTGCAATTTTTGCAACGTAAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                                       <  1:534321/68‑1 (MQ=255)
   ttCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACTCGAt                                                     >  1:3275744/1‑67 (MQ=255)
   ttCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATc                                                    >  1:2897323/1‑68 (MQ=255)
        aaCGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                                              <  1:3085975/69‑1 (MQ=255)
        aaCGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                                              <  1:424284/69‑1 (MQ=255)
        aaCGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                                              <  1:1388878/69‑1 (MQ=255)
        aaCGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                                              <  1:3274201/69‑1 (MQ=255)
         aCGTCACATTGCAATTTTTGCAACGTCAACATCTAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                                              <  1:3058280/68‑1 (MQ=255)
                  tGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCa                                    <  1:986648/69‑1 (MQ=255)
                   gCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCa                                    <  1:1080302/68‑1 (MQ=255)
                   gCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCAt                                   >  1:1921154/1‑69 (MQ=255)
                      aTTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTg                                <  1:444953/69‑1 (MQ=255)
                      aTTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTg                                <  1:1704924/69‑1 (MQ=255)
                      aTTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTg                                <  1:924086/69‑1 (MQ=255)
                      aTTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTg                                <  1:1536382/69‑1 (MQ=255)
                      aTTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTg                                <  1:498594/69‑1 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGAt                                                     >  1:891146/1‑41 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATc                                                    >  1:1798189/1‑42 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:2945288/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:3382536/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:3103045/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:376058/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:881288/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:2493854/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:2354022/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:2344025/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:923339/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:200558/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1858880/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1814280/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1726577/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1664822/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1486522/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1221992/1‑43 (MQ=255)
                             cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                                                   >  1:1212658/1‑43 (MQ=255)
                               aCGTCAACATTAAGGTCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:434076/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCTGAAAgc                       <  1:3404/69‑1 (MQ=255)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:2705439/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:3211018/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:3057633/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:2701718/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:2649123/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:2149214/69‑1 (MQ=25)
                               aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc                       <  1:1544860/69‑1 (MQ=25)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:307134/1‑67 (MQ=255)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:2897739/1‑67 (MQ=255)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:2777024/1‑67 (MQ=255)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:1016497/1‑67 (MQ=255)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:2548672/1‑67 (MQ=255)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:2044088/1‑67 (MQ=255)
                                             ggCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGc           >  1:1046095/1‑67 (MQ=255)
                                                    ccTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGCCGCAgcgcg  >  1:3255661/1‑69 (MQ=255)
                                                      |                                                                  
ATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGCCGCCATTTTGGCCGCAGCGCG  >  minE/2440738‑2440858

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: