Biostar Beta. Not for public use.
Question: Exclude sites on the basis of the proportion of missing data using BCFTOOLS
Entering edit mode

I want to do something very simple. I want to keep sites that do not have any missing genotypes (i.e. 100% present across all samples).

In vcftools I usually do --max-missing 1.0 as indicated in the manual, but I want to use BCFTOOLS as it is much quicker. However, I cannot find the equivalent flag to do this. I have only found F_MISSING, which I think removes individuals and not sites.

Entering edit mode
$ bcftools view -e 'GT[*] = "mis"' multi-sample.vcf

This would exclude sites where any genotype is missing.

ADD COMMENTlink 8 months ago finswimmer 11k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0