Merge multiple vcfs vgtools
2
0
Entering edit mode
9 months ago

Hi, Is it possible to combine vcfs produced by vg tools, all mapped using the same variant graph? Simply using bcftools etc does not seem to work since the reference calls differ.

vg • 336 views
ADD COMMENT
1
Entering edit mode
8 months ago
glenn.hickey ▴ 240

The -a option can now be used with vg call to add reference calls: https://github.com/vgteam/vg/pull/2966

vg call -a ... -s sample1 > sample1.vcf
vg call -a ... -s sample2 > sample2.vcf

etc.

Then bcftools merge sample1.vcf.gz sample2.vcf.gz -m all to make the combined VCF.

ADD COMMENT
0
Entering edit mode
9 months ago
glenn.hickey ▴ 240

This only works right now in "genotyping mode", where you construct a graph from a VCF with vg construct -a then genotype that same VCF with vg call -v.

I'm hoping to have a solution for general graphs, in the form of an option that makes vg call call every site with consistent positions no matter what, within the next week or so. It's coming up as a problem for many people.

ADD COMMENT

Login before adding your answer.

Traffic: 2799 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6