| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 241,523 | C→T | 61.6% | Q900Q (CAG→CAA) | sbcC ← | exonuclease, dsDNA, ATP‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 241,523 | 0 | C | T | 61.6% | 15.7 / 22.8 | 39 | Q900Q (CAG→CAA) | sbcC | exonuclease, dsDNA, ATP‑dependent |
| Reads supporting (aligned to +/- strand): ref base C (6/9); new base T (16/8); total (22/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.84e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01 | |||||||||||
CAGCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGAATTCAGATATC > minE/241469‑241581 | cAGCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGc > 1:2334302/1‑68 (MQ=255)cAGCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGc > 1:2867380/1‑68 (MQ=255)cAGCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGc > 1:2612896/1‑68 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:713878/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:3143835/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:152294/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:1525843/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:2868670/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:1557342/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:1560688/69‑1 (MQ=255) gCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCTTGGGCAAACTTGCGGa < 1:1662754/69‑1 (MQ=255) ttGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGcc < 1:2881282/69‑1 (MQ=255) aGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGa < 1:1566493/69‑1 (MQ=255) aGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGa > 1:2126591/1‑69 (MQ=255) cGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGaa < 1:1882757/68‑1 (MQ=255) cGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGaa < 1:2184157/68‑1 (MQ=255) cGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGaa < 1:1448257/68‑1 (MQ=255) agaTGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCt < 1:2683016/69‑1 (MQ=255) cTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCg < 1:2594832/69‑1 (MQ=255) cTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCg < 1:555981/69‑1 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTa > 1:2617316/1‑39 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:2628746/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:3211762/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:457167/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:2502711/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:2472504/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:2353652/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:1814016/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:1804281/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:826890/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:1538622/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:153568/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:1505207/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGaa > 1:1117845/1‑66 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGa > 1:964918/1‑65 (MQ=255) ttATCCAGCGTCAGCCCTTGGGCAAACTTGCGGAATTTATCGCCCTCTTCGGAACCTATTAGCGaa > 1:2609072/1‑66 (MQ=255) cAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGAATTCAGata > 1:1676636/1‑69 (MQ=255) cAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGAATTCAGata > 1:1144342/1‑69 (MQ=255) gCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGAATTCAGATATc < 1:2273868/69‑1 (MQ=255) | CAGCCGGGTAAGTTGCTGATTAGCGAGATGGACTAAATTATCCAGCGTCAGCCCCTGGGCAAACTTGCGGAATTTATCGCCCTCTTTGGAACCTATTAGCGAATTCAGATATC > minE/241469‑241581 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |