Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,014,611 | A→G | 22.9% | G486G (GGT→GGC) | fumA ← | aerobic Class I fumarate hydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,014,611 | 0 | A | G | 22.9% | 83.8 / 23.8 | 48 | G486G (GGT→GGC) | fumA | aerobic Class I fumarate hydratase |
Reads supporting (aligned to +/- strand): ref base A (20/17); new base G (8/3); total (28/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01 |
TCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGC > minE/1014542‑1014678 | tCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccac > 1:3749888/1‑71 (MQ=255) ccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTTCGCCAATACAGCGGCCGGaccacc > 1:999837/1‑70 (MQ=255) cAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa > 1:998939/1‑71 (MQ=255) aGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa > 1:2439872/1‑70 (MQ=255) ccGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa > 1:3498711/1‑66 (MQ=255) gATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCa < 1:1178934/71‑1 (MQ=255) atatTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc > 1:321020/1‑69 (MQ=255) tCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAgg > 1:1381880/1‑69 (MQ=255) cAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTa < 1:2045148/70‑1 (MQ=255) aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg < 1:3894155/70‑1 (MQ=255) aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg < 1:192785/70‑1 (MQ=255) caTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgcc < 1:2294078/71‑1 (MQ=255) aTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccg > 1:3128097/1‑71 (MQ=255) ttCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccg > 1:1118823/1‑70 (MQ=255) ccAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGccgccg > 1:122986/1‑71 (MQ=255) cAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt < 1:731203/71‑1 (MQ=255) aGGCTCTTAATACTTCCCTGCGCCATTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtg < 1:983636/71‑1 (MQ=255) ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt > 1:780063/1‑71 (MQ=255) ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt > 1:705820/1‑71 (MQ=255) ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt > 1:481533/1‑71 (MQ=255) ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt > 1:2730254/1‑71 (MQ=255) gCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt < 1:3249765/70‑1 (MQ=255) tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGa > 1:304612/1‑57 (MQ=255) tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt < 1:1814657/66‑1 (MQ=255) tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtttt > 1:2830426/1‑71 (MQ=255) tAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa > 1:2726405/1‑71 (MQ=255) aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtt > 1:927201/1‑65 (MQ=255) aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa > 1:522316/1‑70 (MQ=255) aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCAAAGGTAGAAGCCGCCGTGtttttt < 1:1505069/69‑1 (MQ=255) aTACTTCCCTACGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc < 1:1032632/70‑1 (MQ=255) tgcTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc > 1:689321/3‑71 (MQ=255) cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa < 1:1169085/69‑1 (MQ=255) ccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTa > 1:3333425/1‑43 (MQ=255) cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACacg > 1:2586187/1‑63 (MQ=255) cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcg < 1:3625366/67‑2 (MQ=255) cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg > 1:3057193/1‑69 (MQ=255) cagtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct > 1:1661580/4‑71 (MQ=255) aaTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa > 1:1843042/1‑54 (MQ=255) gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctg > 1:283608/2‑70 (MQ=255) tACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgct > 1:3128585/1‑71 (MQ=255) cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg < 1:3721172/70‑1 (MQ=255) cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggctg < 1:2730593/71‑1 (MQ=255) cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggct > 1:3006410/1‑70 (MQ=255) gCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCg > 1:1321545/1‑71 (MQ=255) ccGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCgg < 1:421606/71‑1 (MQ=255) cGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGt < 1:1790629/71‑1 (MQ=255) ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt < 1:367644/71‑1 (MQ=255) ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGt > 1:2382111/1‑70 (MQ=255) gACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtt < 1:3204316/36‑1 (MQ=255) gACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc < 1:1575686/42‑1 (MQ=255) accgccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGttac > 1:2005342/1‑69 (MQ=255) accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc < 1:2117837/71‑1 (MQ=255) accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc > 1:2014914/1‑71 (MQ=255) accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc < 1:1255534/71‑1 (MQ=255) ccaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt > 1:3467307/1‑50 (MQ=255) | TCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGC > minE/1014542‑1014678 |
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: ‑ |