Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,440,071 G→A 51.3% A191A (GCC→GCT 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,0710GA51.3% ‑4.6 / 69.2 76A191A (GCC→GCTilvCketol‑acid reductoisomerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (15/22);  new base A (21/18);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTACG  >  minE/2440011‑2440134
                                                            |                                                               
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:2363101/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:624166/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:2908958/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:500592/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:3279006/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:2134827/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:1354834/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:3223228/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:1570741/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:1458692/69‑1 (MQ=255)
aaCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:2459756/69‑1 (MQ=255)
 aCGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:2315423/68‑1 (MQ=255)
   gAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCAt                                                      <  1:764380/69‑1 (MQ=255)
   gAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCAt                                                      <  1:863405/69‑1 (MQ=255)
   gAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCAt                                                      <  1:2359557/69‑1 (MQ=255)
                  cacaCCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         <  1:347918/51‑1 (MQ=255)
                      cccGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATcac                                                         >  1:1052004/1‑45 (MQ=255)
                      cccGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         >  1:1605963/1‑47 (MQ=255)
                      cccGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGc                                                         >  1:1240751/1‑47 (MQ=255)
                              gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                           >  1:365460/1‑69 (MQ=255)
                              gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                           >  1:2085448/1‑69 (MQ=255)
                              gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                           >  1:3432359/1‑69 (MQ=255)
                              gTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCgg                           >  1:694804/1‑69 (MQ=255)
                                 accaccGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTg                        <  1:2298248/69‑1 (MQ=255)
                                 accaccGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTg                        <  1:366580/69‑1 (MQ=255)
                                 accaccGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTg                        <  1:2903982/69‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGTATCGTTTTcc                             <  1:217738/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:3361443/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2556854/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2588084/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:273202/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:3036141/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2790277/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:3168866/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2915635/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2319254/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:1025155/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:1104670/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:167913/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:1747582/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:1922125/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2083932/55‑1 (MQ=255)
                                          tGCAGCCGCCCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:2446943/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:411462/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:700080/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:828330/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:874156/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:909501/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:1964393/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:1047305/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:1251398/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:1354105/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:1409718/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:1865606/1‑69 (MQ=255)
                                               ccgccCATGCTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCg          >  1:2897450/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:2806572/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)
                                                        cTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTAct  >  1:2597829/1‑67 (MQ=255)
                                                        cTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTACg  >  1:2059647/1‑68 (MQ=255)
                                                        cTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTACg  >  1:273255/1‑68 (MQ=255)
                                                        cTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTACg  >  1:1582161/1‑68 (MQ=255)
                                                        cTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTACg  >  1:3190514/1‑68 (MQ=255)
                                                        cTTTAGCGATCGCCATGCCTTCGCCTTTCGGATCGTTTTcc                             <  1:673045/41‑1 (MQ=38)
                                                            |                                                               
AACGAAGGACGATTCCAGCACACCCGCACGGTGACCACCGGTTGCAGCCGCCCAGGCTTTGGCAATCGCCATGCCTTCGCCTTTCGGATCGTTTTCCGGGTGAACGGCAATCAGCGTCGGTACG  >  minE/2440011‑2440134

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: