I'm interested in taking .vcf files obtained from TCGA and merging them into a single file for a group of the same samples. TCGA used 4 variant callers (Somatic Sniper, MuSE, MuTect2, and VarScan2). Is there a script that can do this?
•
written
by
j.lunger18 •
1 answer
If you have access to the VCFs, then I assume that you have access to the BAMs? My preference, time-permitting, would be to re-call somatic variants from the BAMs using the same variant caller, but I'm aware that this may pose a logistical challenge (in addition to taking more time).
In any case, you can do this yourself by setting rules about which variants to keep (or not). Questions you need to ask (and answer):
- are all samples profiled by each method? If not, which samples are profiled across more than one method?
- if a sample has results across just 2 methods, should I only retain variants identified by both methods, or just take a Union of the variants identified by both?
BCFtools will be your friend for getting through this, and awk.
Kevin
Log in to answer this question.
-
MuSE calls too few snps
written by yalin •
-
Somatic-Sniper calls very few variants
written by runfreely2 •
-
How to handle VCFs from the same sample but using different aligners and variant callers?
written by aldhairmedico •
-
Understanding TCGA somatic mutation files
written by j.lunger18 •
-
Merge vcf files of structural variants
written by ma17 •
-
Tools to filter Mutect2 output
written by nanana
-
Qustions about Somatic Variant Caller
written by Phoe •
-
annotating TCGA VCF file
written by jan
-
about GATK VariantAnnotator and VCF files
written by Bogdan
-
TCGA germline variants status
written by jan
