bcftools merge creates high missingness rate
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?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
PCA Genotyping rate 1000genome and own data
written by curious_butterfly •
-
Merge only bim files with plink
written by mgdrnl •
-
How do I use PLINK to filter out SNPs?
written by melissachua90 •
-
The result of plink --freq is filled with NA
written by yoshifumimiya •
-
PCA for combined set
written by raalsuwaidi •
-
plink : Batch effect issues after merge of two datasets
written by Nicolas Rosewick
-
bcftools error duplicate samples
written by abyousaf •
-
bcftools merge duplicate names
written by evelyn
-
How to force missing values for genotype?
written by aritra90
-
vcftools vcf-merge error
written by biogirl
Hello,
could you please provide an example of what you mean by this?
Is
bcftools mergethe right command for you? This is used to combine multiple vcf files from different sample(names). Why do you use the--force-samplesparameter?fin swimmer