Question: Merging VCFs: Illumina/GATK and LoFreq
0
gravatar for ari.nazarian
23 months ago by
ari.nazarian10
ari.nazarian10 wrote:

I have VCFs from a number of different tumor samples sourced mostly from two different variant callers that I'm trying to merge onto one sheet: GATK (via an Illumina MiSeq machine) and LoFreq. The end goal is to upload the merged VCF to CRAVAT for easier visualization/interpretation of variant information.

I'm struggling to figure out how to merge variants, though---neither vcf-merge nor CombineVariants have been able to do this successfully.

Does anyone have any experience with merging these VCFs generated by these specific variant callers?

miseq vcf variant calling • 942 views
ADD COMMENTlink written 23 months ago by ari.nazarian10

I'm struggling to figure out how to merge variants, though---neither vcf-merge nor CombineVariants have been able to do this successfully.

so, what is the problem here.

ADD REPLYlink written 23 months ago by Pierre Lindenbaum118k

I've just run into so many different errors that I was hoping I didn't have to go into the weeds but here's a brief overview of the initial issue I'm having (not even considering the errors when I try to manipulate the vcfs, like adding columns with sed and so forth:

I ran vcf-merge on two vcfs, one generated from GATK: GATK VCF ...another from LoFreq: LOFREQ VCF ...and this was the output merge.vcf: BADMERGE VCF ...and the corresponding output in Terminal: MERGE OUTPUT I want them to look more like this: GOOD MERGE

I've tried adding an extra column and Sample ID field to the LoFreq one but to no avail. I'm looking for pointers on how I can add columns and a sample ID to LoFreq vcfs to successfully merge this with GATK vcfs.

ADD REPLYlink written 23 months ago by ari.nazarian10

It's clearly not the same 'dictionary'; one set of VCF have a 'chr' prefix in the chromosome names, the other not. They should all have the very same names than defined in the REFerence file

ADD REPLYlink modified 23 months ago • written 23 months ago by Pierre Lindenbaum118k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2458 users visited in the last hour