# HG changeset patch # User modencode-dcc # Date 1358793181 18000 # Node ID 1ceeba7b95a498f9216c6cc7371f1da59a75e510 # Parent 7a51d316635a3a5bad59e20586bcd72ee95ceff5 Uploaded diff -r 7a51d316635a -r 1ceeba7b95a4 bamedit.py --- a/bamedit.py Fri Jan 18 15:55:46 2013 -0500 +++ b/bamedit.py Mon Jan 21 13:33:01 2013 -0500 @@ -12,6 +12,7 @@ def main(): options = simplejson.load( open( sys.argv[1] ) ) + script_path = sys.argv[2] #experiment_name = '_'.join( options['bamout'] ) @@ -20,9 +21,9 @@ if('input3' in options): cmdline = "samtools merge %s %s %s %s" % ( options['bamout'], options['input1'], options['input2'], options['input3'] ) elif (options['action'] == "split"): - cmdline = "bash /mnt/galaxyTools/galaxy-central/tools/modENCODE_DCC_tools/bamedit/split.sh %s %s %s" % ( options['bamout'], options['bamout2'], options['input1'] ) + cmdline = "bash %s/split.sh %s %s %s" % ( script_path, options['bamout'], options['bamout2'], options['input1'] ) elif (options['action'] == "pileup"): - cmdline = "perl /mnt/galaxyTools/galaxy-central/tools/modENCODE_DCC_tools/bamedit/pileup.pl %s %s %s %s %s" % ( options['input1'], options['input2'], options['bamout'], options['bamname'], options['refname'] ) + cmdline = "perl %s/pileup.pl %s %s %s %s %s" % ( script_path, options['input1'], options['input2'], options['bamout'], options['bamname'], options['refname'] ) elif (options['action'] == "filter"): cmdline = "samtools view -q %s %s -bo %s" % ( options['quality'], options['input1'], options['bamout'] )