| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 598,260 | C→T | 41.0% | Y201Y (TAC→TAT) | dmsA → | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 598,260 | 0 | C | T | 41.0% | 12.8 / 24.2 | 39 | Y201Y (TAC→TAT) | dmsA | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Reads supporting (aligned to +/- strand): ref base C (7/16); new base T (8/8); total (15/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.18e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01 | |||||||||||
CCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCA > minE/598211‑598320 | cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCcgcg > 1:224691/1‑68 (MQ=255)cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCcgcg > 1:1849776/1‑68 (MQ=255)cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCcgc > 1:2168812/1‑67 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:976480/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:70364/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:2833711/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:2573647/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:1743985/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:151911/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:1425882/1‑69 (MQ=255) tGGTCGCGCGGCTAATGAACAGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa > 1:2870499/1‑69 (MQ=255) gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGGACCATTACGGCGACTACTCCTCCGCGCAAAtt > 1:1643637/1‑69 (MQ=255) gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAAtt > 1:1678160/1‑69 (MQ=255) gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAAtt > 1:2261557/1‑69 (MQ=255) gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAAACATTACGGCGACTACTCCTCCGCGCAAAtt > 1:1765870/1‑69 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:1725796/67‑1 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:788063/67‑1 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:1927173/67‑1 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:536917/67‑1 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:216511/67‑1 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:1475076/67‑1 (MQ=255) gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg < 1:2653944/67‑1 (MQ=255) aaCTGCTGCGGCGGCTATCTGCACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc < 1:1035475/69‑1 (MQ=255) aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc < 1:2520035/69‑1 (MQ=255) aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc < 1:3130460/69‑1 (MQ=255) aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc < 1:611823/69‑1 (MQ=255) aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc < 1:153515/69‑1 (MQ=255) ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:2499329/69‑1 (MQ=255) ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:1457156/69‑1 (MQ=255) ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:2332049/69‑1 (MQ=255) ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:626833/69‑1 (MQ=255) ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:1885844/69‑1 (MQ=255) ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:909878/69‑1 (MQ=255) tgctgcGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:2123543/68‑1 (MQ=255) tgctgcGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta < 1:1050960/68‑1 (MQ=255) ccATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCa < 1:2493107/67‑1 (MQ=255) ccATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCa < 1:2135890/67‑1 (MQ=255) ccATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCa < 1:579514/67‑1 (MQ=255) ccATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCa < 1:958436/67‑1 (MQ=255) | CCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCA > minE/598211‑598320 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |