[fc35802]

It worked

Thanks a lot for your help

java -Xmx4G -jar /GenoStorage/Software/snpEff_2_0_5d/snpEff.jar eff -c /GenoStorage/Software/snpEff_2_0_5d/snpEff.config -v -onlyCoding true -i vcf -o txt hg19 varSample02187ATESTE.flt.vcf > snpEff_sample02187A.txt
00:00:00.000 Reading configuration file '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
00:00:00.226 done
00:00:00.226 Reading database for genome 'hg19' (this might take a while)
00:00:13.202 done
00:00:13.203 Building interval forest
00:00:27.282 done.
00:00:27.282 Predicting variants
00:00:33.278 Creating summary file: snpEff_summary.html
00:00:34.303 Creating genes file: snpEff_genes.txt
00:00:39.008 done.
00:00:39.010 Finishing up

[gringer]

I tried to do that in the snpEff.config file but i won't let me save

This is the config file that was owned by root that you enabled read access to? If you want to modify that, you'll need to enable write access globally as well, which means that anyone else using your computer system can modify that file (i.e. a bad idea).

You could also set the ownership of that config file to your user name, which would mean only you (and root) could modify it, or possibly match the ownership of the data directory, then set group permissions for write access (so people in that group can modify the config file):

chown pcadmin:geneusers <location>/snpEff.config
chmod ug+w <location>/snpEff.config

[fc35802]

I tried to do that in the snpEff.config file but i won't let me save

[gringer]

That's a different error. You now have another file not found error because it's looking for '/GenoStorage/Software/snpEff_2_0_5d/./data/GRCh37.65/snpEffectPredictor.bin' rather than '/GenoStorage/BasesDados/GRCh37.65/snpEffectPredictor.bin'. You may need to change the data_dir parameter in the config file, because it's looking in the same directory that you mentioned your config file was in.

http://snpeff.sourceforge.net/faq.html#Can_you_explain_the_configuration_file?

[fc35802]

It doesn't work either

java -Xmx4G -jar /GenoStorage/Software/snpEff_2_0_5d/snpEff.jar eff -c /GenoStorage/Software/snpEff_2_0_5d/snpEff.config -v -onlyCoding true -i vcf -o txt GRCh37.65 varSample02187ATESTE.flt.vcf > snpEff_sample02187A.txt
00:00:00.000 Reading configuration file '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
00:00:00.225 done
00:00:00.225 Reading database for genome 'GRCh37.65' (this might take a while)
java.lang.RuntimeException: java.io.FileNotFoundException: /GenoStorage/Software/snpEff_2_0_5d/./data/GRCh37.65/snpEffectPredictor.bin (No such file or directory)
at ca.mcgill.mcb.pcingola.Gpr.readFileSerializedGz(Gpr.java:426)
at ca.mcgill.mcb.pcingola.snpEffect.SnpEffectPredictor.load(SnpEffectPredictor.java:48)
at ca.mcgill.mcb.pcingola.snpEffect.Config.loadSnpEffectPredictor(Config.java:101)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:536)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:244)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:72)
Caused by: java.io.FileNotFoundException: /GenoStorage/Software/snpEff_2_0_5d/./data/GRCh37.65/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.Gpr.readFileSerializedGz(Gpr.java:421)
... 5 more
00:00:00.235 Finishing up

[gringer]

What about 'GRCh37.65', as mentioned in the 'Latest' column on this page:

http://snpeff.sourceforge.net/download.html

[fc35802]

java -Xmx4G -jar /GenoStorage/Software/snpEff_2_0_5d/snpEff.jar eff -c /GenoStorage/Software/snpEff_2_0_5d/snpEff.config -v -onlyCoding true -i vcf -o txt hg37 varSample02187ATESTE.flt.vcf > snpEff_sample02187A.txt
00:00:00.000 Reading configuration file '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
java.lang.RuntimeException: No such genome 'hg37'
at ca.mcgill.mcb.pcingola.snpEffect.Config.<init>(Config.java:48)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:523)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:244)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:72)
00:00:00.225 Finishing up

That would be nice, but didn't work

[gringer]

Disclaimer: I've never used snpEff, so can only guess at these things

Looking at the snpEff manual, it looks like you need the code name for the genome version, rather than a link to a particular file:

java -Xmx4G -jar /GenoStorage/Software/snpEff_2_0_5d/snpEff.jar eff -c /GenoStorage/Software/snpEff_2_0_5d/snpEff.config -v -onlyCoding true -i vcf -o txt hg37

[fc35802]

The problem of snpEff.config is solved, now I have another

java -Xmx4G -jar /GenoStorage/Software/snpEff_2_0_5d/snpEff.jar eff -c /GenoStorage/Software/snpEff_2_0_5d/snpEff.config -v -onlyCoding true -i vcf -o txt /GenoStorage/BasesDados/GRCh37.65 varSample02187ATESTE.flt.vcf > snpEff_sample02187A.txt
00:00:00.000 Reading configuration file '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
java.lang.RuntimeException: No such genome '/GenoStorage/BasesDados/GRCh37.65'
at ca.mcgill.mcb.pcingola.snpEffect.Config.<init>(Config.java:48)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:523)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:244)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:72)
00:00:00.226 Finishing up
[meusebio@IMMGene1 meusebio]$ ls -lah /GenoStorage/BasesDados/GRCh37.65/
total 67M
drwxr-xr-x 2 pcadmin geneusers 43 Feb 17 10:48 .
drwxrwx---. 6 pcadmin geneusers 84 Feb 17 10:53 ..
-rw-r--r-- 1 root root 67M Feb 17 10:48 snpEffectPredictor.bin

[fc35802]

I don't have permission to do that, but I already asked for it.

Now I have to wait to see if it worked

[gringer]

I think it didn't fail

ls -lah  '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
-rw------- 1 root root 45K Jan 20 02:10 /GenoStorage/Software/snpEff_2_0_5d/snpEff.config

This file can only be read by root (oh, and thanks for making it long format so that I noticed this problem). If you're running Java as root, that's a bad idea, but I'm guessing based on the error that this is being run as a more normal user. You can change the file modes to permit reading with chmod:

chmod go+r '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'

[because this file is owned by root, you'll only be able to do that as root, or as a suitably privileged user]

[fc35802]

I think it didn't fail

ls -lah '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
-rw------- 1 root root 45K Jan 20 02:10 /GenoStorage/Software/snpEff_2_0_5d/snpEff.config

[gringer]

The error mentions that java is unable to find a file at that location. You can check to make sure the file doesn't exist with ls:

ls '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'

If ls fails, then the Java command will also fail. Check the capitalisation of names, and make sure that your 'GenoStorage' directory is actually based off the root of the file system (which seems like a strange place to put it).

[fc35802]

I tried, but it didin't work

java -Xmx4G -jar /GenoStorage/Software/snpEff_2_0_5d/snpEff.jar eff -c /GenoStorage/Software/snpEff_2_0_5d/snpEff.config -v -onlyCoding true -i vcf -o txt /GenoStorage/BasesDados/GRCh37.64 varSample02187ATESTE.flt.vcf > snpEff_sample02187A.txt
00:00:00.000 Reading configuration file '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
java.lang.RuntimeException: Cannot find config file '/GenoStorage/Software/snpEff_2_0_5d/snpEff.config'
at ca.mcgill.mcb.pcingola.snpEffect.Config.read(Config.java:117)
at ca.mcgill.mcb.pcingola.snpEffect.Config.<init>(Config.java:46)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEffCmdEff.run(SnpEffCmdEff.java:523)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.run(SnpEff.java:244)
at ca.mcgill.mcb.pcingola.snpEffect.commandLine.SnpEff.main(SnpEff.java:72)
00:00:00.008 Finishing up

[mehc]

Yes try it out