| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,858,283 | A→G | 38.6% | A433A (GCT→GCC) | cysI ← | sulfite reductase, beta subunit, NAD(P)‑binding |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,858,283 | 0 | A | G | 38.6% | 17.7 / 26.7 | 44 | A433A (GCT→GCC) | cysI | sulfite reductase, beta subunit, NAD(P)‑binding |
| Reads supporting (aligned to +/- strand): ref base A (13/14); new base G (7/10); total (20/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.60e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01 | |||||||||||
TAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCG > minE/1858219‑1858340 | tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCTAGTCGGGAATGACACGCAAGcc > 1:2490360/1‑68 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:1137864/1‑69 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:477983/1‑69 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:1372485/1‑69 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:2196522/1‑69 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:1733425/1‑69 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:1714941/1‑69 (MQ=255)tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa > 1:1719343/1‑69 (MQ=255) aaaGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCaa < 1:25406/63‑1 (MQ=255) cGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGtt < 1:2164409/69‑1 (MQ=255) cGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGtt < 1:663837/69‑1 (MQ=255) ggAAACGCTCCGCTTCCGGCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:2182302/1‑69 (MQ=255) ggAAACGCTCCGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:1794268/1‑69 (MQ=255) ggAAACGCTCCGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:863203/1‑69 (MQ=255) ggAAACGCTCCGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:271950/1‑69 (MQ=255) ggAAACGCTCCGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:1880694/1‑69 (MQ=255) ggAAACGCTCCGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:1804275/1‑69 (MQ=255) ggAAACGCTCCGCTTACGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCAc > 1:2087980/1‑69 (MQ=255) cTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt < 1:1819649/67‑1 (MQ=255) cTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt < 1:2724104/67‑1 (MQ=255) cTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt < 1:1258867/67‑1 (MQ=255) tCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGAc > 1:584772/1‑68 (MQ=255) tCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGAc > 1:2126085/1‑68 (MQ=255) tCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGAc > 1:556367/1‑68 (MQ=255) tCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGAc > 1:253790/1‑68 (MQ=255) aTCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattc < 1:3338301/69‑1 (MQ=255) tCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattc < 1:2066682/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCCTTa < 1:1610000/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:2350290/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:2223308/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:1540757/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:2997480/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:781019/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:448516/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:1228562/68‑1 (MQ=255) cAGCGGGCAAGTCGGGAATGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:2104003/68‑1 (MQ=255) aaGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGctct > 1:1145654/1‑69 (MQ=255) aaTGACACGCAGGCCATTGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa < 1:2130482/52‑1 (MQ=39) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:1912655/69‑1 (MQ=255) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:1509064/69‑1 (MQ=255) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:650841/69‑1 (MQ=255) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:436609/69‑1 (MQ=255) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:362741/69‑1 (MQ=255) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:1063068/69‑1 (MQ=255) aaTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCg < 1:3075014/69‑1 (MQ=255) | TAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCG > minE/1858219‑1858340 |
| 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: ‑ |