Question: Merge multiple vcfs vgtools
0
gravatar for arunkumarramrsh
3 months ago by
arunkumarramrsh40 wrote:

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 • 209 views
ADD COMMENTlink modified 3 months ago by glenn.hickey160 • written 3 months ago by arunkumarramrsh40
1
gravatar for glenn.hickey
3 months ago by
glenn.hickey160
glenn.hickey160 wrote:

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 COMMENTlink modified 3 months ago by RamRS30k • written 3 months ago by glenn.hickey160
0
gravatar for glenn.hickey
3 months ago by
glenn.hickey160
glenn.hickey160 wrote:

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 COMMENTlink written 3 months ago by glenn.hickey160
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: 1570 users visited in the last hour