Biostar Beta. Not for public use.
Script: VCF to 3 base mutation pattern
0
Entering edit mode
13 months ago
Shicheng Guo ♦ 7.5k

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 • 157 views
ADD COMMENTlink
1
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

What about indel data? Discard?

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
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 REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1