| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,235,894 | G→A | 67.8% | V366V (GTG→GTA) | purT → | phosphoribosylglycinamide formyltransferase 2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,235,894 | 0 | G | A | 67.8% | 20.8 / 12.9 | 31 | V366V (GTG→GTA) | purT | phosphoribosylglycinamide formyltransferase 2 |
| Reads supporting (aligned to +/- strand): ref base G (5/5); new base A (8/13); total (13/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.01e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCG > minE/1235827‑1235952 | tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTgg < 1:1608106/69‑1 (MQ=255)tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTgg < 1:2459365/69‑1 (MQ=255) ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa > 1:2771870/1‑69 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:722345/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:1505185/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:3322165/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:2831402/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:2855697/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:2224595/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCagag > 1:2916004/1‑67 (MQ=255) tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCaga > 1:2541436/1‑66 (MQ=255) tttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAgtg < 1:1063167/68‑1 (MQ=255) ttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAgtg < 1:2991621/67‑1 (MQ=255) ggaaGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:2818692/66‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:2990378/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:1248839/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:3137530/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:791359/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:868945/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:2231949/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:2167285/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:1687425/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:1551258/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:1489944/68‑1 (MQ=255) gTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:1274558/68‑1 (MQ=255) ttGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcg < 1:1734144/69‑1 (MQ=255) ttGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcg < 1:3082580/69‑1 (MQ=255) ttGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcg < 1:950459/69‑1 (MQ=255) gtcgtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTg > 1:2962040/1‑51 (MQ=255) gtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGAc > 1:1413345/1‑41 (MQ=255) gtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcgccgctg > 1:2454291/1‑67 (MQ=255) gtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGccgccg > 1:2264490/1‑69 (MQ=255) cTGGGGGTAGCACTGGCTACTGCAGAGAGTGTTGTTGAc < 1:2848056/39‑1 (MQ=255) | TAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCG > minE/1235827‑1235952 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |