annotate macs2_callpeak.xml @ 4:56e104999978 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 9df9b52baf62b70fbcfc3fbe965d7197d4e8738e
author iuc
date Mon, 06 Feb 2017 02:30:37 -0500
parents 6d4babad010f
children beb902da6e5f
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
1 <tool id="macs2_callpeak" name="MACS2 callpeak" version="@VERSION_STRING@.0">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
2 <description>Call peaks from alignment results</description>
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
3 <macros>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
4 <import>macs2_macros.xml</import>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
5 </macros>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
6 <expand macro="requirements">
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
7 <requirement type="package" version="3.1.2">R</requirement>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
8 <requirement type="set_environment">MACS2_SCRIPT_PATH</requirement>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
9 </expand>
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
10 <expand macro="stdio" />
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
11 <expand macro="version_command" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
12 <command>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
13 #set $temp_stderr = 'macs2_stderr'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
14 (macs2 callpeak
fe62ba547975 Uploaded
iuc
parents:
diff changeset
15
fe62ba547975 Uploaded
iuc
parents:
diff changeset
16 --name "MACS2"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
17 -t ${ ' '.join( map( lambda x:'"%s"' % ( x ), $input_treatment_file ) ) }
fe62ba547975 Uploaded
iuc
parents:
diff changeset
18
fe62ba547975 Uploaded
iuc
parents:
diff changeset
19 #if str( $input_control_file ) != 'None':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
20 -c ${ ' '.join( map( lambda x:'"%s"' % ( x ), $input_control_file ) ) }
fe62ba547975 Uploaded
iuc
parents:
diff changeset
21 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
22
fe62ba547975 Uploaded
iuc
parents:
diff changeset
23 #for $ifile in $input_treatment_file:
fe62ba547975 Uploaded
iuc
parents:
diff changeset
24 #if $ifile.ext.upper() == 'BAM' and $bampe:
fe62ba547975 Uploaded
iuc
parents:
diff changeset
25 --format BAMPE
fe62ba547975 Uploaded
iuc
parents:
diff changeset
26 #else
fe62ba547975 Uploaded
iuc
parents:
diff changeset
27 --format='$ifile.ext.upper()'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
28 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
29 #end for
fe62ba547975 Uploaded
iuc
parents:
diff changeset
30
fe62ba547975 Uploaded
iuc
parents:
diff changeset
31 @effective_genome_size@
fe62ba547975 Uploaded
iuc
parents:
diff changeset
32
fe62ba547975 Uploaded
iuc
parents:
diff changeset
33 --bw='$band_width'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
34
fe62ba547975 Uploaded
iuc
parents:
diff changeset
35 ## advanced options
fe62ba547975 Uploaded
iuc
parents:
diff changeset
36 #if $advanced_options.advanced_options_selector == 'on':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
37 $advanced_options.nolambda
fe62ba547975 Uploaded
iuc
parents:
diff changeset
38 $advanced_options.to_large
fe62ba547975 Uploaded
iuc
parents:
diff changeset
39 --ratio $advanced_options.ratio
fe62ba547975 Uploaded
iuc
parents:
diff changeset
40 --slocal $advanced_options.slocal
fe62ba547975 Uploaded
iuc
parents:
diff changeset
41 --llocal $advanced_options.llocal
fe62ba547975 Uploaded
iuc
parents:
diff changeset
42 #if $advanced_options.broad_options.broad_options_selector == 'broad':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
43 --broad
fe62ba547975 Uploaded
iuc
parents:
diff changeset
44 --broad-cutoff='${ advanced_options.broad_options.broad_cutoff }'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
45 #else:
fe62ba547975 Uploaded
iuc
parents:
diff changeset
46 $advanced_options.broad_options.call_summits
fe62ba547975 Uploaded
iuc
parents:
diff changeset
47 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
48
fe62ba547975 Uploaded
iuc
parents:
diff changeset
49 #if str( $advanced_options.keep_dup_options.keep_dup_options_selector ) == "user":
fe62ba547975 Uploaded
iuc
parents:
diff changeset
50 --keep-dup "${ advanced_options.keep_dup_options.user_keepdup }"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
51 #else
fe62ba547975 Uploaded
iuc
parents:
diff changeset
52 --keep-dup "${ advanced_options.keep_dup_options.keep_dup_options_selector }"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
53 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
54
fe62ba547975 Uploaded
iuc
parents:
diff changeset
55 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
56
fe62ba547975 Uploaded
iuc
parents:
diff changeset
57 ## With --bdg two additional output files will be generated.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
58 #if 'bdg' in str($outputs).split(','):
fe62ba547975 Uploaded
iuc
parents:
diff changeset
59 --bdg
fe62ba547975 Uploaded
iuc
parents:
diff changeset
60 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
61
fe62ba547975 Uploaded
iuc
parents:
diff changeset
62 ## cutoff selection
fe62ba547975 Uploaded
iuc
parents:
diff changeset
63 #if str( $cutoff_options.cutoff_options_selector ) == 'qvalue':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
64 --qvalue "${ cutoff_options.qvalue }"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
65 #elif str( $cutoff_options.cutoff_options_selector ) == 'pvalue':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
66 #if str($cutoff_options.pvalue).strip() != '':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
67 --pvalue "${ cutoff_options.pvalue }"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
68 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
69 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
70
fe62ba547975 Uploaded
iuc
parents:
diff changeset
71 ## model options
fe62ba547975 Uploaded
iuc
parents:
diff changeset
72 #if $nomodel_type.nomodel_type_selector == 'nomodel':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
73 --nomodel
fe62ba547975 Uploaded
iuc
parents:
diff changeset
74 ##--shiftsize '$nomodel_type.shiftsize'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
75 --extsize '${ nomodel_type.extsize }'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
76 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
77
fe62ba547975 Uploaded
iuc
parents:
diff changeset
78 2> $temp_stderr)
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
79 #if 'peaks_tabular' in str($outputs).split(','):
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
80 &amp;&amp;
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
81 cp MACS2_peaks.xls "${ output_tabular }"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
82 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
83
fe62ba547975 Uploaded
iuc
parents:
diff changeset
84 ## run R to create pdf from model script
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
85 #if $nomodel_type.nomodel_type_selector == 'create_model' and 'pdf' in str($outputs).split(','):
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
86 &amp;&amp;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
87 Rscript MACS2_model.r > MACS2_model.r.log
fe62ba547975 Uploaded
iuc
parents:
diff changeset
88 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
89
fe62ba547975 Uploaded
iuc
parents:
diff changeset
90 #if 'html' in str($outputs).split(','):
fe62ba547975 Uploaded
iuc
parents:
diff changeset
91 ## if output files exists, move them to the files_path and create a html result page linking to them
fe62ba547975 Uploaded
iuc
parents:
diff changeset
92 &amp;&amp;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
93 (
fe62ba547975 Uploaded
iuc
parents:
diff changeset
94 count=`ls -1 MACS2* 2>/dev/null | wc -l`;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
95 if [ \$count != 0 ];
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
96 then
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
97 mkdir "${ output_extra_files.files_path }";
fe62ba547975 Uploaded
iuc
parents:
diff changeset
98 cp MACS2* "${ output_extra_files.files_path }";
fe62ba547975 Uploaded
iuc
parents:
diff changeset
99 python "\$MACS2_SCRIPT_PATH/dir2html.py" "${ output_extra_files.files_path }" $temp_stderr > "${ output_extra_files }";
fe62ba547975 Uploaded
iuc
parents:
diff changeset
100 fi;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
101 )
fe62ba547975 Uploaded
iuc
parents:
diff changeset
102 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
103 ;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
104 exit_code_for_galaxy=\$?;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
105 cat $temp_stderr 2&gt;&amp;1;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
106 (exit \$exit_code_for_galaxy)
fe62ba547975 Uploaded
iuc
parents:
diff changeset
107 </command>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
108 <inputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
109 <param name="input_treatment_file" type="data" format="bam,sam,bed" multiple="True" label="ChIP-Seq Treatment File" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
110 <param name="input_control_file" type="data" format="bam,sam,bed" multiple="True" optional="True" label="ChIP-Seq Control File" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
111
fe62ba547975 Uploaded
iuc
parents:
diff changeset
112 <param name="bampe" type="boolean" truevalue="--format BAMPE" falsevalue="" checked="False" label="Are your inputs Paired-end BAM files?"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
113 help="The 'Build model step' will be ignored and the real fragments will be used for each template defined by leftmost and rightmost mapping positions. (--format BAMPE)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
114
fe62ba547975 Uploaded
iuc
parents:
diff changeset
115 <expand macro="conditional_effective_genome_size" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
116 <expand macro="band_width" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
117
fe62ba547975 Uploaded
iuc
parents:
diff changeset
118 <conditional name="cutoff_options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
119 <param name="cutoff_options_selector" type="select" label="Peak detection based on" help="default uses q-value">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
120 <option value="qvalue" selected="true">q-value</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
121 <option value="pvalue">p-value</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
122 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
123 <when value="pvalue">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
124 <param name="pvalue" type="float" value="" label="p-value cutoff for peak detection"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
125 help="default: not set (--pvalue)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
126 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
127 <when value="qvalue">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
128 <param name="qvalue" type="float" value="0.05" label="Minimum FDR (q-value) cutoff for peak detection"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
129 help="default: 0.05 (--qvalue)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
130 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
131 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
132
fe62ba547975 Uploaded
iuc
parents:
diff changeset
133 <conditional name="nomodel_type">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
134 <param name="nomodel_type_selector" type="select" label="Build Model">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
135 <option value="nomodel">Do not build the shifting model (--nomodel)</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
136 <option value="create_model" selected="true">Build the shifting model</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
137 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
138 <when value="create_model"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
139 <when value="nomodel">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
140 <!--<param name="shiftsize" type="integer" label="Arbitrary shift size in bp" value="100" help="(shiftsize)"/>-->
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
141 <param name="extsize" type="integer" value="100" label="The arbitrary extension size in bp"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
142 help="MACS will use this value as fragment size to extend each read towards 3' end, then pile them up. It's exactly twice the number of legacy shiftsize. In previous language, each read is moved 3' direction to middle of fragment by 1/2 d, then extended to both direction with 1/2 d. This is equivalent to say each read is extended towards 3' into a d size fragment. DEFAULT: 200 (--extsize)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
143 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
144 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
145
4
56e104999978 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 9df9b52baf62b70fbcfc3fbe965d7197d4e8738e
iuc
parents: 3
diff changeset
146 <param name="outputs" type="select" display="checkboxes" multiple="True" optional="false" label="Outputs" help="PDF only created when model is build">
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
147 <option value="peaks_tabular" selected="True">Peaks as tabular file</option>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
148 <!--<option value="narrow">narrow Peaks</option>-->
fe62ba547975 Uploaded
iuc
parents:
diff changeset
149 <option value="summits" selected="true">summits</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
150 <option value="bdg" selected="true">Scores in bedGraph files (--bdg)</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
151 <option value="html">Summary page (html)</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
152 <option value="pdf">Plot in PDF</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
153 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
154
fe62ba547975 Uploaded
iuc
parents:
diff changeset
155 <conditional name="advanced_options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
156 <param name="advanced_options_selector" type="select" label="Advanced options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
157 <option value="off" selected="true">Hide advanced options</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
158 <option value="on">Display advanced options</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
159 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
160 <when value="on">
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
161 <param name="to_large" type="boolean" truevalue="--to-large" falsevalue="" checked="False"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
162 label="When set, scale the small sample up to the bigger sample"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
163 help="By default, the bigger dataset will be scaled down towards the smaller dataset, which will lead to smaller p/qvalues and more specific results. Keep in mind that scaling down will bring down background noise more. (--to-large)"/>
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
164 <param name="nolambda" type="boolean" truevalue="--nolambda" falsevalue="" checked="False"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
165 label="Use fixed background lambda as local lambda for every peak region" help="up to 9X more time consuming (--nolambda)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
166 <param name="ratio" type="float" value="1.0"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
167 label="When set, use a custom scaling ratio of ChIP/control (e.g. calculated using NCIS) for linear scaling"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
168 help="(--ratio)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
169 <param name="slocal" value="1000" type="integer" label="The small nearby region in basepairs to calculate dynamic lambda"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
170 help="This is used to capture the bias near the peak summit region. Invalid if there is no control data. If you set this to 0, MACS will skip slocal lambda calculation. *Note* that MACS will always perform a d-size local lambda calculation. The final local bias should be the maximum of the lambda value from d, slocal, and llocal size windows. (--slocal)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
171 <param name="llocal" value="10000" type="integer" label="The large nearby region in basepairs to calculate dynamic lambda"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
172 help="This is used to capture the surround bias. If you set this to 0, MACS will skip llocal lambda calculation. *Note* that MACS will always perform a d-size local lambda calculation. The final local bias should be the maximum of the lambda value from d, slocal, and llocal size windows. (--llocal)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
173 <conditional name="broad_options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
174 <param name="broad_options_selector" type="select"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
175 label="Composite broad regions" help="by putting nearby highly enriched regions into a broad region with loose cutoff (--broad)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
176 <option value="nobroad" selected="true">No broad regions</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
177 <option value="broad">broad regions</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
178 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
179 <when value="broad">
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
180 <param name="broad_cutoff" type="float" label="Cutoff for broad region" value="0.1"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
181 help="value is either p-value or q-value as specified above (--broad-cutoff)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
182 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
183 <when value="nobroad">
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
184 <param name="call_summits" type="boolean" truevalue="--call-summits" falsevalue="" checked="False"
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
185 label="Use a more sophisticated signal processing approach to find subpeak summits in each enriched peak region"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
186 help="(--call-summits)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
187 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
188 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
189 <expand macro="keep_duplicates" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
190 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
191 <when value="off" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
192 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
193 </inputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
194 <outputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
195 <!--callpeaks output-->
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
196 <data name="output_tabular" format="tabular" label="${tool.name} on ${on_string} (Peaks in tabular format)">
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
197 <filter>'peaks_tabular' in outputs</filter>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
198 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
199 <data name="output_broadpeaks" format="bed" from_work_dir="MACS2_peaks.broadPeak" label="${tool.name} on ${on_string} (broad Peaks)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
200 <filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
201 ((
fe62ba547975 Uploaded
iuc
parents:
diff changeset
202 advanced_options['advanced_options_selector'] == "on" and
fe62ba547975 Uploaded
iuc
parents:
diff changeset
203 advanced_options['broad_options']['broad_options_selector'] == "broad"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
204 ))
fe62ba547975 Uploaded
iuc
parents:
diff changeset
205 </filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
206 </data>
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
207 <data name="output_gappedpeaks" format="bed" from_work_dir="MACS2_peaks.gappedPeak" label="${tool.name} on ${on_string} (gapped Peaks)">
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
208 <filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
209 ((
fe62ba547975 Uploaded
iuc
parents:
diff changeset
210 advanced_options['advanced_options_selector'] == "on" and
fe62ba547975 Uploaded
iuc
parents:
diff changeset
211 advanced_options['broad_options']['broad_options_selector'] == "broad"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
212 ))
fe62ba547975 Uploaded
iuc
parents:
diff changeset
213 </filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
214 </data>
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
215 <data name="output_narrowpeaks" format="bed" from_work_dir="MACS2_peaks.narrowPeak" label="${tool.name} on ${on_string} (narrow Peaks)">
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
216 <filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
217 (
fe62ba547975 Uploaded
iuc
parents:
diff changeset
218 advanced_options['advanced_options_selector'] == "off" or
fe62ba547975 Uploaded
iuc
parents:
diff changeset
219 (
fe62ba547975 Uploaded
iuc
parents:
diff changeset
220 advanced_options['advanced_options_selector'] == "on" and
fe62ba547975 Uploaded
iuc
parents:
diff changeset
221 advanced_options['broad_options']['broad_options_selector'] == "nobroad"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
222 ))
fe62ba547975 Uploaded
iuc
parents:
diff changeset
223 </filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
224 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
225 <data name="output_summits" format="bed" from_work_dir="MACS2_summits.bed" label="${tool.name} on ${on_string} (summits in BED)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
226 <filter>'summits' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
227 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
228 <data name="output_plot" format="pdf" from_work_dir="MACS2_model.pdf" label="${tool.name} on ${on_string} (plot)">
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
229 <filter>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
230 ((
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
231 'pdf' in outputs and
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
232 nomodel_type['nomodel_type_selector'] == "create_model"
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
233 ))
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
234 </filter>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
235 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
236 <data name="output_treat_pileup" format="bedgraph" from_work_dir="MACS2_treat_pileup.bdg" label="${tool.name} on ${on_string} (Bedgraph Treatment)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
237 <filter>'bdg' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
238 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
239 <data name="output_control_lambda" format="bedgraph" from_work_dir="MACS2_control_lambda.bdg" label="${tool.name} on ${on_string} (Bedgraph Control)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
240 <filter>'bdg' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
241 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
242 <data name="output_extra_files" format="html" label="${tool.name} on ${on_string} (html report)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
243 <filter>'html' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
244 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
245 </outputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
246 <tests>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
247 <test>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
248 <param name="input_control_file" value="Control_200K.bed" ftype="bed"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
249 <param name="input_treatment_file" value="ChIP_200K.bed" ftype="bed"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
250 <param name="cutoff_options_selector" value="qvalue"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
251 <param name="qvalue" value="0.05"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
252 <param name="band_width" value="300"/>
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
253 <param name="outputs" value="peaks_tabular,bdg"/>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
254 <param name="effective_genome_size_options_selector" value="user_defined" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
255 <param name="gsize" value="3300000000" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
256 <output name="output_control_lambda" compare="contains" file="callpeak_control_part.bdg"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
257 <output name="output_treat_pileup" compare="contains" file="callpeak_treatment_part.bdg"/>
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
258 <output name="output_tabular" compare="contains" file="callpeak_part.tabular"/>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
259 <output name="output_summits" compare="contains" file="callpeak_summits_part.bed"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
260 </test>
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
261 <test>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
262 <param name="input_control_file" value="Control_200K.bed" ftype="bed"/>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
263 <param name="input_treatment_file" value="ChIP_200K.bed" ftype="bed"/>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
264 <param name="cutoff_options_selector" value="qvalue"/>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
265 <param name="qvalue" value="0.05"/>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
266 <param name="band_width" value="300"/>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
267 <param name="outputs" value="pdf"/>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
268 <param name="effective_genome_size_options_selector" value="user_defined" />
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
269 <param name="gsize" value="3300000000" />
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
270 <output name="output_plot" file="magic.pdf" ftype="pdf" compare="contains" />
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
271 </test>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
272 </tests>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
273 <help>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
274 **What it does**
fe62ba547975 Uploaded
iuc
parents:
diff changeset
275
fe62ba547975 Uploaded
iuc
parents:
diff changeset
276 With the improvement of sequencing techniques, chromatin immunoprecipitation followed by high throughput sequencing (ChIP-Seq)
fe62ba547975 Uploaded
iuc
parents:
diff changeset
277 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
fe62ba547975 Uploaded
iuc
parents:
diff changeset
278 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
fe62ba547975 Uploaded
iuc
parents:
diff changeset
279 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.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
280
fe62ba547975 Uploaded
iuc
parents:
diff changeset
281 View the original MACS2 documentation: https://github.com/taoliu/MACS/blob/master/README
fe62ba547975 Uploaded
iuc
parents:
diff changeset
282
fe62ba547975 Uploaded
iuc
parents:
diff changeset
283 ------
fe62ba547975 Uploaded
iuc
parents:
diff changeset
284
fe62ba547975 Uploaded
iuc
parents:
diff changeset
285 **Usage**
fe62ba547975 Uploaded
iuc
parents:
diff changeset
286
fe62ba547975 Uploaded
iuc
parents:
diff changeset
287 **Peak Calling**: Main MACS2 Function to Call peaks from alignment results.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
288
fe62ba547975 Uploaded
iuc
parents:
diff changeset
289 If you choose "Scores in bedGraph files" MACS will output the fragment pileup, control lambda, -log10-pvalue and -log10-qvalue scores in bedGraph files.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
290 The peaks in BED format contain the following colomns: chr end length abs_summit pileup -log10(pvalue) fold_enrichment -log10(qvalue) name
fe62ba547975 Uploaded
iuc
parents:
diff changeset
291
fe62ba547975 Uploaded
iuc
parents:
diff changeset
292 **Compare .bdg files**: Deduct noise by comparing two signal tracks in bedGraph.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
293
fe62ba547975 Uploaded
iuc
parents:
diff changeset
294
fe62ba547975 Uploaded
iuc
parents:
diff changeset
295 @citation@
fe62ba547975 Uploaded
iuc
parents:
diff changeset
296 </help>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
297 <expand macro="citations" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
298 </tool>