Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 1,151,925 | A→C | 21.2% | 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 | 21.2% | 43.1 / 17.7 | 33 | P207P (CCA→CCC) | fabD | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
Reads supporting (aligned to +/- strand): ref base A (17/9); new base C (2/5); total (19/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT > W3110S.gb/1151861‑1151994 | cgggcgCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGcc < 1:130251/71‑1 (MQ=255) cgcAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACa < 1:3724683/71‑1 (MQ=255) cgcAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACa < 1:4671555/71‑1 (MQ=255) aCGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACca > 1:2519904/1‑55 (MQ=255) gcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGAAGCCGACAAACTGGCAGTAg < 1:3317719/71‑1 (MQ=255) gcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGAAGCCGACAAACTGGCAGTAg < 1:3311985/71‑1 (MQ=255) gCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACa > 1:4169102/1‑56 (MQ=255) aCCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCg < 1:250927/70‑1 (MQ=255) cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa < 1:2442649/71‑1 (MQ=255) cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaa > 1:2853910/1‑70 (MQ=255) tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATc < 1:180685/71‑1 (MQ=255) tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATc < 1:3049954/71‑1 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:1617685/1‑71 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:4212783/1‑71 (MQ=255) gTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCtt < 1:1744139/71‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt < 1:1434667/71‑1 (MQ=255) ccGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa > 1:4191669/1‑71 (MQ=255) ccGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa > 1:1911013/1‑71 (MQ=255) ccGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa > 1:412501/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg < 1:386327/71‑1 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGc > 1:2919793/1‑54 (MQ=255) cTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:2390446/1‑57 (MQ=255) tgtgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa > 1:2056495/1‑70 (MQ=255) gtgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa > 1:3939416/1‑69 (MQ=255) tgCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAg > 1:1799611/1‑71 (MQ=255) cTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt < 1:3442940/54‑1 (MQ=255) cTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt < 1:656823/54‑1 (MQ=255) cTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTcc > 1:4452597/1‑71 (MQ=255) cTGAAGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTcc > 1:1615507/1‑71 (MQ=255) atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGt > 1:4400683/1‑65 (MQ=255) atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgt > 1:432010/1‑71 (MQ=255) tgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:4122144/1‑46 (MQ=255) gaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttg > 1:4434840/1‑71 (MQ=255) aaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAAt > 1:741028/1‑40 (MQ=255) cagcagCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAt < 1:3753915/71‑1 (MQ=255) | CGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT > W3110S.gb/1151861‑1151994 |
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: ‑ |