| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | DH1 | 236,086 | A→G | intergenic (+28/+26) | ECDH1_RS01120 → / ← gor | hypothetical protein/glutathione reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | DH1 | 236,086 | 0 | A | G | 84.6% | 26.5 / ‑0.2 | 13 | intergenic (+28/+26) | ECDH1_RS01120/gor | hypothetical protein/glutathione reductase |
| Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (5/6); total (7/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.62e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01 | |||||||||||
CGCATGGAGAAAACGGGTGATTGATAAAGCAATCATCGTTCTAGGGGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAACACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTGGTTGCCCCCATCTTCAGCGCCACCGCGAAGCCCTGCAACATTTCGTCCAT > DH1/235959‑236231 | cGCATGGAGAAAACGGGTGATTGATAAAGCAATCATCGTTCTAGGGGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTaa < 3:272039/150‑1 (MQ=255) atcgagaAAACGGGTGATTGATAAAGCAATCATCGTTCTAGGCGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACACTTAGCCCTTTAACAtt > 3:214151/4‑151 (MQ=255) gattgatAAAGCAATCATCGTTCTAGGGGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACg > 7:672691/1‑151 (MQ=255) aGCAATCATCGTTCTAGGGGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTCAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTCCTCTTAGCCCTTTAACATTAAACGCATTGTCACGAACtcttct > 3:1209273/1‑151 (MQ=255) atcGTGCTAGGGGCGTTAATTGCGCTGCTGGAACTGCTCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGccgccg < 7:792868/151‑1 (MQ=255) aGGGGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAACACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTg > 1:806188/1‑150 (MQ=255) cGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGc < 5:1112187/151‑1 (MQ=255) tCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAACACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCtttttt > 1:794283/1‑151 (MQ=255) ccGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTg < 3:461378/151‑1 (MQ=255) cTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTGGTTGCCCCCATCTTCAgcg > 2:1151712/1‑151 (MQ=255) aTTAAGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTGGTTGCCCCCATCTTCAGCGCCACCGCGAAGCCCTGCAACAt < 6:762973/151‑1 (MQ=255) aaGGGCTAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTGGTTGCCCCCATCTTCAGCGCCACCGCGAAGCCCTGCAACATTTc < 8:672329/151‑1 (MQ=255) tAAGAGCACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTGGTTGCCCCCATCTTCAGCGCCACCGCGCAGCCCTGCACCATTTCGTCCAt > 8:1071900/1‑151 (MQ=255) | CGCATGGAGAAAACGGGTGATTGATAAAGCAATCATCGTTCTAGGGGCGTTAATTGCGCTGCTGGAACTGATCCGCTTTCTGCTTCAGCTTCTGAACTGATAGCGGAAACGTAATTAAGGGCTAAGAACACACTACTCTTAGCCCTTTAACATTTAACGCATTGTCACGAACTCTTCTGCCGCCGTTGGGTGAATGGCGACGGTATTGTCGAAGTCTTTTTTGGTTGCCCCCATCTTCAGCGCCACCGCGAAGCCCTGCAACATTTCGTCCAT > DH1/235959‑236231 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |