Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,967,130 C→T 81.0% A500A (GCG→GCA tktA ← transketolase 1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,967,1300CT81.0% 66.9 / 8.4 47A500A (GCG→GCAtktAtransketolase 1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (21/17);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGC  >  minE/1967085‑1967192
                                             |                                                              
tGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACg                                          <  1:125759/68‑1 (MQ=255)
    gTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGc                                     <  1:3253154/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTTGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:1402776/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:100618/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:2450057/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:2459884/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:2797177/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:2943917/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:3183661/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:328527/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:369489/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:1615258/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:1464835/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:13890/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:72749/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:1355442/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:1133656/69‑1 (MQ=255)
                      tCAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGGCACACGGACGCCATGTAGACATGTTCggg                   <  1:1980384/69‑1 (MQ=255)
                       cAACGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   <  1:1453415/68‑1 (MQ=255)
                       cAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   >  1:51638/1‑68 (MQ=255)
                       cAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   >  1:697245/1‑68 (MQ=255)
                       cAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   >  1:3068782/1‑68 (MQ=255)
                       cAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCggg                   >  1:2326213/1‑68 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTcac                                             >  1:824840/1‑41 (MQ=37)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:287647/1‑69 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:2963549/1‑69 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:595338/1‑69 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:522313/1‑69 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:3232132/1‑69 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:1956869/1‑69 (MQ=255)
                        aaCGCCGTATTTCCACGCTACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGt                 >  1:1898502/1‑69 (MQ=255)
                           acCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACg                                          <  1:2433625/41‑1 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:1131486/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:813678/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:1385953/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:678950/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:1757655/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:430837/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:420151/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:2072890/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:3160833/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:2131289/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:2272553/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:2287582/1‑69 (MQ=255)
                                       cgctACTGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTCACGCGCAGAGAAGc  >  1:269419/1‑69 (MQ=255)
                                       cgcgACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGc  <  1:2023194/69‑1 (MQ=255)
                                       cgcgACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGc  <  1:1498204/69‑1 (MQ=255)
                                             |                                                              
TGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGC  >  minE/1967085‑1967192

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: