picard CreateSequenceDictionary errors when trying to make a .dict file
1
0
Entering edit mode
8 weeks ago

I'm trying to use picard to index my reference file, but I get the error below. I can't seem to find a solution online by searching the errors.

REF=~/refs/hg38/chr4.fa
ACC=chr4
picard CreateSequenceDictionary REFERENCE="$REF" OUTPUT="~/refs/hg38/${ACC}.dict"

Throws this error:

[Fri Feb 23 11:32:26 GMT 2024] picard.sam.CreateSequenceDictionary done. Elapsed time: 0.01 minutes.
Runtime.totalMemory()=538968064
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" picard.PicardException: File /Users/michaelflower/~/refs/hg38/chr4.dict not found
    at picard.sam.CreateSequenceDictionary.doWork(CreateSequenceDictionary.java:221)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:305)
    at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:103)
    at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:113)
picard dict • 238 views
ADD COMMENT
0
Entering edit mode
8 weeks ago

Try using a full path, I think the tilde ~ is confusing the tool.

#instead of this
picard CreateSequenceDictionary REFERENCE="$REF" OUTPUT="~/refs/hg38/${ACC}.dict"
#use
picard CreateSequenceDictionary REFERENCE="$REF" OUTPUT="/my/home/refs/hg38/${ACC}.dict"
ADD COMMENT

Login before adding your answer.

Traffic: 1277 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6