Mercurial > repos > modencode-dcc > macs

changeset 6:3ef6fc38df62 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author modencode-dcc
date Mon, 21 Jan 2013 17:34:16 -0500
parents 8fe679d31b82
children 613f8c6c539a
files macs_wrapper.py
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- 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'] )