Question: snpEff build database
0
gravatar for cgenome
2.7 years ago by
cgenome0
cgenome0 wrote:

Hello, I'm trying to use the snpEff tool lately but it has been more dificult than I thought.

build snpEff databases

mkdir /usr/local/Cellar/snpeff/4.2/share/snpeff/data mkdir /usr/local/Cellar/snpeff/4.2/share/snpeff/data/SBW25_pQBR103 cat ./ref/AM181176.embl ./ref/AM235768.embl > /usr/local/Cellar/snpeff/4.2/share/snpeff/data/SBW25_pQBR103/genes.embl

edit snpEff.config file according to documentation at the directory /usr/local/Cellar/snpeff/4.2/share/snpeff

Pseudomonas_fluorescens_with_pQBR103

SBW25_pQBR103.genome : Pseudomonas_fluorescens_with_pQBR103 SBW25_pQBR103.chromosomes : AM181176, AM235768 SBW25_pQBR103.AM181176.codonTable : Bacterial_and_Plant_Plastid SBW25_pQBR103.AM235768.codonTable : Bacterial_and_Plant_Plastid

then run the command

snpEff build -embl -v SBW25_pQBR103

And everything is fine so far. However when I'm trying to run this command: java -jar ~/bin/snpEff/snpEff.jar -no-upstream -no-downstream ./ref/SBW25_pQBR103 ./raw_vcf/SBW25_hap.vcf > ./snpEff/SBW25_hap.snpEff.vcf

I still get the error: java.lang.RuntimeException: Property: './ref/SBW25_pQBR103.genome' not found at ca.mcgill.mcb.pcingola.interval.Genome.<init>(Genome.java:106) at ca.mcgill.mcb.pcingola.snpEffect.Config.readGenomeConfig(Config.java:619) at ca.mcgill.mcb.pcingola.snpEffect.Config.readConfig(Config.java:587) at ca.mcgill.mcb.pcingola.snpEffect.Config.init(Config.java:434) at ca.mcgill.mcb.pcingola.snpEffect.Config.<init>(Config.java:110) at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.loadConfig(SnpEff.java:273) at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:955) at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:939) at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:978) at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:136)

Any thoughts what could be wrong? Do I need to create this file (SBW25_pQBR103.genome) ? Or is somehow saved somewhere? Because I couldn't find it anywhere except from the fact that it's written at the configuration file.

Thanks.

snpeff • 1.3k views
ADD COMMENTlink written 2.7 years ago by cgenome0
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: 798 users visited in the last hour