| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,627,110 | C→T | 37.4% | R235R (CGC→CGT) | 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,110 | 0 | C | T | 37.4% | 24.2 / 22.0 | 40 | R235R (CGC→CGT) | rpoH | RNA polymerase, sigma 32 (sigma H) factor |
| Reads supporting (aligned to +/- strand): ref base C (15/10); new base T (9/6); total (24/16) | |||||||||||
| 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.87e-01 | |||||||||||
GAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGAACT > minE/2627042‑2627172 | gAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGc > 1:1086681/1‑69 (MQ=255)gAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGc > 1:1044407/1‑69 (MQ=255) aGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCa < 1:1601106/68‑1 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:910011/1‑69 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:353816/1‑69 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:710920/1‑69 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:3000510/1‑69 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:2322959/1‑69 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:231656/1‑69 (MQ=255) gatAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAg > 1:1439100/1‑69 (MQ=255) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:820287/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:3425612/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:722421/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:781461/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:2703424/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:1649677/36‑1 (MQ=16) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:91287/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2438659/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:219576/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:1401763/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2820685/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:301195/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:3115050/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:1100802/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2231095/69‑1 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTg > 1:2247009/1‑63 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:513453/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:423151/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:399328/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:3405870/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:1293498/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:939131/1‑68 (MQ=255) cGAGCGTAGTCAGGATATCATCCGCGCGCGCTGGCTGGATGAAGACAACAAGTCCACGTTGCAGGAAc > 1:3386689/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) | GAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGAACT > minE/2627042‑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: ‑ |