| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 863,516 | G→C | 59.2% | S155S (TCG→TCC) | yciO → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 863,516 | 0 | G | C | 59.2% | 21.1 / 46.4 | 49 | S155S (TCG→TCC) | yciO | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (15/5); new base C (7/22); total (22/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.04e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01 | |||||||||||
TGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAC > minE/863450‑863570 | tGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGa < 1:1712829/69‑1 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaaga > 1:2289408/1‑69 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaaga > 1:1858419/1‑69 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaaga > 1:2642613/1‑69 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaaga > 1:2792729/1‑69 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaaga > 1:1451939/1‑69 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaaga > 1:2975685/1‑69 (MQ=255) gggCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGATCCGgaag > 1:1522780/1‑68 (MQ=255) cGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAAt < 1:2108472/69‑1 (MQ=255) cGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAAt < 1:1836805/69‑1 (MQ=255) cGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAAt < 1:1437190/69‑1 (MQ=255) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGTATCCGa < 1:435361/54‑1 (MQ=38) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:3285009/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:466335/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:2865530/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:2743101/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:549036/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:2526451/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:2377969/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:742072/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:761902/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:116368/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:787887/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1519599/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1164029/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:120289/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1256091/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1389228/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1496978/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:175715/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1760269/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1895606/54‑1 (MQ=39) aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGTGAATTTACCGAATCCGa < 1:1768986/54‑1 (MQ=39) cGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAg < 1:357499/69‑1 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:1769075/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:1826067/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:1784119/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:509249/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:1508529/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:1963385/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:3191097/1‑67 (MQ=255) tCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGaaa > 1:3095644/1‑67 (MQ=255) ctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGa > 1:1766796/1‑68 (MQ=255) ctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGAt > 1:622808/1‑69 (MQ=255) ctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGAt > 1:297511/1‑69 (MQ=255) gCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGatta > 1:574778/1‑69 (MQ=255) gCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGatt > 1:684366/1‑68 (MQ=255) gAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAc > 1:447184/1‑69 (MQ=255) gAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAc > 1:254742/1‑69 (MQ=255) | TGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAC > minE/863450‑863570 |
| 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: ‑ |