Hi,
I' am running galaxy locally and I have installed Bismark and its dependencies from the Galaxy main tool shed. I receive multiple errors when using a small test dataset.
(1) My first error occurs when I try running Bismark bisulfite mapper (bowtie2):
Error indexing reference sequence
Which I can fix by editing the env.sh file in the dependencies directory. I change:
(2) My second error occurs when I try running Bismark methylation extractor:
usage: bismark_methylation_extractor.py [-h] [--infile INFILE] [--single-end]
[--paired-end]
[--report-file REPORT_FILE]
[--comprehensive] [--merge-non-cpg]
[--no-overlap] [--compress COMPRESS]
[--ignore-bps IGNORE_BPS]
[--cpg_ot CPG_OT] [--chg_ot CHG_OT]
[--chh_ot CHH_OT]
[--cpg_ctot CPG_CTOT]
[--chg_ctot CHG_CTOT]
[--chh_ctot CHH_CTOT]
[--cpg_ob CPG_OB]
[--chg_ob CHG_OB]
[--chh_ob CHH_OB]
[--cpg_ctob CPG_CTOB]
[--chg_ctob CHG_CTOB]
[--chh_ctob CHH_CTOB]
[--cpg_context CPG_CONTEXT]
[--chg_context CHG_CONTEXT]
[--chh_context CHH_CONTEXT]
[--non_cpg_context NON_CPG_CONTEXT]
[--non_cpg_context_ot NON_CPG_CONTEXT_OT]
[--non_cpg_context_ctot NON_CPG_CONTEXT_CTOT]
[--non_cpg_context_ob NON_CPG_CONTEXT_OB]
[--non_cpg_context_ctob NON_CPG_CONTEXT_CTOB]
Seems to me the python wrapper lacks the code to tell galaxy where the bismark_methylation_extractor executable is. I can go around the error by editing bismark_methylation_extractor.py. I add the following line:
parser.add_argument( '--bismark_path' )
Then on line 72 I change:
cmd = 'bismark_methylation_extractor --no_header -o %s %s %s'
To:
cmd = args.bismark_path + '/bismark_methylation_extractor --no_header -o %s %s %s'
I suspect I' am not doing something correctly. Let me know if anything is unclear or if you need more information. Thanks for any help!
Charlie