| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,151,925 | A→C | 30.0% | P207P (CCA→CCC) | fabD → | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,151,925 | 0 | A | C | 30.0% | 8.8 / 17.7 | 20 | P207P (CCA→CCC) | fabD | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
| Reads supporting (aligned to +/- strand): ref base A (5/9); new base C (2/4); total (7/13) | |||||||||||
| 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 | |||||||||||
GCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAA > W3110S.gb/1151866‑1151993 | gcAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACaa < 1:1042162/71‑1 (MQ=255) cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa < 1:1557573/71‑1 (MQ=255) cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:722734/70‑1 (MQ=255) tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATc < 1:2190718/71‑1 (MQ=255) gTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCtt < 1:2494604/71‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt < 1:1776667/71‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa > 1:2476635/1‑63 (MQ=255) ccGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa < 1:480052/71‑1 (MQ=255) tcACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:1656868/1‑60 (MQ=255) aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAAt > 1:1007090/1‑45 (MQ=255) aCTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGaa < 1:2102396/44‑1 (MQ=255) aCTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCAcc < 1:187784/70‑1 (MQ=255) aCTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCAcc > 1:2297680/1‑70 (MQ=255) aCTGTGAGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCAcc > 1:816/1‑70 (MQ=255) cTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa < 1:1998566/71‑1 (MQ=255) tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAAc < 1:1560225/69‑1 (MQ=255) cgcTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt < 1:1136249/56‑1 (MQ=255) tgatgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:661879/43‑1 (MQ=255) atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa > 1:1758181/1‑43 (MQ=255) aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt > 1:2440652/1‑71 (MQ=255) cagcagCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGaa > 1:1111165/1‑70 (MQ=255) cagcagCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGaa > 1:1049487/1‑70 (MQ=255) agcagcCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgtt > 1:1875320/1‑64 (MQ=255) | GCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAA > W3110S.gb/1151866‑1151993 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |