I'm working on making my first workflow in Galaxy, using a
local server. A high level overview of the workflow is:
1) take a fasta file as input with dbkey set to organism1
2) blast input fasta file against organism2 database
3) manipulate the blast output into bed format
4) extract genomic DNA from organism2 based on the bed data
I can't find a way to set
the dbkey on the bed data to organism2. This results in the
extract genomic DNA step failing, as it seems to default to
the only dbkey that exists in my workflow, the dbkey from the
input fasta file. (The names and intervals are not present in
the organism1 genome, and thus the extract step fails.)
I have found that if I set
my the dbkey of my input fasta file to be organism2 (even
though the sequences contained in it are actually from
organism1), then my workflow successfully completes with the
desired output. But this is a bit of a hack. I'd like to share
my workflow with some coworkers and make it as easy to use as
possible. Therefore I don't want to have to explain this hack
to them.
Is it possible to set the
dbkey within a workflow?
thanks in advance for your
help!
-Will