Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,780,693 | G→A | 26.8% | A556A (GCG→GCA) | 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,780,693 | 0 | G | A | 26.8% | 88.3 / 28.9 | 78 | A556A (GCG→GCA) | metH | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent |
Reads supporting (aligned to +/- strand): ref base G (46/11); new base A (17/4); total (63/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATC > minE/2780630‑2780745 | caacaaCTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAt > 1:578325/1‑69 (MQ=255) caacaaCTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAt > 1:293516/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGAttcc > 1:1532143/1‑67 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:2275530/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:1630963/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:1578533/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:3146101/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:1987524/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:526380/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:570884/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:2004006/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:1090750/1‑69 (MQ=255) caatTACGCACAGGACTTTATCGGCGCGTATGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTc > 1:1145402/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:1955838/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:2106008/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:2333671/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:1008746/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:1879428/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:2988190/1‑69 (MQ=255) acTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcg > 1:1482367/1‑69 (MQ=255) aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTgcta > 1:2318153/1‑59 (MQ=255) aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt > 1:167108/1‑62 (MQ=255) aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt > 1:2217370/1‑62 (MQ=255) aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt > 1:2625434/1‑62 (MQ=255) aCGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt > 1:3424644/1‑62 (MQ=255) cgcATGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg < 1:1706355/67‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg < 1:2877195/67‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg < 1:1147613/67‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg < 1:3168474/67‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCcggcgg < 1:646857/67‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg < 1:1993706/69‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg < 1:1215262/69‑1 (MQ=255) cgcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg < 1:306019/69‑1 (MQ=255) gcAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg < 1:2142107/68‑1 (MQ=255) aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:2864585/1‑69 (MQ=255) aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:2511895/1‑69 (MQ=255) aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:2342857/1‑69 (MQ=255) aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:1044215/1‑69 (MQ=255) aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:1076592/1‑69 (MQ=255) aCTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:1318247/1‑69 (MQ=255) cTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTa > 1:633209/1‑68 (MQ=255) gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCcggcgg > 1:1524191/1‑50 (MQ=255) gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc > 1:3140688/1‑69 (MQ=255) gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc > 1:3175312/1‑69 (MQ=255) gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc > 1:2073586/1‑69 (MQ=255) gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc > 1:1694831/1‑69 (MQ=255) gcgTGTGAAGACATCAAACGTGAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGtgtcgttc > 1:370677/1‑69 (MQ=255) gtgtGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGgcggcg < 1:2772742/50‑1 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGGGTATCTAACGTTTCTTTCTCGTTcc > 1:224672/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:3426869/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:1930512/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:2925311/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:19044/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:2654285/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:259632/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:1262110/1‑69 (MQ=255) aaGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTcc > 1:1649691/1‑69 (MQ=255) tCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCa < 1:554944/69‑1 (MQ=255) aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCt > 1:2839176/1‑42 (MQ=255) aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc > 1:2115145/1‑69 (MQ=255) aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc > 1:1553536/1‑69 (MQ=255) aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc > 1:2388940/1‑69 (MQ=255) aaaCGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAAc > 1:3007848/1‑69 (MQ=255) gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc < 1:646523/67‑1 (MQ=255) gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc < 1:1229258/67‑1 (MQ=255) gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc < 1:517945/67‑1 (MQ=255) gtgAACTGCCGCACGCACTGATTTCCGGCGGCGTTTCTAACGTGTCGTTCTCGTTTCGTGGCAACGATc < 1:2103134/67‑1 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:699103/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:2383541/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:1123011/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:2461658/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:376749/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:268059/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:198392/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:3299863/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:3255811/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:3253766/1‑59 (MQ=255) gCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACg > 1:3028480/1‑59 (MQ=255) | CAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATC > minE/2780630‑2780745 |
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: ‑ |