Script: VCF to 3 base mutation pattern

0

Entering edit mode

5.5 years ago

Shicheng Guo ★ 9.6k

Hi All,

Who can share a script to transfer mutation VCF files to 3-base mutation pattern:

Bekow is the output i want to obtain. I want the counts for each mutation category.

enter image description here

Thanks.

vcf • 1.2k views

ADD COMMENT • link 5.5 years ago by Shicheng Guo ★ 9.6k

1

Entering edit mode

Provide example input, expected output, and what have you tried so far.

ADD REPLY • link 5.5 years ago by zx8754 12k

0

Entering edit mode

What about indel data? Discard?

ADD REPLY • link 5.5 years ago by MatthewP ★ 1.4k

0

Entering edit mode

Yes，discard indel and complex mutation type. another problem maybe for one mutation, it may have 2 annation.

ADD REPLY • link 5.5 years ago by Shicheng Guo ★ 9.6k

0

Entering edit mode

What are the numbers that are attached to each change? These 3-base changes - are they codon changes or just +/- 1 nucleotide? Are you looking for changes in the reference strand only or is this strand specific?

More importantly, why are you doing this exercise? What do you hope to accomplish from it?

ADD REPLY • link 5.5 years ago by Ram 44k

Login before adding your answer.