Mercurial > repos > zzhou > bamedit
diff bamedit.py @ 81:500d54c2edf9 draft
Uploaded
| author | zzhou |
|---|---|
| date | Mon, 26 Nov 2012 15:31:07 -0500 |
| parents | 0d3c69699f0a |
| children | 784e546849d6 |
line wrap: on
line diff
--- a/bamedit.py Mon Nov 26 15:31:01 2012 -0500 +++ b/bamedit.py Mon Nov 26 15:31:07 2012 -0500 @@ -12,18 +12,19 @@ def main(): options = simplejson.load( open( sys.argv[1] ) ) + script_path = sys.argv[2] + print script_path - testing = sys.argv - print testing if(options['action'] == "merge"): cmdline = "samtools merge %s %s %s" % ( options['bamout'], options['input1'], options['input2'] ) 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/bamedit/split.sh %s %s %s" % ( options['bamout'], options['bamout2'], options['input1'] ) - cmdline = "bash 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/bamedit/pileup.pl %s %s %s %s %s" % ( options['input1'], options['input2'], options['bamout'], options['bamname'], options['refname'] ) + #cmdline = "perl /mnt/galaxyTools/galaxy-central/tools/modENCODE_DCC/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'] )
