Question: Variant Effect Predictor parameters for custom gff and fasta
1
gravatar for R.Blues
2.4 years ago by
R.Blues80
Somewhere
R.Blues80 wrote:

Hello everyone,

I am learning to use Ensemble VEP (Variant Effect Predictor) in command line mode. I am interested in predicting the effect of variants in an not-standard genome, with my own annotation.

I am using the VEP Github pre-release version (it is a new version, https://github.com/Ensembl/ensembl-vep). In that page it says that for a custom annotation and fasta you can use:

perl vep.pl -i input.vcf -gff data.gff.gz -fasta genome.fa.gz

This is not working for me (it has an error when inferring the transcripts from the GFF). If I do it as with the earlier VEP versions for a custom GFF ("-custom myFeatures.gff.gz,myFeatures,gff,overlap,0"), it creates an output, but I have the impression that this output does not use the fasta file for anything (if I remove the " -fasta genome.fa.gz" option, I obtain the same results).

So I am a bit clueless (there is a lot of information, but I have not found enough information for this case). How should I input these parameters? What am I doing wrong??

Thank you.

snp vep gff fasta • 1.2k views
ADD COMMENTlink modified 12 months ago by nicolas.dussex30 • written 2.4 years ago by R.Blues80

Given the lack of reply here, you may want to create an issue on GitHub as suggested in their README.

ADD REPLYlink written 2.4 years ago by Denise - Open Targets5.0k

CooVar was developed for annotating non-standard genomes. The paper contains an in-depth comparison with VEP results. (Disclaimer: I'm co-author of this paper).

ADD REPLYlink modified 2.4 years ago • written 2.4 years ago by Christian2.8k

You can also try variant_effect_predictor.pl from Ensembl (the antecesor of VEP). You can build a "cache" by providing a GFF and the reference genome sequence. Then you can run variant_effect_predictor.pl specifying that cache. This worked for me.

ADD REPLYlink written 2.4 years ago by abascalfederico1.1k
0
gravatar for nicolas.dussex
12 months ago by
New Zealand
nicolas.dussex30 wrote:

Hi R.Blues,

Did you manage to fix your issue? I am trying to do the same but I am having issues creating the cache and I think that it is due to the fact that my annotation has gene, mRNA, CDSs but not exons.

Cheers, Nic

ADD COMMENTlink written 12 months ago by nicolas.dussex30
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: 1049 users visited in the last hour