| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,235,971 | C→T | 13.0% | R373C (CGT→TGT) | purH → | fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,235,971 | 0 | C | T | 13.0% | 152.1 / 12.7 | 54 | R373C (CGT→TGT) | purH | fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (29/18); new base T (4/3); total (33/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01 | |||||||||||
AGAAGCCCTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGT > minE/2235904‑2236041 | agaagCCCTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGtgtg > 1:2537105/1‑71 (MQ=255) gaagCCCTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCgtgt < 1:1498413/71‑1 (MQ=255) aagCCCTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGtt < 1:1270626/71‑1 (MQ=255) gCCCTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGTGTGTTcc < 1:2766950/71‑1 (MQ=255) ccTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCgg > 1:3034787/1‑71 (MQ=255) cTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCggg < 1:3398766/71‑1 (MQ=255) cTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCggg > 1:1139204/1‑71 (MQ=255) cTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCggg < 1:2709292/71‑1 (MQ=255) cTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCggg < 1:2820495/71‑1 (MQ=255) gAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGTGTGTTCCGGGcc > 1:1887310/1‑71 (MQ=255) gAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGcc > 1:2853103/1‑71 (MQ=255) gAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGcc > 1:927885/1‑71 (MQ=255) gAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGc > 1:3306199/1‑70 (MQ=255) aaaaTCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGcc < 1:2997584/70‑1 (MQ=255) aaaaTCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCt > 1:1337858/1‑71 (MQ=255) aaaaTCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCt > 1:2206090/1‑71 (MQ=255) aaTCACCGCCGCCAAACAGAACGTACGCGTTGTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTc > 1:285113/1‑70 (MQ=255) aaTCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCgt > 1:682060/1‑58 (MQ=255) aaTCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTc > 1:336210/1‑70 (MQ=255) aCCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGAttt > 1:3092058/1‑71 (MQ=255) cgccAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAAc > 1:4115013/1‑71 (MQ=255) ccAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCaa > 1:1060366/1‑67 (MQ=255) cAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGa < 1:3367503/60‑1 (MQ=255) aaaCAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGTGTGTTCCGGGCCTCGATTTCAAACGCGt > 1:2886104/1‑71 (MQ=255) aaCAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGTGTGTTCCGGGCCTCGATTTCAAACGCGt > 1:1800101/1‑70 (MQ=255) aCAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTg > 1:1162029/1‑70 (MQ=255) gAACGTACGCGTTCTGACCTGCGGTCTGTGGGGCGAGTGTGTTCCGGGCCTCGATTTCAAACGCGTGAAcg > 1:2159452/1‑71 (MQ=255) aaCGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGTGTGTTCCGGGCCTCGATTTCAAACGCGTGAAcgg < 1:2858028/71‑1 (MQ=255) aaCGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAAcgg > 1:3316258/1‑71 (MQ=255) aaCGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAAcg > 1:1731217/1‑70 (MQ=255) aCGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAAcgg > 1:1989165/1‑70 (MQ=255) cGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAAcggcg > 1:173329/1‑71 (MQ=255) aCGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCa > 1:2174206/1‑53 (MQ=255) gcgTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTctg > 1:3481278/1‑71 (MQ=255) cTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGt > 1:1625716/1‑63 (MQ=255) cTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTct > 1:3803670/1‑65 (MQ=255) tGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGAttt < 1:323046/43‑1 (MQ=255) gACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGtt > 1:2363631/1‑70 (MQ=255) gACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGtt > 1:2315163/1‑70 (MQ=255) aCCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTc < 1:3751408/36‑1 (MQ=255) aCCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTc < 1:960575/70‑1 (MQ=255) cTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCa < 1:2079191/41‑1 (MQ=255) cTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCa < 1:1344913/41‑1 (MQ=255) gTCAGTGGGGCGAGTGTGTCCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGt < 1:2231104/71‑1 (MQ=255) cAGTGGGGCGAGCGTGTTCCTGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGt > 1:1044277/1‑69 (MQ=255) cAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGa > 1:1768448/1‑71 (MQ=255) cAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGa > 1:358968/1‑71 (MQ=255) cAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGa > 1:1723121/1‑71 (MQ=255) cAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGa > 1:1094720/1‑71 (MQ=255) gTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGAcc < 1:3504624/71‑1 (MQ=255) gTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGAcc < 1:3909231/71‑1 (MQ=255) gggCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTgg < 1:3996814/71‑1 (MQ=255) cGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGa < 1:1005295/71‑1 (MQ=255) gCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGa < 1:144121/54‑1 (MQ=255) gCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGa < 1:3985982/54‑1 (MQ=255) cGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGt < 1:3520467/71‑1 (MQ=255) | AGAAGCCCTGAAAATCACCGCCGCCAAACAGAACGTACGCGTTCTGACCTGCGGTCAGTGGGGCGAGCGTGTTCCGGGCCTCGATTTCAAACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGT > minE/2235904‑2236041 |
| 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: ‑ |