This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Vcf-compare two multi-sample vcf files

Hi,

I have two large multi-sample vcf files. The samples in both files are the same except the snp's were called using different programs. i want to see if they match. I am not sure if vcf-compare works. I gave it a shot and it returned all zeros. Is there any other way i can compare two multi-sample vcf files?

Thanks

vcf-compare vcftools
• 5,695 views • link
updated by maduh17 • maduh17 • written by AB AB

2 answers

Is there any other way i can compare two multi-sample vcf files?

you could use picard/GenotypeConcordance

Evaluate genotype concordance between callsets.This tool evaluates the concordance between genotype calls for samples in different callsets where one is being considered as the truth (aka standard, or reference) and the other as the call that is being evaluated for accuracy.

java -jar picard.jar GenotypeConcordance \
      CALL_VCF=input.vcf \
      CALL_SAMPLE=sample_name \
      O=gc_concordance.vcf \
      TRUTH_VCF=truth_set.vcf \
      TRUTH_SAMPLE=truth_sample#
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

!!!!!!!!!!!!!!!!! I found a link saying that GATK's GenotypeConcordance is being outphased: https://gatkforums.broadinstitute.org/gatk/discussion/5795/genotype-concordance-output

So might wanna just use picard as Pierre suggested

• 0 views • link
written by maduh17 • maduh17

Log in to answer this question.

More posts like this