Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

how to run vcf2maf with clinvar db #327

Open
songeric1107 opened this issue Aug 25, 2022 · 2 comments
Open

how to run vcf2maf with clinvar db #327

songeric1107 opened this issue Aug 25, 2022 · 2 comments
Assignees

Comments

@songeric1107
Copy link

I have vcf files which are generated using freebay, how could I add the clinvar annotation to the vcf and convert to maf?

I try to run as following

vep --input_file JH020.soma.freebayes.vcf --output_file JH020.soma.clinc.freebayes.vcf --format vcf --species homo_sapiens --custom /vep_ref/hg19/new/clinvar.vcf.gz,ClinVar,vcf,exact,0,CLNSIG,CLNREVSTAT,CLNSIGCONF

all the variant are intergenic_variant, which is not right

@ckandoth ckandoth self-assigned this Aug 25, 2022
@ckandoth
Copy link
Collaborator

You are using VEP, not vcf2maf. But I'll try to help. If all the variants are showing as intergenic, then it is likely that you are using a transcript/reference that is a different human assembly than your VCF. E.g. hg19 VCF with hg38 VEP cache.

@songeric1107
Copy link
Author

yes, I have tried to add the option to vcf2maf with the customed database which points to clinvar, but always failed. What is the best way to add the clinvar in vcf2maf?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants