Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Variant calling step failing when using process scratch #3

Closed
pditommaso opened this issue Oct 9, 2017 · 1 comment
Closed

Variant calling step failing when using process scratch #3

pditommaso opened this issue Oct 9, 2017 · 1 comment

Comments

@pditommaso
Copy link
Collaborator

When using local scratch folder, the step 5_rnaseq_call_variants returns an error.

This happens because the genome.dict input file contains a reference to a file created in temporary folder not accessible to the task, for example:

# cat genome.dict 
@HD	VN:1.5
@SQ	SN:chr22	LN:51304566	M5:a718acaa6135fdca8357d5bfe94211dd	UR:file:/tmp/nxf.Mzz6eisI1J/genome.fa
@pditommaso
Copy link
Collaborator Author

pditommaso commented Oct 9, 2017

Tried both:

sed -i 's@UR:file:.*genome.fa@UR:file:genome.fa@g' genome.dict

and

sed -i 's@UR:file:.*genome.fa@UR:file:'"$PWD"'/genome.fa@g' genome.dict

But still reporting:

##### ERROR
##### ERROR MESSAGE: Invalid command line: Cannot process the provided BAM/CRAM file(s) because they were not indexed.  The GATK does offer limited processing of unindexed BAM/CRAMs in --unsafe mode, but this feature is unsupported -- use it at your own risk!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant