| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,627,107 | A→G | 35.7% | E234E (GAA→GAG) | 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,107 | 0 | A | G | 35.7% | 26.5 / 22.3 | 42 | E234E (GAA→GAG) | rpoH | RNA polymerase, sigma 32 (sigma H) factor |
| Reads supporting (aligned to +/- strand): ref base A (17/10); new base G (9/6); total (26/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.63e-01 | |||||||||||
ATTGAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGAACT > minE/2627039‑2627172 | aTTGAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGaa > 1:706786/1‑69 (MQ=255) 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:2703424/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:722421/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:781461/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:820287/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:3425612/36‑1 (MQ=16) cgatgcgatgCAAGGTCTGGACGAGCGTAGTCAGGa < 1:1649677/36‑1 (MQ=16) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2438659/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:91287/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2820685/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:219576/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:1100802/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:301195/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:3115050/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:1401763/69‑1 (MQ=255) gCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCAc < 1:2231095/69‑1 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTg > 1:2247009/1‑63 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:423151/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:399328/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:1293498/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:513453/1‑69 (MQ=255) gACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGaa > 1:3405870/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) | ATTGAAGATGATAACTGGGAAGAGCAGGCGGCAAACCGTCTGACCGACGCGATGCAGGGTCTGGACGAACGCAGCCAGGACATCATCCGTGCGCGCTGGCTGGACGAAGACAACAAGTCCACGTTGCAGGAACT > minE/2627039‑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: ‑ |