# HG changeset patch # User modencode-dcc # Date 1358807656 18000 # Node ID 3ef6fc38df62f13ffa9f4bcf7a765864d01dd149 # Parent 8fe679d31b82d8b7f524ed27821fc53c0f6f5b35 Uploaded diff -r 8fe679d31b82 -r 3ef6fc38df62 macs_wrapper.py --- a/macs_wrapper.py Mon Jan 21 17:28:42 2013 -0500 +++ b/macs_wrapper.py Mon Jan 21 17:34:16 2013 -0500 @@ -50,7 +50,7 @@ output_peaks = sys.argv[5] experiment_name = '_'.join( options['experiment_name'].split() ) #save experiment name here, it will be used by macs for filenames (gzip of wig files will fail with spaces - macs doesn't properly escape them)..need to replace all whitespace, split makes this easier - cmdline = "macs -t %s" % ",".join( options['input_chipseq'] ) + cmdline = "macs14 -t %s" % ",".join( options['input_chipseq'] ) if options['input_control']: cmdline = "%s -c %s" % ( cmdline, ",".join( options['input_control'] ) ) cmdline = "%s --format='%s' --name='%s' --gsize='%s' --tsize='%s' --bw='%s' --pvalue='%s' --mfold='%s' %s --lambdaset='%s' %s" % ( cmdline, options['format'], experiment_name, options['gsize'], options['tsize'], options['bw'], options['pvalue'], options['mfold'], options['nolambda'], options['lambdaset'], options['futurefdr'] )