Question: bcftools query missense_variant SNPs from gnomad VCF files
2
gravatar for Shicheng Guo
23 months ago by
Shicheng Guo8.5k
Shicheng Guo8.5k wrote:

Hi All,

I downloaded vcf files from gnomad , I want to filter some SNPs which are 'missense_variant' in vep annotation.

vep=A|missense_variant|MODERATE|ABCB1|ENSG00000085563|Transcript|ENST00000265724|

However, I cannot find any effective way to filter and extract these SNPs with bcftools

bcftools view -f PASS -i 'INFO/AF[0] > 0.1 & vep{1} =="missense_variant"' gnomad.exomes.r2.1.sites.chr13.vcf.bgz

Any suggestions?

Thanks.

vep bcftools • 1.5k views
ADD COMMENTlink modified 23 months ago by Emily_Ensembl21k • written 23 months ago by Shicheng Guo8.5k
1

Can you share the header and first few lines of the file?

ADD REPLYlink written 23 months ago by Arup Ghosh2.7k
7
gravatar for finswimmer
23 months ago by
finswimmer14k
Germany
finswimmer14k wrote:

Hello Shicheng Guo ,

try this:

$ bcftools view -f PASS -i 'INFO/AF[0] > 0.1 & INFO/vep ~ "missense_variant"' gnomad.exomes.r2.1.sites.chr13.vcf.bgz

fin swimmer

ADD COMMENTlink written 23 months ago by finswimmer14k

Thanks finswimmer. it works.

ADD REPLYlink written 23 months ago by Shicheng Guo8.5k

Please consider accepting as an answer: "Thanks" == "upvote"; "it works" == "accept"

ADD REPLYlink written 23 months ago by zx87549.9k

Hi Finswimmer, how to avoid multi-allelic SNPs with bcftools view or query? Thanks.

ADD REPLYlink written 23 months ago by Shicheng Guo8.5k

It looks I find it: bcftools view -m2 -M2 -v snps gnomad.exomes.r2.1.sites.chr13.vcf.bgz

ADD REPLYlink written 23 months ago by Shicheng Guo8.5k
6
gravatar for Emily_Ensembl
23 months ago by
Emily_Ensembl21k
EMBL-EBI
Emily_Ensembl21k wrote:

You can just use the filter options that ship with the VEP.

ADD COMMENTlink written 23 months ago by Emily_Ensembl21k
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: 1311 users visited in the last hour