How do I call somatic indels from a VCF file?
Entering edit mode
6.2 years ago


I'm doing an analysis to find Indels in a cancer sample. I've used Platypus to call variants and VCF tools to filter out all variations except indels. I have matched normal exome data, and the recommendation is to call the indels on both the tumor and normal sample at the same time against the reference genome (hg38 in my case), then use the script to call the actual mutations based on the VCF output from Platypus. In trying to run that extension script, I get the error:

Could not find normal sample normaldna in input VCF header

I get the same error with the tumor sample. I've tried the original file name (normaldna.bam), the file name without the extension (normaldna), and the full file path (/dir/that/contains/normaldna.bam). I already looked through the VCF input file, and I can't find where it is supposed to specify the original BAM a particular indel was called for. Does anyone know of any other post-processing tool of a VCF I could do to accomplish this, or have any experience with this script in particular? I know there are dozens of indel callers, but I'm pretty set on Platypus.

VCF cancer mutation indel platypus • 1.9k views
Entering edit mode

Update: You can see Andy Rimmer's answer to this question on the Google Group for Platypus here:!topic/platypus-users/1kIfOBDfkfA


Login before adding your answer.

Traffic: 1615 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6