Question: Working with SNPEFF for Annotating Variant Calling
gravatar for ajayi
4.6 years ago by
South Africa
ajayi0 wrote:

Hi all,

Is it possible to run snpEFF on VCF file directly without the database? What is the exactly command line for variant annotation using snpEFF? I am kind of confuse with the command-line I saw on snpEFF webpage

Here is the command I used:

java -jar -Xmx29g .../snpEff/snpeff.jar .../reference/HumanTest500k_g1k_H37Rv_decoy.fasta ../fastqs/H37Rv1116.raw.snps.indels.g.vcf > ../H37Rv1116.vcf

Here is my error:

java.lang.RuntimeException: Property: '/usr/people/ajayi/test/enter code hereexample/fastqs/H37Rv_Raw_reads/H37Rv1116.raw.snps.indels.g.vcf.genome' not found at ca.mcgill.mcb.pcingola.interval.Genome.<init>( at ca.mcgill.mcb.pcingola.snpEffect.Config.readGenomeConfig( at ca.mcgill.mcb.pcingola.snpEffect.Config.readConfig( at ca.mcgill.mcb.pcingola.snpEffect.Config.init( at ca.mcgill.mcb.pcingola.snpEffect.Config.<init>( at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.loadConfig( at at at at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(

ADD COMMENTlink modified 4.6 years ago • written 4.6 years ago by ajayi0

In case you don't solve your problems with snpEff, I would recommend the VEP tool (variant effect predictor): It is very versatile and relatively easy to use. For standard genomes (human, etc) is my preferred choice.

ADD REPLYlink written 4.6 years ago by abascalfederico1.1k
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: 2290 users visited in the last hour