Question: bcftools merge creates high missingness rate
0
gravatar for biogirl
14 months ago by
biogirl190
European Union
biogirl190 wrote:

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?

vcftools plink bcftools • 277 views
ADD COMMENTlink written 14 months ago by biogirl190

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 REPLYlink written 14 months ago by finswimmer14k
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: 1669 users visited in the last hour