Mercurial > repos > pjbriggs > macs21

changeset 0:fdad0c8c0957 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded initial version to test toolshed.
author pjbriggs
date Wed, 21 Jan 2015 11:07:37 -0500
parents
children d0986d2be693
files macs21_wrapper.py macs21_wrapper.xml tool_dependencies.xml
diffstat 3 files changed, 629 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/macs21_wrapper.py	Wed Jan 21 11:07:37 2015 -0500
@@ -0,0 +1,208 @@
+#purpose: macs2 python wrapper
+#author: Ziru Zhou
+#date: November, 2012
+
+import sys, subprocess, tempfile, shutil, glob, os, os.path, gzip
+from galaxy import eggs
+import pkg_resources
+pkg_resources.require( "simplejson" )
+import simplejson
+
+CHUNK_SIZE = 1024
+
+#==========================================================================================
+#functions
+#==========================================================================================
+def gunzip_cat_glob_path( glob_path, target_filename, delete = False ):
+    out = open( target_filename, 'wb' )
+    for filename in glob.glob( glob_path ):
+        fh = gzip.open( filename, 'rb' )
+        while True:
+            data = fh.read( CHUNK_SIZE )
+            if data:
+                out.write( data )
+            else:
+                break
+        fh.close()
+        if delete:
+            os.unlink( filename )
+    out.close()
+
+def xls_to_interval( xls_file, interval_file, header = None ):
+    out = open( interval_file, 'wb' )
+    if header:
+        out.write( '#%s\n' % header )
+    wrote_header = False
+    #From macs readme: Coordinates in XLS is 1-based which is different with BED format.
+    #PJB: updated to keep original 'start' coordinate i.e. do not convert to BED
+    for line in open( xls_file ):
+        #keep all existing comment lines
+        if line.startswith( '#' ):
+            out.write( line )
+	#added for macs2 since there is an extra newline 
+        ##elif line.startswith( '\n' ):
+        ##    out.write( line )
+        elif not wrote_header:
+            out.write( '#%s' % line )
+	    print line
+            wrote_header = True
+        else:
+            fields = line.split( '\t' )
+            if len( fields ) > 1:
+                try:
+                    # Try to convert 'start' to int and shift
+                    ##fields[1] = str( int( fields[1] ) - 1 )
+                    # PJB don't correct, just test to see if it's an integer
+                    int( fields[1] )
+                except ValueError:
+                    # Integer conversion failed so comment out
+                    # "bad" line instead
+                    fields[0] = "#%s" % fields[0]
+            out.write( '\t'.join( fields ) )
+    out.close()
+
+#==========================================================================================
+#main
+#==========================================================================================
+def main():
+    #take in options file and output file names
+    options = simplejson.load( open( sys.argv[1] ) )
+    outputs = simplejson.load( open( sys.argv[2] ) )
+
+    #=================================================================================
+    #parse options and execute macs2
+    #=================================================================================
+    #default inputs that are in every major command
+    experiment_name = '_'.join( options['experiment_name'].split() ) #save experiment name here, it will be used by macs for some file names
+    cmdline = "macs2 %s -t %s" % ( options['command'], ",".join( options['input_chipseq'] ) )
+    if options['input_control']:
+        cmdline = "%s -c %s" % ( cmdline, ",".join( options['input_control'] ) )
+
+    #=================================================================================
+    if (options['command'] == "callpeak"):
+    	output_summits_bed = outputs['output_summits_bed_file']
+    	output_extra_html = outputs['output_extra_file']
+    	output_extra_path = outputs['output_extra_file_path']
+    	output_peaks =  outputs['output_peaks_file']
+   	output_narrowpeaks = outputs['output_narrowpeaks_file']
+   	output_broadpeaks = outputs['output_broadpeaks_file']
+   	output_gappedpeaks = outputs['output_gappedpeaks_file']
+	output_xls_to_interval_peaks_file = outputs['output_xls_to_interval_peaks_file']
+        output_treat_pileup = outputs['output_treat_pileup_file']
+        output_lambda_bedgraph = outputs['output_lambda_bedgraph_file']
+
+        cmdline = "%s --format='%s' --name='%s' --gsize='%s' --bw='%s'" % ( cmdline, options['format'], experiment_name, options['gsize'], options['bw'] )
+        if (options['broad'] == 'broad'):
+            cmdline = "%s --broad --broad-cutoff='%s'" % ( cmdline, options['broad_cutoff'] )
+	if 'pvalue' in options:
+            cmdline = "%s --pvalue='%s'" % ( cmdline, options['pvalue'] )
+	elif 'qvalue' in options:
+            cmdline = "%s --qvalue='%s'" % ( cmdline, options['qvalue'] )
+        cmdline = "%s --mfold %s %s %s %s %s %s --keep-dup %s" % (cmdline, options['mfoldlo'], options['mfoldhi'], options['nolambda'], options['bdg'], options['spmr'], options['call_summits'], options['keep_dup'] )
+		
+	if 'nomodel' in options:
+        	cmdline = "%s --nomodel --extsize='%s'" % ( cmdline, options['nomodel'] )
+    #=================================================================================
+    if (options['command'] == "bdgcmp"):
+	output_bdgcmp = outputs['output_bdgcmp_file']
+
+	cmdline = "%s -m %s -p %s -o bdgcmp_out.bdg" % ( cmdline, options['m'], options['pseudocount'] )
+    #=================================================================================
+
+    tmp_dir = tempfile.mkdtemp() #macs makes very messy output, need to contain it into a temp dir, then provide to user
+    stderr_name = tempfile.NamedTemporaryFile().name # redirect stderr here, macs provides lots of info via stderr, make it into a report
+    proc = subprocess.Popen( args=cmdline, shell=True, cwd=tmp_dir, stderr=open( stderr_name, 'wb' ) )
+    proc.wait()
+    #We don't want to set tool run to error state if only warnings or info, e.g. mfold could be decreased to improve model, but let user view macs log
+    #Do not terminate if error code, allow dataset (e.g. log) creation and cleanup
+    if proc.returncode:
+        stderr_f = open( stderr_name )
+        while True:
+            chunk = stderr_f.read( CHUNK_SIZE )
+            if not chunk:
+                stderr_f.close()
+                break
+            sys.stderr.write( chunk )
+    
+    #=================================================================================
+    #copy files created by macs2 to appripriate directory with the provided names
+    #=================================================================================
+
+    #=================================================================================
+    #move files generated by callpeak command 
+    if (options['command'] == "callpeak"):
+    	#run R to create pdf from model script
+    	if os.path.exists( os.path.join( tmp_dir, "%s_model.r" % experiment_name ) ):
+     	   	cmdline = 'R --vanilla --slave < "%s_model.r" > "%s_model.r.log"' % ( experiment_name, experiment_name )
+     	   	proc = subprocess.Popen( args=cmdline, shell=True, cwd=tmp_dir )
+		proc.wait()
+        
+    	#move bed out to proper output file
+    	created_bed_name =  os.path.join( tmp_dir, "%s_summits.bed" % experiment_name )
+    	if os.path.exists( created_bed_name ):
+        	shutil.move( created_bed_name, output_summits_bed )
+
+    	#OICR peak_xls file
+    	created_peak_xls_file =  os.path.join( tmp_dir, "%s_peaks.xls" % experiment_name )
+    	if os.path.exists( created_peak_xls_file ):
+        	# shutil.copy( created_peak_xls_file, os.path.join ( "/mnt/galaxyData/tmp/", "%s_peaks.xls" % ( os.path.basename(output_extra_path) )))
+		shutil.copyfile( created_peak_xls_file, output_peaks )   
+
+    	#peaks.encodepeaks (narrowpeaks) file
+    	created_narrowpeak_file = os.path.join (tmp_dir, "%s_peaks.narrowPeak" % experiment_name )
+    	if os.path.exists( created_narrowpeak_file ):
+		shutil.move (created_narrowpeak_file, output_narrowpeaks )
+
+        #peaks.broadpeaks file
+    	created_broadpeak_file = os.path.join (tmp_dir, "%s_peaks.broadPeak" % experiment_name )
+    	if os.path.exists( created_broadpeak_file ):
+		shutil.move (created_broadpeak_file, output_broadpeaks )
+
+        #peaks.gappedpeaks file
+    	created_gappedpeak_file = os.path.join (tmp_dir, "%s_peaks.gappedPeak" % experiment_name )
+    	if os.path.exists( created_gappedpeak_file ):
+		shutil.move (created_gappedpeak_file, output_gappedpeaks )
+
+ 	#parse xls files to interval files as needed
+   	#if 'xls_to_interval' in options:
+    	if (options['xls_to_interval'] == "True"):
+        	create_peak_xls_file = os.path.join( tmp_dir, '%s_peaks.xls' % experiment_name )
+        	if os.path.exists( create_peak_xls_file ):
+            		xls_to_interval( create_peak_xls_file, output_xls_to_interval_peaks_file, header = 'peaks file' )
+
+        #bedgraph bedgraph files
+        if (options['bdg']):
+            created_treat_pileup_file = os.path.join (tmp_dir, "%s_treat_pileup.bdg" % experiment_name )
+            if os.path.exists( created_treat_pileup_file ):
+		shutil.move (created_treat_pileup_file, output_treat_pileup )
+            created_lambda_bedgraph_file = os.path.join (tmp_dir, "%s_control_lambda.bdg" % experiment_name )
+            if os.path.exists( created_lambda_bedgraph_file ):
+		shutil.move (created_lambda_bedgraph_file, output_lambda_bedgraph )
+
+    	#move all remaining files to extra files path of html file output to allow user download
+    	out_html = open( output_extra_html, 'wb' )
+    	out_html.write( '<html><head><title>Additional output created by MACS (%s)</title></head><body><h3>Additional Files:</h3><p><ul>\n' % experiment_name )
+    	os.mkdir( output_extra_path )
+    	for filename in sorted( os.listdir( tmp_dir ) ):
+    		shutil.move( os.path.join( tmp_dir, filename ), os.path.join( output_extra_path, filename ) )
+        	out_html.write( '<li><a href="%s">%s</a></li>\n' % ( filename, filename ) )
+		#out_html.write( '<li><a href="%s">%s</a>peakxls %s SomethingDifferent tmp_dir %s path %s exp_name %s</li>\n' % ( created_peak_xls_file, filename, filename, tmp_dir, output_extra_path, experiment_name ) )
+    	out_html.write( '</ul></p>\n' )
+    	out_html.write( '<h3>Messages from MACS:</h3>\n<p><pre>%s</pre></p>\n' % open( stderr_name, 'rb' ).read() )
+    	out_html.write( '</body></html>\n' )
+    	out_html.close()
+
+    #=================================================================================
+    #move files generated by bdgcmp command
+    if (options['command'] == "bdgcmp"):
+    	created_bdgcmp_file = os.path.join (tmp_dir, "bdgcmp_out.bdg" )
+    	if os.path.exists( created_bdgcmp_file ):
+		shutil.move (created_bdgcmp_file, output_bdgcmp )
+     
+    #=================================================================================    
+    #cleanup
+    #=================================================================================    
+    os.unlink( stderr_name )
+    os.rmdir( tmp_dir )
+
+if __name__ == "__main__": main()
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/macs21_wrapper.xml	Wed Jan 21 11:07:37 2015 -0500
@@ -0,0 +1,387 @@
+<tool id="fls_modencode_peakcalling_macs2.1" name="MACS2.1.0" version="2">
+  <requirements>
+    <requirement type="package" version="2.7">python</requirement>
+    <requirement type="package" version="1.8.1">numpy</requirement>
+    <requirement type="package" version="2.1.0.20140616">macs2</requirement>
+  </requirements>
+  <description>Model-based Analysis of ChIP-Seq [golem]</description>
+  <command interpreter="python">macs21_wrapper.py $options_file $outputs_file</command>
+  <inputs>
+    <!--experiment name and option of selecting paired or single end will always be present-->
+    <param name="experiment_name" type="text" value="MACS2.1.0 in Galaxy" size="50"
+	   label="Experiment Name"/>
+    <!--select one of the 7 major commands offered by macs2-->
+    <conditional name="major_command">
+      <param name="major_command_selector" type="select" label="Select action to be performed">
+	<option value="callpeak">Peak Calling</option>
+	<!--<option value="filterdup">filterdup</option>
+	<option value="randsample">randsample</option>-->
+	<option value="bdgcmp">Compare .bdg Files</option>
+	<!--<option value="bdgdiff">bdgdiff</option>
+	<option value="bdgpeakcall">bdgpeakcall</option>
+	<option value="bdgbroadcall">bdgbroadcall</option>-->
+      </param>
+      <!--callpeak option of macs2-->
+      <when value="callpeak">
+	<!--choose 'broad' or 'narrow' regions-->
+	<conditional name="broad_options">
+	  <param name="broad_regions" type="select" label="Type of region to call"
+		 help="Broad regions are formed by linking nearby enriched regions">
+	    <option value="" selected="true">Narrow regions</option>
+	    <option value="broad">Broad regions</option>
+	  </param>
+	  <when value="broad">
+	    <param name="broad_cutoff" type="float"
+		   label="Cutoff for broad regions"
+		   value="0.1" help="default: 0.1 (--broad-cutoff)"/>
+	  </when>
+	</conditional>
+	<!--may need to add a few more formats at later time-->
+        <param name="input_chipseq_file1" type="data" format="bed,sam,bam"
+	       label="ChIP-seq read file" />
+        <param name="input_control_file1" type="data" format="bed,sam,bam" optional="True"
+	       label="ChIP-seq control read file" />
+	<conditional name="genome_size">
+	  <param name="gsize" type="select" label="Effective genome size"
+		 help="Either pre-defined (for common organisms), or user-defined (--gsize)">
+	    <option value="hs" selected="true">Human (2.7e9)</option>
+	    <option value="mm">Mouse (1.87e9)</option>
+	    <option value="ce">C. elegans (9e7)</option>
+	    <option value="dm">Fruitfly (1.2e8)</option>
+	    <option value="">User-defined</option>
+	  </param>
+	  <when value="">
+	    <!-- User-defined effective genome size -->
+	    <param name="user_defined_gsize" type="float" value=""
+		   label="Enter effective genome size (number of bases)"
+		   help="e.g. '1.0e+9' or '1000000000'" />
+	  </when>
+	</conditional>
+	<param name="bw" type="integer" label="Band width" value="300" help="(--bw)"/>
+	<param name="xls_to_interval" label="Include XLS file from MACS"
+	       type="boolean" truevalue="True" falsevalue="False" checked="True"
+	       help="MACS2 XLS file will be output to the history in 'interval' format (suitable for subsequent analysis in Galaxy). Note that start positions are 1-based."/>
+
+	<conditional name="bdg_options">
+	  <param name="bdg"
+		 label="Save fragment pileup, control lambda, -log10pvalue/qvalue in bedGraph"
+		 type="boolean" truevalue="-B" falsevalue="" checked="False" />
+	  <when value="-B">
+	    <param name="spmr"
+		   type="boolean" truevalue="--SPMR" falsevalue="" checked="False"
+		   label="Save signal per million reads for fragment pileup profiles"
+		   help="(--SPMR)" />
+	  </when>
+	  <when value="">
+	    <!-- Display nothing -->
+	  </when>
+	</conditional>
+
+	<conditional name="pq_options">
+	  <param name="pq_options_selector" type="select"
+		 label="Select p-value or q-value" help="default uses q-value">
+	    <option value="qvalue">q-value</option>
+	    <option value="pvalue">p-value</option>
+	  </param>
+	  <when value="pvalue">
+	    <param name="pvalue" type="float"
+		   label="p-value cutoff for binding region detection"
+		   value="1e-2" help="default: 1e-2 (--pvalue)"/>
+	  </when>
+	  <when value="qvalue">
+	    <param name="qvalue" type="float"
+		   label="q-value cutoff for binding region detection"
+		   value="0.01" help="default: 0.01 (--qvalue)"/>
+	  </when>
+	</conditional>
+	<conditional name="advanced_options">
+	  <param name="advanced_options_selector" type="select"
+		 label="Display advanced options">
+	    <option value="off">Hide</option>
+	    <option value="on">Display</option>
+	  </param>
+	  <when value="on">
+            <param name="mfoldlo" type="integer"
+		   label="Select the regions with MFOLD high-confidence enrichment ratio against background to build model (lower-limit)"
+		   value="10" help="(--mfold)"/>
+	    <param name="mfoldhi" type="integer"
+		   label="Select the regions with MFOLD high-confidence enrichment ratio against background to build model (upper-limit)"
+		   value="30" help="(--mfold)"/>
+	    <param name="nolambda"
+		   label="Use fixed background lambda as local lambda for every binding region"
+		   type="boolean" truevalue="--nolambda" falsevalue="" checked="False"
+		   help="(--nolambda)"/>
+	    <param name="call_summits"
+		   label="Detect subpeaks within binding region"
+		   type="boolean" truevalue="--call-summits" falsevalue="" checked="False"
+		   help="(--call-summits)"/>
+	    <conditional name="keep_duplicates">
+	      <param name="keep_dup" type="select"
+		     label="Use of duplicate reads">
+		<option value="auto">Automatically calculate maximum number of duplicates to keep (auto)</option>
+		<option value="all">Use all duplicates (all)</option>
+		<option value="" selected="true">Manually specify maxium number of duplicates</option>
+	      </param>
+	      <when value="">
+		<param name="maximum_tags" type="integer" value="1"
+		       label="Maxium number of duplicated tags to keep at each location"/>
+	      </when>
+	    </conditional>
+	  </when>
+	  <when value="off">
+	    <!--display nothing-->
+	  </when>
+	</conditional>
+    	<conditional name="nomodel_type">
+          <param name="nomodel_type_selector" type="select" label="Build Model">
+	   <option value="nomodel">Do not build the shifting model (--nomodel enabled)</option>
+           <option value="create_model" selected="true">Build the shifting model (--nomodel disabled)</option>
+          </param>
+          <when value="nomodel">
+            <param name="extsize" type="integer" label="Arbitrary extension size in bp" value="200" help="Used as fragment size to extend each read towards 3' end (--extsize)"/>
+          </when>
+        </conditional>
+      </when>
+
+      <!--callpeak option of macs2-->
+      <when value="bdgcmp">
+        <param name="input_chipseq_file1" type="data" format="bed,sam,bam"
+	       label="ChIP-seq read file" />
+        <param name="input_control_file1" type="data" format="bed,sam,bam" optional="True"
+	       label="ChIP-seq control read file" />
+	<param name="pseudocount" type="float" label="Set pseudocount" value="0.00001"
+	       help="default: 0.00001 (-p)"/>
+        <conditional name="bdgcmp_options">
+          <param name="bdgcmp_options_selector" type="select"
+		 label="Select action to be performed">
+	    <option value="ppois">ppois</option>
+	    <option value="qpois">qpois</option>
+	    <option value="subtract">subtract</option>
+	    <option value="logFE">logFE</option>
+	    <option value="FE">FE</option>
+	    <option value="logLR">logLR</option>
+          </param>
+	</conditional>
+      </when>
+    </conditional>
+  </inputs>
+
+  <outputs>
+    <!--callpeaks output-->
+    <data name="output_extra_files" format="html"
+	  label="${tool.name}: callpeak on ${on_string} (html report)">
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+    </data>
+    <data name="output_summits_bed_file" format="bed"
+	  label="${tool.name}: callpeak on ${on_string} (summits: bed)">
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+    </data>
+    <data name="output_peaks_file" format="xls"
+	  label="${tool.name}: callpeak on ${on_string} (peaks: xls)">
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+      <filter>major_command['xls_to_interval'] is False</filter>
+    </data>
+    <data name="output_narrowpeaks_file" format="interval"
+	  label="${tool.name}: callpeak on ${on_string} (peaks: narrowPeak)">
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+      <filter>major_command['broad_options']['broad_regions'] == ''</filter>
+    </data>
+    <data name="output_broadpeaks_file" format="interval"
+	  label="${tool.name}: callpeak on ${on_string} (peaks: broadPeak)">
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+      <filter>major_command['broad_options']['broad_regions'] == 'broad'</filter>
+    </data>
+    <data name="output_gappedpeaks_file" format="interval"
+	  label="${tool.name}: callpeak on ${on_string} (peaks: gappedPeak)">
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+      <filter>major_command['broad_options']['broad_regions'] == 'broad'</filter>
+    </data>
+    <data name="output_xls_to_interval_peaks_file" format="interval"
+	  label="${tool.name}: callpeak on ${on_string} (peaks: interval)">
+      <filter>major_command['xls_to_interval'] is True</filter>
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+    </data>
+    <data name="output_treat_pileup_file" format="bedgraph"
+	  label="${tool.name}: callpeak on ${on_string} (treat pileup: bedGraph)">
+      <filter>major_command['bdg_options']['bdg'] is True</filter>
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+    </data>
+    <data name="output_lambda_bedgraph_file" format="bedgraph"
+	  label="${tool.name}: callpeak on ${on_string} (control lambda: bedGraph)">
+      <filter>major_command['bdg_options']['bdg'] is True</filter>
+      <filter>major_command['major_command_selector'] == 'callpeak'</filter>
+    </data>
+    <!--bdgcmp output-->
+    <data name="output_bdgcmp_file" format="bdg"
+	  label="${tool.name}: bdgcmp on ${on_string} (bdg)">
+      <filter>major_command['major_command_selector'] == 'bdgcmp'</filter>
+    </data>
+  </outputs>
+  <configfiles>
+    <configfile name="outputs_file">&lt;%
+import simplejson
+%&gt;
+##=======================================================================================
+#set $__outputs = { 'command':str( $major_command.major_command_selector ) }
+#if str( $major_command.major_command_selector ) == 'callpeak':
+	#set $__outputs['output_summits_bed_file'] = str( $output_summits_bed_file )
+	#set $__outputs['output_extra_file'] = str( $output_extra_files )
+	#set $__outputs['output_extra_file_path'] = str( $output_extra_files.files_path )
+	#set $__outputs['output_peaks_file'] = str( $output_peaks_file )
+	#set $__outputs['output_narrowpeaks_file'] = str( $output_narrowpeaks_file )
+	#set $__outputs['output_broadpeaks_file'] = str( $output_broadpeaks_file )
+	#set $__outputs['output_gappedpeaks_file'] = str( $output_gappedpeaks_file )
+	#set $__outputs['output_xls_to_interval_peaks_file'] = str( $output_xls_to_interval_peaks_file )
+	#set $__outputs['output_treat_pileup_file'] = str( $output_treat_pileup_file )
+	#set $__outputs['output_lambda_bedgraph_file'] = str( $output_lambda_bedgraph_file )
+#end if
+##=======================================================================================
+#if str( $major_command.major_command_selector ) == 'bdgcmp':
+	#set $__outputs['output_bdgcmp_file'] = str( $output_bdgcmp_file )
+#end if
+
+${ simplejson.dumps( __outputs ) }
+    </configfile>
+    <configfile name="options_file">&lt;%
+import simplejson
+%&gt;
+##=======================================================================================
+#set $__options = { 'experiment_name':str( $experiment_name ) }
+##treatment/tag input files and format
+#set $__options['input_chipseq'] = [ str( $major_command.input_chipseq_file1 ) ]
+#set $__options['format'] = $major_command.input_chipseq_file1.extension.upper()
+
+##control/input files
+#set $__options['input_control'] = []
+#if str( $major_command.input_control_file1 ) != 'None':
+	#set $_hole = __options['input_control'].append( str( $major_command.input_control_file1 ) )
+#end if
+
+#if str( $major_command.major_command_selector ) == 'callpeak':
+	#set $__options['command'] = str( "callpeak" )
+	#set $__options['bw'] = str( $major_command.bw )
+	#set $__options['xls_to_interval'] = str( $major_command.xls_to_interval )
+
+	##bdg options
+	#if $major_command.bdg_options.bdg == True:
+		#set $__options['bdg'] = str( "-B" )
+		#set $__options['spmr'] = str( $major_command.bdg_options.spmr )
+	#else:
+		#set $__options['bdg'] = str( "" )
+		#set $__options['spmr'] = str( "" )
+	#end if
+
+	##broad_options
+	#if str( $major_command.broad_options.broad_regions ) == 'broad':
+		#set $__options['broad'] = str( $major_command.broad_options.broad_regions )
+		#set $__options['broad_cutoff'] = str( $major_command.broad_options.broad_cutoff )
+	#else:
+		#set $__options['broad'] = str( "" )
+		#set $__options['broad_cutoff'] = str( "" )
+	#end if
+
+	##genome sizes
+	#if str( $major_command.genome_size.gsize ) == '':
+		#set $__options['gsize'] = int( $major_command.genome_size.user_defined_gsize )
+	#else:
+		#set $__options['gsize'] = str( $major_command.genome_size.gsize )
+	#end if
+
+	##advanced options
+	#if str( $major_command.advanced_options.advanced_options_selector ) == 'on':
+		#set $__options['mfoldlo'] = int( $major_command.advanced_options.mfoldlo )
+		#set $__options['mfoldhi'] = int( $major_command.advanced_options.mfoldhi )
+		#set $__options['nolambda'] = str( $major_command.advanced_options.nolambda )
+		#set $__options['call_summits'] = str( $major_command.advanced_options.call_summits )
+		#if str( $major_command.advanced_options.keep_duplicates.keep_dup ) == '':
+			#set $__options['keep_dup'] = int( $major_command.advanced_options.keep_duplicates.maximum_tags )
+		#else:
+			#set $__options['keep_dup'] = str( $major_command.advanced_options.keep_duplicates.keep_dup )
+		#end if
+	#else:
+		#set $__options['mfoldlo'] = int( "5" )
+		#set $__options['mfoldhi'] = int( "50" )
+		#set $__options['nolambda'] = str( "" )
+		#set $__options['call_summits'] = str( "" )
+		#set $__options['keep_dup'] = int( "1" )
+	#end if
+
+	##enable xls file options
+	##if str( $major_command.xls_to_interval ) == 'create':
+		##set $__options['xls_to_interval'] = { 'peaks_file': str( $output_xls_to_interval_peaks_file ), 'negative_peaks_file': str( $output_xls_to_interval_negative_peaks_file ) }
+	##end if
+	
+	##pq value select options
+	#if str( $major_command.pq_options.pq_options_selector ) == 'qvalue':
+		#set $__options['qvalue'] = str( $major_command.pq_options.qvalue )
+	#else:
+		#set $__options['pvalue'] = str( $major_command.pq_options.pvalue )
+	#end if
+	
+	##model options
+	#if str( $major_command.nomodel_type.nomodel_type_selector ) == 'nomodel':
+		#set $__options['nomodel'] = str( $major_command.nomodel_type.extsize )
+	#end if
+#end if
+##=======================================================================================
+#if str( $major_command.major_command_selector ) == 'bdgcmp':
+	#set $__options['command'] = str( "bdgcmp" )
+	#set $__options['pseudocount'] = float( str( $major_command.pseudocount ) )
+	#set $__options['m'] = str( $major_command.bdgcmp_options.bdgcmp_options_selector )
+#end if
+##=======================================================================================
+
+${ simplejson.dumps( __options ) }
+    </configfile>
+  </configfiles>
+  <tests>
+	<!--none yet for macs2-->
+  </tests>
+  <help>
+
+.. class:: warningmark
+
+**This is a modified version of the standard Galaxy toolshed "MACS2" tool,
+which has been customised for users at the University of Manchester to work
+with MACS 2.1.0.**
+
+It is based on the 16:14f378e35191 revision of the tool at
+
+ *  http://toolshed.g2.bx.psu.edu/view/modencode-dcc/macs2 
+
+------
+
+**What it does**
+
+With the improvement of sequencing techniques, chromatin immunoprecipitation
+followed by high throughput sequencing (ChIP-Seq) is getting popular to study
+genome-wide protein-DNA interactions. To address the lack of powerful ChIP-Seq
+analysis method, we present a novel algorithm, named Model-based Analysis of
+ChIP-Seq (MACS), for identifying transcript factor binding sites. MACS captures
+the influence of genome complexity to evaluate the significance of enriched
+ChIP regions, and MACS improves the spatial resolution of binding sites through
+combining the information of both sequencing tag position and orientation. MACS
+can be easily used for ChIP-Seq data alone, or with control sample with the
+increase of specificity.
+
+View the original MACS2 documentation:
+https://github.com/taoliu/MACS/blob/master/README.rst
+
+------
+
+**Usage**
+
+**Peak Calling**: Main MACS2 Function to Call peaks from alignment results.
+
+**Compare .bdg files**: Deduct noise by comparing two signal tracks in bedGraph.
+
+
+------
+
+**Citation**
+
+For the underlying tool, please cite Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
+
+Integration of MACS2 with Galaxy performed by Ziru Zhou ( ziruzhou@gmail.com ). Please send your comments/questions to modENCODE DCC at help@modencode.org.
+  </help>
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml	Wed Jan 21 11:07:37 2015 -0500
@@ -0,0 +1,34 @@
+<?xml version="1.0"?>
+<tool_dependency>
+  <package name="python" version="2.7">
+    <repository changeset_revision="74d009a2a4cb" name="package_python_2_7" owner="iuc" prior_installation_required="True" toolshed="https://testtoolshed.g2.bx.psu.edu" />
+  </package>
+  <package name="numpy" version="1.8.1">
+    <repository changeset_revision="0f9f634dec8a" name="package_numpy_1_8" owner="iuc" prior_installation_required="True" toolshed="https://testtoolshed.g2.bx.psu.edu" />
+  </package>
+  <package name="macs2" version="2.1.0.20140616">
+      <install version="1.0">
+            <actions>
+              <action type="download_by_url">https://pypi.python.org/packages/source/M/MACS2/MACS2-2.1.0.20140616.tar.gz</action>
+              <action type="set_environment_for_install">
+		<repository changeset_revision="74d009a2a4cb" name="package_python_2_7" owner="iuc" toolshed="https://testtoolshed.g2.bx.psu.edu">
+		  <package name="python" version="2.7" />
+		</repository>
+		<repository changeset_revision="0f9f634dec8a" name="package_numpy_1_8" owner="iuc" toolshed="https://testtoolshed.g2.bx.psu.edu">
+		  <package name="numpy" version="1.8.1" />
+		</repository>
+	      </action>
+              <action type="make_directory">$INSTALL_DIR/lib/python</action>
+              <action type="shell_command">
+                export PYTHONPATH=$PYTHONPATH:$INSTALL_DIR/lib/python &amp;&amp; 
+                python setup.py install --install-lib $INSTALL_DIR/lib/python --install-scripts $INSTALL_DIR/bin
+              </action>
+              <action type="set_environment">
+                <environment_variable action="prepend_to" name="PYTHONPATH">$INSTALL_DIR/lib/python</environment_variable>
+                <environment_variable action="prepend_to" name="PATH">$INSTALL_DIR/bin</environment_variable>
+              </action>
+            </actions>
+      </install>
+      <readme>Macs2.1 depends on having python2.7 and numpy 1.8 installed on all nodes of the work cluster</readme>
+    </package>
+</tool_dependency>