Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 34,427 | C→T | 25.0% | R364H (CGC→CAC) | caiC ← | predicted crotonobetaine CoA ligase:carnitine CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 34,427 | 0 | C | T | 25.0% | 60.2 / 18.4 | 36 | R364H (CGC→CAC) | caiC | predicted crotonobetaine CoA ligase:carnitine CoA ligase |
Reads supporting (aligned to +/- strand): ref base C (14/13); new base T (2/7); total (16/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.45e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTT > minE/34362‑34496 | tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGAt < 1:3302761/69‑1 (MQ=255) tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGa < 1:2037027/68‑1 (MQ=255) aaGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATtgtg < 1:2150807/71‑1 (MQ=255) aaGATGGTGTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATtgtg < 1:2769180/71‑1 (MQ=255) ttttCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGa < 1:1298080/65‑1 (MQ=255) tCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTcgc > 1:3873042/1‑69 (MQ=255) ccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGa < 1:1055756/70‑1 (MQ=255) gTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGa < 1:1147074/50‑1 (MQ=255) tATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCg < 1:3799434/71‑1 (MQ=255) tGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCt < 1:886354/71‑1 (MQ=255) ccTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTc < 1:3586117/71‑1 (MQ=255) cTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCg > 1:2089950/1‑71 (MQ=255) ttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTcc > 1:2658201/1‑63 (MQ=255) ttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCg > 1:3727501/1‑69 (MQ=255) tGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa < 1:158468/71‑1 (MQ=255) gCAGATTTCACCGATCTCCCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa > 1:2519795/1‑70 (MQ=255) cAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaaaa < 1:1856596/71‑1 (MQ=255) cAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaaa < 1:1205126/70‑1 (MQ=255) aTTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAc > 1:2196250/1‑69 (MQ=255) ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCg > 1:2720048/1‑53 (MQ=255) ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa < 1:2254911/64‑1 (MQ=255) cACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCAcc > 1:2761246/1‑71 (MQ=255) cGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGa > 1:3130442/1‑71 (MQ=255) cGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGa > 1:3396342/1‑71 (MQ=255) cGATCTCACCAGCCGGGAGCGCGCGATTGTGATCGTCGCGGATCTCCGCt > 1:3198314/1‑50 (MQ=255) ctcACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGAtct < 1:3164973/41‑1 (MQ=255) ctcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCAGCTTCGTAGCAAAACCCCAccc < 1:750135/65‑1 (MQ=255) cACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATc > 1:212119/1‑71 (MQ=255) aCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCaa > 1:416049/1‑68 (MQ=255) aCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCa < 1:2889345/67‑1 (MQ=255) ccAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGa < 1:3683040/71‑1 (MQ=255) cAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGAc < 1:3800268/71‑1 (MQ=255) cAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGAc > 1:3153023/1‑71 (MQ=255) cAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGAc < 1:2025730/71‑1 (MQ=255) cGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGc > 1:2355506/1‑70 (MQ=255) gCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGtt > 1:937090/1‑71 (MQ=255) | TGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTT > minE/34362‑34496 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |