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
0 answers
No answers yet.
Log in to answer this question.
-
Filtering VCF to divide with equal sizes
written by avelarbio46 •
-
Counting the specific number of a type of variant per sample in a multi-sample VCF
written by tacrolimus
-
bcftools split-vep -- how to split the INFO column up and also assign the ensembl vep headers
written by hello_bioinf •
-
Count number of occurrences an annotation in a multisample vcf split into multiple chunks
written by tacrolimus
-
Extract samples with alternative genotype from a multi-sample VCF
written by tacrolimus
-
Finding unique variants by sample ID in a multisample vcf
written by tacrolimus
-
Counting allele counts/frequencies from a vcf file with both cases and controls
written by tacrolimus
-
Faster filtering of vcfs using split-vep over filter_vep
written by tacrolimus
-
Filtering a multi-person structural variant VCF file by type of SV
written by tacrolimus
-
Bcftools merge taking too much time and producing large file
written by BAGeno
Do you have a solution?
Afraid not no sorry!