| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 327,386 | C→T | 31.4% | S288F (TCC→TTC) | hemH → | ferrochelatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 327,386 | 0 | C | T | 31.4% | 44.8 / 25.4 | 35 | S288F (TCC→TTC) | hemH | ferrochelatase |
| Reads supporting (aligned to +/- strand): ref base C (7/17); new base T (6/5); total (13/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.58e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
ATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGA > minE/327316‑327452 | aTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGtctc > 1:3286475/1‑71 (MQ=255) gtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTccc > 1:2696635/1‑71 (MQ=255) tgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTccc < 1:3984159/70‑1 (MQ=255) tgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCt < 1:3158107/71‑1 (MQ=255) cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTTGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTc < 1:561288/70‑1 (MQ=255) cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCg < 1:3508492/71‑1 (MQ=255) cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCg < 1:1504084/71‑1 (MQ=255) cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCg < 1:2229476/71‑1 (MQ=255) gggCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGt < 1:3759384/70‑1 (MQ=255) tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggcg > 1:1364823/1‑71 (MQ=255) tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggcg > 1:2359589/1‑71 (MQ=255) tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCcggc < 1:1969683/70‑1 (MQ=255) tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggc < 1:2246116/67‑1 (MQ=255) gagaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAata < 1:3905056/71‑1 (MQ=255) agaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAatat < 1:2053202/71‑1 (MQ=255) agaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAatat < 1:3641510/71‑1 (MQ=255) aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatata < 1:632896/70‑1 (MQ=255) aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAatatat < 1:3525931/71‑1 (MQ=255) gAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgc > 1:1116514/1‑71 (MQ=255) gAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAata < 1:1160236/62‑1 (MQ=255) agaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAata > 1:2701764/1‑58 (MQ=255) agaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAata > 1:4002922/1‑58 (MQ=255) agaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCtt < 1:334994/71‑1 (MQ=255) gaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGaaaaa < 1:385849/48‑1 (MQ=255) ttGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAt > 1:1398696/1‑71 (MQ=255) tGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATg > 1:2388395/1‑71 (MQ=255) tGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATg > 1:1831623/1‑71 (MQ=255) tGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGa > 1:1683775/1‑51 (MQ=255) aGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAata < 1:4097318/49‑1 (MQ=255) gCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAata > 1:1276015/1‑42 (MQ=255) gCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGc < 1:3630080/71‑1 (MQ=255) cAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGcc > 1:29166/1‑71 (MQ=255) aaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCgg < 1:3870204/71‑1 (MQ=255) aaaCCGTGAGGTCTCCCTCGGTGCCGGCGGGaaaaaa < 1:657563/37‑1 (MQ=255) aaaCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCg > 1:92389/1‑70 (MQ=255) aaaCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCg > 1:3319611/1‑70 (MQ=255) aaCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGa < 1:1415788/71‑1 (MQ=255) gTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGa < 1:1264050/71‑1 (MQ=255) | ATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGA > minE/327316‑327452 |
| 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: ‑ |