| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 734,251 | C→T | pseudogene (32/1434 nt) | ybfO → | conserved protein, rhs‑like YbfO |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 734,251 | 0 | C | T | 92.9% | 22.0 / ‑3.6 | 14 | pseudogene (32/1434 nt) | ybfO | conserved protein, rhs‑like YbfO |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (8/5); total (9/5) | |||||||||||
| 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.39e-01 | |||||||||||
CTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCA > NC_000913/734214‑734328 | cTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACg > 2:151339/1‑78 (MQ=14) cAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGgaga > 2:151042/1‑80 (MQ=34) aGCATGTGGCCGGATAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGgaga < 1:151305/79‑1 (MQ=21) ggATAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt > 1:151331/1‑81 (MQ=18) ggATAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt > 1:151350/1‑81 (MQ=18) ggATAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt > 2:151327/1‑81 (MQ=18) aTAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc < 2:151306/81‑1 (MQ=18) tAACCGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc < 1:151307/81‑1 (MQ=18) cGTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa < 1:151308/81‑1 (MQ=18) gTATCGCCCGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAg > 2:151338/1‑81 (MQ=18) cGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATc < 2:151309/81‑1 (MQ=18) gTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc > 2:151340/1‑81 (MQ=18) gACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCg > 2:151351/1‑80 (MQ=21) aCGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCa > 2:151337/1‑81 (MQ=18) | CTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCA > NC_000913/734214‑734328 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |