This is a test version of Biostars. For the public version, visit https://www.biostars.org.
snpEff build database

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

0 answers

No answers yet.

Log in to answer this question.