| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,627,113 | C→T | 38.4% | S236S (AGC→AGT) | rpoH → | RNA polymerase, sigma 32 (sigma H) factor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,627,113 | 0 | C | T | 38.4% | 21.5 / 21.9 | 39 | S236S (AGC→AGT) | rpoH | RNA polymerase, sigma 32 (sigma H) factor |
| Reads supporting (aligned to +/- strand): ref base C (15/9); new base T (9/6); total (24/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01 | |||||||||||
GATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGAACT > minE/2627048‑2627172 | gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:353816/1‑69 (MQ=255)gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:910011/1‑69 (MQ=255)gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:2322959/1‑69 (MQ=255)gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:231656/1‑69 (MQ=255)gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:3000510/1‑69 (MQ=255)gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:1439100/1‑69 (MQ=255)gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:710920/1‑69 (MQ=255) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:1649677/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:722421/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:3425612/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:820287/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:2703424/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:781461/36‑1 (MQ=16) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:3115050/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:301195/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:1100802/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2820685/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2438659/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:91287/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2231095/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:219576/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:1401763/69‑1 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTg > 1:2247009/1‑63 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:513453/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:3405870/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:1293498/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:399328/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:423151/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGa > 1:3253080/1‑68 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGa > 1:493650/1‑68 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGctggctgg > 1:1109095/1‑38 (MQ=12) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAAc > 1:3386689/1‑68 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAAc > 1:939131/1‑68 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAACt > 1:3310525/1‑69 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAACt > 1:2933225/1‑69 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAACt > 1:280812/1‑69 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAACt > 1:2477361/1‑69 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAACt > 1:2436709/1‑69 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAACt > 1:2305190/1‑69 (MQ=255) | GATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGAACT > minE/2627048‑2627172 |
| 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: ‑ |