Question: Filtering large VCFs by allele freqency using bcftools split-vep function
gravatar for omid.alavijeh
8 months ago by
omid.alavijeh40 wrote:

Dear Biostars community,

I have a large VCF file that has been annotated using VEP. The intrinsic filtering function built into VEP whilst effective is very slow across large VCFs from my experience.

I have been playing around with the spli-vep plug-in on bcftools and found it to be much faster, however, I am struggling to get it to work with allele frequency.

The VCF is a standard format v4.1 VCF. I want to pull the alleles out that are <0.01 by MAX_AF. I am running the following script:

bcftools +split-vep vcf_file.vcf -i 'MAX_AF<0.01' -f -A -s worst

when running this I get an error saying 'MAX_AF' not defined in the vcf header but when running

bcftools +split-vep test/split-vep.vcf -l | less

It appears in the list.

Has anyone else used this tool to filter on allele_frequency?

Many thanks

vep bcftools split-vep • 308 views
ADD COMMENTlink written 8 months ago by omid.alavijeh40
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1495 users visited in the last hour