| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,781,149 | T→C | 8.7% | F708F (TTT→TTC) | metH → | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,781,149 | 0 | T | C | 8.7% | 133.9 / 4.2 | 46 | F708F (TTT→TTC) | metH | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent |
| Reads supporting (aligned to +/- strand): ref base T (13/29); new base C (2/2); total (15/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.87e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01 | |||||||||||
AGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGCCTAC > minE/2781079‑2781218 | agGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGttt < 1:1363359/71‑1 (MQ=255)agGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGttt < 1:2590779/71‑1 (MQ=255) ggCTACTCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTg < 1:2212796/71‑1 (MQ=255) ggCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTg < 1:2001757/71‑1 (MQ=255) aCGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCg > 1:3354002/1‑70 (MQ=255) gcgcCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAg > 1:2644149/1‑71 (MQ=255) gcCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAggg < 1:2483070/71‑1 (MQ=255) ccGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaa < 1:1603137/71‑1 (MQ=255) gATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaaaa > 1:3321140/1‑71 (MQ=255) aTTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaaaa > 1:1372411/1‑70 (MQ=255) aTTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAAt > 1:2516836/1‑71 (MQ=255) ttGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATg < 1:2256499/71‑1 (MQ=255) aaGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGa < 1:1301317/58‑1 (MQ=255) aaGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGa < 1:49677/58‑1 (MQ=255) aaGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGtt < 1:2777277/70‑1 (MQ=255) aaGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGtt < 1:2783293/70‑1 (MQ=255) gATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGttt < 1:2017253/49‑1 (MQ=255) atcgAAGGTCCGCTGATGGACGGCATGAACATAGTCGGCGACCTGTTCGGCCAGGGCAAGATGTTCCTGc > 1:884485/4‑70 (MQ=255) aTTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGc < 1:875379/70‑1 (MQ=255) ttGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGttt < 1:435021/47‑1 (MQ=255) ttGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGttt < 1:139670/47‑1 (MQ=255) ttGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGttt < 1:1313610/47‑1 (MQ=255) ttGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCca < 1:1020590/71‑1 (MQ=255) aaGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCcaca < 1:3406984/70‑1 (MQ=255) gCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAggtg > 1:1462404/1‑71 (MQ=255) gTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTgg > 1:1731850/1‑39 (MQ=255) gTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTgg > 1:1076102/1‑39 (MQ=255) gTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAgg < 1:952714/65‑1 (MQ=255) gTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAgg < 1:1527793/65‑1 (MQ=255) gATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCgaa < 1:1490680/70‑1 (MQ=255) gATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCaaa < 1:466091/70‑1 (MQ=255) gATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAAt > 1:2337072/1‑71 (MQ=255) aTGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAAt < 1:1896038/70‑1 (MQ=255) aTGGACGGCATGAACATCGTCGGCGACCTGTTCGGCCAGGGCAAGATGTTCCTGCCGCAGGTGGTCAAGTc > 1:3052628/1‑71 (MQ=255) tGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCg < 1:1788755/71‑1 (MQ=255) tGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCg < 1:1894212/71‑1 (MQ=255) gACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCgg > 1:812178/1‑70 (MQ=255) gCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGgcgcgc < 1:826985/71‑1 (MQ=255) gCATGAACATCGTCGGCGACCTGTTCGGCCAGGGCAAGATGTTCCTGCCGCAGGTGGTCAAGTCGgcgcgc < 1:1740951/71‑1 (MQ=255) gCATGAACATCGTCGGCGACCTGTTCGGCCAGGGCAAGATGTTCCTGCCGCAGGTGGTCAAGTCGgcgcgc < 1:1034951/71‑1 (MQ=255) tgtgGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAgg < 1:2642609/47‑1 (MQ=255) tgtgGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGa < 1:1450233/71‑1 (MQ=255) tgtgGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGa < 1:1900722/71‑1 (MQ=255) tgtgGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTC‑‑‑CCACAGGTGGTCAAATCGGCGCGCGTCATGaaa < 1:1843971/70‑1 (MQ=255) tgGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGaaa < 1:1982279/71‑1 (MQ=255) ggCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAgg > 1:2538665/1‑40 (MQ=255) ggCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGc < 1:3303601/71‑1 (MQ=255) gCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAAt > 1:1767415/1‑48 (MQ=255) gCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAAt > 1:3273683/1‑48 (MQ=255) gACCTGTTTGGCGAAGTGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGa < 1:2693425/61‑1 (MQ=255) gACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGt < 1:2160797/71‑1 (MQ=255) gTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTg < 1:271049/67‑1 (MQ=255) tttGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCaaa > 1:1837275/1‑39 (MQ=255) tttGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCaaa > 1:547654/1‑39 (MQ=255) tttGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCaaa > 1:612072/1‑39 (MQ=255) ttGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCAAGAAACAGGCGGTGGCCTAc > 1:2891284/1‑71 (MQ=255) tGGCGAATGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGCCTAc < 1:775866/70‑1 (MQ=255) tGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGaaa > 1:2994018/1‑55 (MQ=255) tGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGaaa > 1:943210/1‑55 (MQ=255) | AGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGCCTAC > minE/2781079‑2781218 |
| 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: ‑ |