Hello all, I have a set of VCF files of normal and tumor samples. Are there any tools available to get the somatic only variants from the vcf files? I have checked the post here but found it inconclusive.
•
written
by
Sandeep
1 answer
Sounds like you have unpaired data and just want to filter for variants that are only called in tumour?
I would merge them all together (bcftools merge), or just merge them into their respective tumour and normal groupings; so, you would have 2 final VCFs: tumour.vcf and normal.vcf. I would then create a unique key that identifiers each variant found in each group (CHR:POS:REF:ALT), and then cross-compare with a simple script like awk. You could also set the ID field to be this unique key.
Another way is to merge all files together and then use SnpSift CaseControl, but only checking the allele tallies at the end.
Kevin
Log in to answer this question.
-
Allele Frequency >0.5 for cancer somatic mutations?
written by Rozita •
-
VarScan.v2.4.3 generating incomplete or partial VCF file
written by Sandeep
-
The different results for the same tumor sample between VarScan2 somatic and mpileup2snp
written by rayliu0115 •
-
identifying TP53 germline mutations which are also in tumor cells.
written by lait
-
converting wannovar generated csv annotation data to a vcf=4.1 file.
written by sktbanerjee1 •
-
how to interprete the somatic mutations with low frequency in tumor sample compared with the normal…
written by Ontheway •
-
Understanding SomaticSniper output
written by ttom
-
Calling variants in blood tumor samples without matched normal samples
written by umn_bist •
-
Calling germline variants from Tumor/Clone pairs to find the somatic mutations of only tumor and on…
written by ivivek_ngs
-
Breakdancer Somatic Sv Input Order Of Tumor And Normal Bam Files
written by g.skullcap •
