Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,440,065 G→T 48.3% A193A (GCC→GCA ilvC ← ketol‑acid reductoisomerase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,440,0650GT48.3% 9.1 / 84.2 78A193A (GCC→GCAilvCketol‑acid reductoisomerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (14/26);  new base T (20/17);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

TTTCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGG  >  minE/2439999‑2440130
                                                                  |                                                                 
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 >  1:91032/1‑69 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 >  1:3212241/1‑69 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 <  1:284641/69‑1 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 <  1:2108704/69‑1 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 >  1:2057751/1‑69 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 <  1:1592315/69‑1 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 <  1:142451/69‑1 (MQ=255)
tttCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGc                                                                 >  1:755028/1‑69 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:3223228/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:1354834/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:2134827/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:2363101/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:2459756/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:3279006/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:500592/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:1570741/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:2908958/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:624166/69‑1 (MQ=255)
            aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:1458692/69‑1 (MQ=255)
             aCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:2315423/68‑1 (MQ=255)
               gAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCAt                                                  <  1:764380/69‑1 (MQ=255)
               gAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCAt                                                  <  1:2359557/69‑1 (MQ=255)
               gAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCAt                                                  <  1:863405/69‑1 (MQ=255)
                              cacaCCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     <  1:347918/51‑1 (MQ=255)
                                  cccGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATcac                                                     >  1:1052004/1‑45 (MQ=255)
                                  cccGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     >  1:1240751/1‑47 (MQ=255)
                                  cccGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                     >  1:1605963/1‑47 (MQ=255)
                                          gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                       >  1:365460/1‑69 (MQ=255)
                                          gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                       >  1:694804/1‑69 (MQ=255)
                                          gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                       >  1:2085448/1‑69 (MQ=255)
                                          gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                       >  1:3432359/1‑69 (MQ=255)
                                             accaccGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTg                    <  1:2298248/69‑1 (MQ=255)
                                             accaccGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTg                    <  1:2903982/69‑1 (MQ=255)
                                             accaccGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTg                    <  1:366580/69‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGTATCGTTTTcc                         <  1:217738/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:3361443/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:3168866/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:3036141/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2915635/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:167913/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:1747582/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2790277/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:273202/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2588084/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2556854/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:1025155/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:1922125/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2446943/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2083932/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:1104670/55‑1 (MQ=255)
                                                      tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                         <  1:2319254/55‑1 (MQ=255)
                                                      tGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCa             >  1:2590146/1‑67 (MQ=255)
                                                      tGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCa             >  1:2552856/1‑67 (MQ=255)
                                                      tGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCa             >  1:1262030/1‑67 (MQ=255)
                                                           ccgccCATGCTTTAGCTATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1890051/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGCTCGCCATGCCTTCGCCTTTAGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:250339/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCa             >  1:1645262/1‑62 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:2806572/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1047305/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:909501/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:874156/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1251398/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:828330/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1354105/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1409718/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:700080/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:411462/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1865606/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1964393/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1985826/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:2127307/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:25988/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:2768774/1‑69 (MQ=255)
                                                           ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:2897450/1‑69 (MQ=255)
                                                           ccgccCAAGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg      >  1:1566953/1‑69 (MQ=255)
                                                               ccAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCgg  <  1:6598/69‑1 (MQ=255)
                                                               ccAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCgg  <  1:1436140/69‑1 (MQ=255)
                                                               ccAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCgg  <  1:3340005/69‑1 (MQ=255)
                                                                  |                                                                 
TTTCACTTCCGCAACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGG  >  minE/2439999‑2440130

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: