| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,898,833 | C→T | 60.7% | Q222Q (CAG→CAA) | recD ← | exonuclease V (RecBCD complex), alpha chain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,898,833 | 0 | C | T | 60.7% | 25.2 / 40.4 | 61 | Q222Q (CAG→CAA) | recD | exonuclease V (RecBCD complex), alpha chain |
| Reads supporting (aligned to +/- strand): ref base C (14/10); new base T (28/9); total (42/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 | |||||||||||
GGTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCA > minE/1898775‑1898893 | ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg < 1:2753446/69‑1 (MQ=255)ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg < 1:744643/69‑1 (MQ=255)ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg < 1:2698583/69‑1 (MQ=255)ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg < 1:2123207/69‑1 (MQ=255) tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCTCAAAGc < 1:911/68‑1 (MQ=255) tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc < 1:3274919/68‑1 (MQ=255) tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc < 1:2954325/68‑1 (MQ=255) tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc < 1:668408/68‑1 (MQ=255) tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc < 1:660176/68‑1 (MQ=255) tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc < 1:174238/68‑1 (MQ=255) gTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTg > 1:1821496/1‑66 (MQ=255) cTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGcc > 1:1705326/1‑65 (MQ=255) aTCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTc > 1:1436245/1‑69 (MQ=255) tcctCCGGAATGCGTTTCTTTTGTTCAT‑GGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2784127/4‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTGAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:2106968/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:2159918/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:219899/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:2556396/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:593703/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:35772/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:1766312/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:1874844/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc > 1:783627/1‑68 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2564719/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:107544/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:3057903/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2974366/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:115358/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2897097/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:1165060/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:1555350/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:1943514/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:1058904/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2437037/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2435756/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:2387649/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:845494/1‑69 (MQ=255) tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGAAATTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:3080494/1‑69 (MQ=255) ttCCGGAATGCGTTTCTTTTGTTCATCGGTTAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCgg > 1:135571/1‑68 (MQ=255) ttCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCgg > 1:2365755/1‑68 (MQ=255) ttCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCg > 1:658260/1‑67 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCgagag > 1:2120042/1‑48 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTaa > 1:2168670/1‑58 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg > 1:795182/1‑69 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg > 1:2443657/1‑69 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg > 1:2266127/1‑69 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg > 1:2230096/1‑69 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg > 1:1966511/1‑69 (MQ=255) tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg > 1:1141885/1‑69 (MQ=255) tCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGtt > 1:257100/1‑69 (MQ=255) tCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGtt > 1:213213/1‑69 (MQ=255) tCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGtt > 1:914267/1‑69 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:3423337/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:271726/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:509073/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:2704542/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:2314643/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:714867/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:718865/69‑1 (MQ=255) aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:2044983/69‑1 (MQ=255) gCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa < 1:3380302/68‑1 (MQ=255) | GGTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCA > minE/1898775‑1898893 |
| 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: ‑ |