bcftools merged vcf file assigns all variants to one sample
Entering edit mode
4 months ago

I've made one vcf file for each of three samples. I then combined them using bcftools, like so:

# Make a list of vcf files to merge
cat "${OUT}/results/variants/vcf_list"


Then merge the list:

bcftools merge -l "${OUT}/results/variants/vcf_list" -Ov -o "${OUT}/results/variants/merge_individuals.vcf"

And index it:

tabix -p vcf "${OUT}/results/variants/merge_individuals.vcf.gz"

The resulting merged vcf appears to have three columns, one for each sample. When I open it in IGV all the variants are assigned to sample 3a7a (see image)

enter image description here

But when I open it in GenomeBrowse I see them assigned correctly to each of the three samples:

enter image description here

I can't work out what's going on? I've uploaded the individual and merged vcf files here. Any help much appreciated


vcf merge gatk igv bcftools • 132 views

Login before adding your answer.

Traffic: 2034 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6