Does anyone have any experience running snpEff on Galaxy?

We installed snpEff from the toolshed and now have the latest changeset (c052639fa666). We modified the config files to correct the pathnames but now we're running into a series of errors when we try to run 'Compute SNP effect' with a file using the hg19 dataset

The first error was:

java.lang.RuntimeException: Property: 'hg37.genome' not found
	at ca.mcgill.mcb.pcingola.interval.Genome.<init>(Genome.java:78)
	at ca.mcgill.mcb.pcingola.snpEffect.Config.readGenomeConfig(Config.java:334)
	at ca.mcgill.mcb.pcingola.snpEffect.Config.read(Config.java:301)
	at ca.mcgill.mcb.pcingola.snpEffect.Config.<init>(Config.java:67)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:582)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:261)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:77
This first error was confusing because the input file was for hg19, so it's not clear why there was a reference to hg37.

Next, we then edited snpEff.config and added the line

hg37.genome : Homo_sapiens

with the other hg37.XX : Homo_sapiens keypairs.

This led to the second error:

ERROR: Cannot read file '/home/galaxy//snpEff/data//hg37/snpEffectPredictor.bin'.
	You can try to download the database by running the following command:
		java -jar snpEff.jar download hg37

java.lang.RuntimeException: java.io.FileNotFoundException: /home/galaxy/snpEff/data/hg37/snpEffectPredictor.bin (No such file or directory)
	at ca.mcgill.mcb.pcingola.snpEffect.SnpEffectPredictor.load(SnpEffectPredictor.java:58)
	at ca.mcgill.mcb.pcingola.snpEffect.Config.loadSnpEffectPredictor(Config.java:240)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:595)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:261)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:77)
Caused by: java.io.FileNotFoundException: /home/galaxy/snpEff/data/hg37/snpEffectPredictor.bin (No such file or directory)
	at java.io.FileInputStream.open(Native Method)
	at java.io.FileInputStream.<init>(FileInputStream.java:137)
	at java.io.FileInputStream.<init>(FileInputStream.java:96)
	at ca.mcgill.mcb.pcingola.util.Gpr.readFileSerializedGzThrow(Gpr.java:468)
	at ca.mcgill.mcb.pcingola.snpEffect.SnpEffectPredictor.load(SnpEffectPredictor.java:55)
	... 4 more
The problem now is that is is looking for the hg37 directory, which doesn't exist  but GRCh37.65 does exist. So we then change the directory name to hg37 and reran, resulting in the third error:

VcfFileIterator.parseVcfLine(276):	Fatal error reading file '/mnt/galaxyData/files/000/dataset_432.dat' (line: 1):
1	2338204	2338204	G/T	+
java.lang.RuntimeException: java.lang.RuntimeException: WARNING: Unkown IUB code for SNP '+'
	at ca.mcgill.mcb.pcingola.fileIterator.VcfFileIterator.parseVcfLine(VcfFileIterator.java:277)
	at ca.mcgill.mcb.pcingola.fileIterator.VcfFileIterator.readNext(VcfFileIterator.java:331)
	at ca.mcgill.mcb.pcingola.fileIterator.VcfFileIterator.readNext(VcfFileIterator.java:1)
	at ca.mcgill.mcb.pcingola.fileIterator.FileIterator.hasNext(FileIterator.java:69)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.iterateVcf(SnpEffCmdEff.java:215)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.runAnalysis(SnpEffCmdEff.java:698)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:648)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:261)
	at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:77)
Caused by: java.lang.RuntimeException: WARNING: Unkown IUB code for SNP '+'
	at ca.mcgill.mcb.pcingola.vcf.VcfEntry.parseAlts(VcfEntry.java:530)
	at ca.mcgill.mcb.pcingola.vcf.VcfEntry.parse(VcfEntry.java:445)
	at ca.mcgill.mcb.pcingola.vcf.VcfEntry.<init>(VcfEntry.java:66)
	at ca.mcgill.mcb.pcingola.fileIterator.VcfFileIterator.parseVcfLine(VcfFileIterator.java:274)
	... 8 more

I'm going to assume that this is all related to the first issue, but I'm not sure how to fix the issue. 

Cheers!
Rob


Rob Leclerc, PhD
P: (US) +1-(917)-873-3037
P: (Shanghai) +86-1-(861)-612-5469
Personal Email: rob.leclerc@aya.yale.edu