Question: Combining two VCFs each with own set of mutations
0
gravatar for vctrm67
16 months ago by
vctrm6720
vctrm6720 wrote:

I have two VCF files, each containing somatic mutation calls on each chromosome (ie. the VCFs are not separated by chromosome number or mutation position - they each contain calls on every chromosome) . I want to combine these two files into one. What is the best tool to do this? I tried bedtools concat but I get this message:

Different number of samples in /path/to/vcf. Perhaps "bcftools merge" is what you are looking for?
[W::bcf_hdr_merge] Trying to combine "AD" tag definitions of different lengths

Should I be using the merge function?

vcf • 569 views
ADD COMMENTlink modified 16 months ago by jared.andrews078.6k • written 16 months ago by vctrm6720
0
gravatar for jared.andrews07
16 months ago by
jared.andrews078.6k
Memphis, TN
jared.andrews078.6k wrote:

bcftools merge is likely what you want. It will preserve the VCF format and genotype information for each sample appropriately.

ADD COMMENTlink written 16 months ago by jared.andrews078.6k

They should be from the same sample. Does this make a difference?

To provide some clarification, these are VCFs that contain calls made by two different callers on the same sample. I want to combine them into one VCF and annotate that.

ADD REPLYlink written 16 months ago by vctrm6720

That's fine, it can still merge them and will even help you determine which caller called which variant (or where they overlap).

ADD REPLYlink written 16 months ago by jared.andrews078.6k
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: 1493 users visited in the last hour
_