bcftools merge creates high missingness rate
0
0
Entering edit mode
20 months ago
biogirl ▴ 190

Hi (again),

Posting this as a new topic. I've been trying to run an association test in plink, but it seems the rate of missingness in the merge vcf is high, meaning all copies of the second allele have been replaced with missing genotypes. Therefore, all of my variants have a MAF of 0.

This is the command I've used to create a merged vcf:

bcftools merge --force-samples *vcf.gz -o all_isolates.vcf

Which I then convert to a bed using --make-bed in plink.

Any ideas as to why I'm getting such a high level of missingness in the merged vcf though?

bcftools plink vcftools • 398 views
ADD COMMENT
0
Entering edit mode

Hello,

meaning all copies of the second allele have been replaced with missing genotypes

could you please provide an example of what you mean by this?

Is bcftools merge the right command for you? This is used to combine multiple vcf files from different sample(names). Why do you use the --force-samples parameter?

fin swimmer

ADD REPLY

Login before adding your answer.

Traffic: 2551 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