Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,440,780 G→A 41.6% intergenic (‑137/‑13) 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,7800GA41.6% 25.2 / 47.0 68intergenic (‑137/‑13)ilvC/ilvYketol‑acid reductoisomerase, NAD(P)‑binding/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (8/31);  new base A (19/9);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 1.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AAATTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGC  >  minE/2440717‑2440837
                                                               |                                                         
aaaTTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGc                                                      <  1:9167/69‑1 (MQ=255)
aaaTTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGc                                                      <  1:1932812/69‑1 (MQ=255)
aaaTTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGc                                                      <  1:2774790/69‑1 (MQ=255)
aaaTTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGc                                                      <  1:1688405/69‑1 (MQ=255)
aaaTTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGc                                                      <  1:1562806/69‑1 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  >  1:1637257/1‑69 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  >  1:1922940/1‑69 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  <  1:1861678/69‑1 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  >  1:3240125/1‑69 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  >  1:2276420/1‑69 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  >  1:2380069/1‑69 (MQ=255)
    ttATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTc                                                  <  1:71414/69‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGATGGCTGTCCCTGTGGATTTAcgcg                                  <  1:1067610/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcta                                 <  1:2848497/69‑3 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:2554182/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:2305615/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:2777371/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:2249599/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:2052505/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:268196/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:3051439/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:527617/68‑1 (MQ=255)
                     atatTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTAcgcg                                  <  1:701379/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:3274201/69‑1 (MQ=255)
                             aaCGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                         <  1:424284/69‑1 (MQ=255)
                             aaCGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGaaaa                         <  1:1388878/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: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)
                                           aTTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTg           <  1:444953/69‑1 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGAt                                >  1:891146/1‑41 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATc                               >  1:1798189/1‑42 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1664822/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:2493854/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:923339/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1212658/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:881288/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1221992/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1486522/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:2354022/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:2344025/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:200558/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:376058/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:3382536/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1726577/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:3103045/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1814280/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:1858880/1‑43 (MQ=255)
                                                  cAACGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCt                              >  1:2945288/1‑43 (MQ=255)
                                                    aCGTCAACATTAAGGTCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:434076/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCTGAAAgc  <  1:3404/69‑1 (MQ=255)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:2649123/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:1544860/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:2705439/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:3211018/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:3057633/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:2149214/69‑1 (MQ=25)
                                                    aCGTCAACATTAAGGGCTGCGCCATTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAgc  <  1:2701718/69‑1 (MQ=25)
                                                               |                                                         
AAATTTATTGCGGAAATTGATATATTCACAACGTCACATTGCAATTTTTGCAACGTCAACATCGAGGGCTGTCCCTGTGGATTTACGCGATCTGAAAACCTTCCTGCATCTGGCGGAAAGC  >  minE/2440717‑2440837

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