Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,968,732 G→T 54.6% intergenic (‑103/‑175) tktA ← / → yggG transketolase 1, thiamin‑binding/predicted peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,968,7320GT54.6% 3.8 / 98.2 86intergenic (‑103/‑175)tktA/yggGtransketolase 1, thiamin‑binding/predicted peptidase
Reads supporting (aligned to +/- strand):  ref base G (9/30);  new base T (32/15);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 3.63e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCA  >  minE/1968669‑1968798
                                                               |                                                                  
gACTTGACGACAGCGCGTTTTGGGCTACTCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTc                                                                <  1:963186/68‑1 (MQ=255)
gACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTc                                                                <  1:1880135/68‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCTGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:2554962/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:791793/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:580313/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:2612807/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:1638993/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:3002481/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:1737188/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:2856754/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:1765200/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:1866459/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:2819091/69‑1 (MQ=255)
    tGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTa                                                           <  1:1834992/69‑1 (MQ=255)
     gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTAc                                                          >  1:2823584/1‑69 (MQ=255)
     gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTAc                                                          >  1:2338340/1‑69 (MQ=255)
     gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTAc                                                          >  1:1786264/1‑69 (MQ=255)
     gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTAc                                                          >  1:496686/1‑69 (MQ=255)
     gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTAc                                                          >  1:817464/1‑69 (MQ=255)
             cgcgTTTTGGGCTACGCCGGAAAATTTGCGAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACa                                                  >  1:2805984/1‑69 (MQ=255)
             cgcgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACa                                                  >  1:2404373/1‑69 (MQ=255)
                   ttGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCttt                                            >  1:3191358/1‑69 (MQ=255)
                   ttGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCttt                                            >  1:2258616/1‑69 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACAt                                                        >  1:577543/1‑50 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATgg                                                      >  1:2471384/1‑52 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACAt                                                 >  1:976858/1‑57 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:259463/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2874901/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:1737660/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:1721459/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:1076992/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:462502/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2794500/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:293685/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:883651/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2937500/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:1660536/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2583036/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:3212841/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2515959/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2506229/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2478212/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:386989/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:1065056/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2426834/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:3416975/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2392178/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:2221622/1‑68 (MQ=255)
                          aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:3315421/1‑68 (MQ=255)
                          aCGCCGCAAAATTTGCCAACAATTTACCGCAAGCCGCTCGTCATGTACATGGAACATCCTTTTGCCGc                                      >  1:7978/1‑68 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:1134656/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:1838362/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:1388549/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:3113810/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:3117060/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:446416/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:700641/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:106509/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:849697/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:1992645/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:935504/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:2145772/39‑1 (MQ=255)
                            gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:2181900/39‑1 (MQ=255)
                             ccGGAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:1776937/38‑1 (MQ=255)
                               ggAAAATTTGCCAACAATTTACCGCAAGCCGCTCGt                                                                 <  1:3013859/36‑1 (MQ=255)
                               ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGa                                <  1:1139424/69‑1 (MQ=255)
                               ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGa                                <  1:2662925/69‑1 (MQ=255)
                               ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGa                                <  1:2658936/69‑1 (MQ=255)
                               ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGa                                <  1:2110596/69‑1 (MQ=255)
                                 aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGaaa                              <  1:1044491/69‑1 (MQ=255)
                                 aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGaaa                              <  1:398577/69‑1 (MQ=255)
                                 aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGaaa                              <  1:3411810/69‑1 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   >  1:2936787/1‑69 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   >  1:3218849/1‑69 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   <  1:3216687/69‑1 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   <  1:398760/69‑1 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   <  1:463189/69‑1 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   <  1:3045340/69‑1 (MQ=255)
                                            acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                   >  1:2757084/1‑69 (MQ=255)
                                              aaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGaaa                              <  1:872279/56‑1 (MQ=255)
                                                  tACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCg             <  1:2795013/69‑1 (MQ=255)
                                                  tACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCg             >  1:3036821/1‑69 (MQ=255)
                                                  tACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCg             >  1:2126672/1‑69 (MQ=255)
                                                  tACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCg             <  1:3190331/69‑1 (MQ=255)
                                                  tACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCg             <  1:2362057/69‑1 (MQ=255)
                                                     cGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCa           <  1:2771482/68‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:2464199/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:866945/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:299329/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:1895785/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:726181/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:2949210/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:520962/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:3113302/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:46039/69‑1 (MQ=255)
                                                             gcgcgTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa  <  1:3292458/69‑1 (MQ=255)
                                                               |                                                                  
GACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCA  >  minE/1968669‑1968798

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: