Question: Consensus VCF file from several variant callers
gravatar for ajuiwl
16 months ago by
ajuiwl30 wrote:

Hi all, I want to do a consensus VCF file with only the variants that have been called by two variant callers. The problem is that for the same variant the INFO and FORMAT columns of the VCF files show different values, this is to say: GATK shows the variant with DP, MQ, FS values that differ from those shown by Freebayes. I understand that the programs have different internal parameters. But in order to obtain a single VCF file with just one value per each property (DP, MQ, FS etc...) what should I do?

I have checked the GIAB consensus vcf file and they just show one value per each property despite the fact that they have use several variant callers, but I have not found the method that they followed in order to merge that.

I would be very thankful if somebody could help in this issue.

ADD COMMENTlink modified 16 months ago by Pierre Lindenbaum127k • written 16 months ago by ajuiwl30
gravatar for Pierre Lindenbaum
16 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum127k wrote:

GATK SelectVariants


 --concordance Freebayes.vcf  

A site is considered concordant if (1) we are not looking for specific samples and there is a variant called in both the variant and concordance tracks or (2) every sample present in the variant track is present in the concordance track and they have the sample genotype call.

ADD COMMENTlink written 16 months ago by Pierre Lindenbaum127k
Please log in to add an answer.


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