|
10
|
1 <macros>
|
|
27
|
2 <xml name="bamCorrelate_mode_actions">
|
|
30
|
3 <param name="region" type="text" value=""
|
|
|
4 label="Region of the genome to limit the operation to"
|
|
|
5 help="This is useful when testing parameters to reduce the computing time. The format is chr:start:end, for example "chr10" or "chr10:456700:891000"" />
|
|
|
6
|
|
10
|
7 <conditional name="advancedOpt">
|
|
|
8 <param name="showAdvancedOpt" type="select" label="Show advanced options" >
|
|
|
9 <option value="no" selected="true">no</option>
|
|
|
10 <option value="yes">yes</option>
|
|
|
11 </param>
|
|
|
12 <when value="no" />
|
|
|
13 <when value="yes">
|
|
|
14 <param name="doNotExtendPairedEnds" type="boolean" truevalue="--doNotExtendPairedEnds" falsevalue=""
|
|
|
15 label="Do not extend paired ends"
|
|
|
16 help="If set, reads are not extended to match the fragment length reported in the BAM file, instead they will be extended to match the fragment length. Default is to extend the reads if paired end information is available."/>
|
|
|
17
|
|
|
18 <param name="ignoreDuplicates" type="boolean" truevalue="--ignoreDuplicates" falsevalue=""
|
|
|
19 label="Ignore duplicates"
|
|
|
20 help="If set, reads that have the same orientation and start position will be considered only once. If reads are paired, the mate position also has to coincide to ignore a read." />
|
|
|
21
|
|
|
22 <param name="minMappingQuality" type="integer" optional="true" value="1" min="1"
|
|
|
23 label="Minimum mapping quality"
|
|
|
24 help= "If set, only reads that have a mapping quality score higher than the given value are considered. *Note* Bowtie's Mapping quality is related to uniqueness: the higher the score, the more unique is a read. A mapping quality defined by Bowtie of 10 or less indicates that there is at least a 1 in 10 chance that the read truly originated elsewhere."/>
|
|
|
25
|
|
|
26 <param name="includeZeros" type="boolean" truevalue="--includeZeros" falsevalue=""
|
|
|
27 label ="Include zeros"
|
|
|
28 help ="If set, then regions with zero counts for *all* BAM files given are included. The default behavior is to ignore those cases." />
|
|
|
29
|
|
17
|
30 <param name="zMin" type="integer" value="" optional="true" label="Minimum value for the heatmap intensities"
|
|
10
|
31 help="If not specified the value is set automatically."/>
|
|
|
32
|
|
17
|
33 <param name="zMax" type="integer" value="" optional="true" label="Maximum value for the heatmap intensities"
|
|
10
|
34 help="If not specified the value is set automatically."/>
|
|
|
35
|
|
|
36 <expand macro="colormap" />
|
|
|
37 </when>
|
|
|
38 </conditional>
|
|
27
|
39 </xml>
|
|
10
|
40 <token name="@THREADS@">--numberOfProcessors "\${GALAXY_SLOTS:-4}"</token>
|
|
27
|
41 <xml name="requirements">
|
|
10
|
42 <requirements>
|
|
34
|
43 <requirement type="binary">@BINARY@</requirement>
|
|
11
|
44 <requirement type="package" >samtools</requirement>
|
|
10
|
45 <requirement type="package" >deepTools</requirement>
|
|
|
46 <requirement type="package" >ucsc_tools</requirement>
|
|
34
|
47 <requirement type="package" version="1.5.4_570648e1d7a44583a72852f2bd2bf408f76f8f6a">deepTools</requirement>
|
|
15
|
48 <requirement type="package" version="0.1">ucsc_tools</requirement>
|
|
|
49 <requirement type="package" version="1.7.1">numpy</requirement>
|
|
|
50 <requirement type="package" version="0.7.7">pysam</requirement>
|
|
22
|
51 <requirement type="package" version="0.12.0">scipy</requirement>
|
|
20
|
52 <requirement type="package" version="1.2.1">matplotlib</requirement>
|
|
15
|
53 <requirement type="package" version="0.1.19">samtools</requirement>
|
|
27
|
54 <requirement type="package" version="12-2013">bx-python</requirement>
|
|
10
|
55 <yield />
|
|
|
56 </requirements>
|
|
34
|
57 <version_command>@BINARY@ --version</version_command>
|
|
27
|
58 </xml>
|
|
|
59 <xml name="stdio">
|
|
22
|
60 <stdio>
|
|
|
61 <exit_code range="1:" />
|
|
|
62 <exit_code range=":-1" />
|
|
|
63 <regex match="Error:" />
|
|
|
64 <regex match="Exception:" />
|
|
|
65 <regex match="EXception:" />
|
|
25
|
66 <regex match="Traceback" />
|
|
22
|
67 </stdio>
|
|
27
|
68 </xml>
|
|
10
|
69 <token name="@REFERENCES@">
|
|
|
70
|
|
|
71 For more informations, please visit our `project site`_.
|
|
|
72
|
|
|
73 If you would like to give us feedback or you run into any trouble, please send an email to deeptools@googlegroups.com
|
|
|
74
|
|
|
75 This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
|
|
|
76
|
|
|
77
|
|
|
78 .. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
|
|
|
79 .. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
|
|
|
80 .. _project site: https://github.com/fidelram/deepTools
|
|
|
81
|
|
|
82 </token>
|
|
|
83
|
|
27
|
84 <xml name="multiple_input_bams">
|
|
21
|
85 <repeat name="input_files" title="BAM files" min="2">
|
|
10
|
86 <param name="bamfile" type="data" format="bam"
|
|
|
87 label="Bam file"
|
|
|
88 help="The BAM file must be sorted."/>
|
|
|
89 <param name="label" type="text" size="30" optional="true" value=""
|
|
|
90 label="Label"
|
|
|
91 help="Label to use in the output. If not given the dataset name will be used instead."/>
|
|
|
92 </repeat>
|
|
27
|
93 </xml>
|
|
10
|
94
|
|
14
|
95 <token name="@multiple_input_bams@">
|
|
29
|
96 #import tempfile
|
|
|
97 #set $temp_dir = os.path.abspath(tempfile.mkdtemp())
|
|
|
98 #set files=[]
|
|
|
99 #set labels=[]
|
|
14
|
100 #for $i in $input_files:
|
|
|
101 #set $temp_input_handle = tempfile.NamedTemporaryFile( dir=$temp_dir )
|
|
|
102 #set $temp_input_path = $temp_input_handle.name
|
|
|
103 #silent $temp_input_handle.close()
|
|
|
104 #silent os.system("ln -s %s %s.bam" % (str($i.bamfile), $temp_input_path))
|
|
|
105 #silent os.system("ln -s %s %s.bam.bai" % (str($i.bamfile.metadata.bam_index), $temp_input_path))
|
|
|
106 #silent $files.append('%s.bam' % $temp_input_path)
|
|
|
107
|
|
|
108 ##set $files += [str($i.bamfile)]
|
|
|
109 #if str($i.label.value) != "":
|
|
|
110 #set $labels += ["\"%s\"" % ($i.label.value)]
|
|
|
111 #else
|
|
|
112 #set $labels += ["\"%s\"" % ($i.bamfile.name)]
|
|
|
113 #end if
|
|
|
114 #end for
|
|
|
115 </token>
|
|
|
116
|
|
27
|
117 <xml name="reference_genome_source">
|
|
10
|
118 <conditional name="source">
|
|
|
119 <param name="ref_source" type="select" label="Reference genome">
|
|
|
120 <option value="cached">locally cached</option>
|
|
|
121 <option value="history">in your history</option>
|
|
|
122 </param>
|
|
|
123 <when value="cached">
|
|
|
124 <param name="input1_2bit" type="select" label="Using reference genome" help="If your genome of interest is not listed, contact the Galaxy team">
|
|
17
|
125 <options from_data_table="deepTools_seqs">
|
|
|
126 <filter type="sort_by" column="1" />
|
|
|
127 <validator type="no_options" message="No indexes are available." />
|
|
|
128 </options>
|
|
10
|
129 </param>
|
|
|
130 </when>
|
|
|
131 <when value="history">
|
|
|
132 <param name="input1" type="data" format="twobit" label="Select a reference dataset in 2bit format" />
|
|
|
133 </when>
|
|
|
134 </conditional>
|
|
27
|
135 </xml>
|
|
10
|
136
|
|
|
137 <token name="@reference_genome_source@">
|
|
|
138 #if $source.ref_source=="history":
|
|
|
139 --genome $source.input1
|
|
|
140 #else:
|
|
|
141 --genome "${source.input1_2bit.fields.path}"
|
|
|
142 #end if
|
|
|
143 </token>
|
|
|
144
|
|
27
|
145 <xml name="effectiveGenomeSize">
|
|
10
|
146 <conditional name="effectiveGenomeSize">
|
|
|
147 <param name="effectiveGenomeSize_opt" type="select" label="Effective genome size"
|
|
|
148 help="The effective genome size is the portion of the genome that is mappable. Large fractions of the genome are stretches of NNNN that should be discarded.
|
|
|
149 Also, if repetitive regions were not included in the mapping of reads, the effective genome size needs to be adjusted accordingly.
|
|
|
150 See Table 2 of http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0030377 or http://www.nature.com/nbt/journal/v27/n1/fig_tab/nbt.1518_T1.html for several effective genome sizes.">
|
|
|
151 <option value="2150570000">mm9</option>
|
|
|
152 <option value="2451960000">hg19</option>
|
|
|
153 <option value="121400000">dm3</option>
|
|
|
154 <option value="93260000">ce10</option>
|
|
|
155 <option value="specific">user specified</option>
|
|
|
156 </param>
|
|
|
157 <when value="specific">
|
|
|
158 <param name="effectiveGenomeSize" type="integer" value="" label="Effective genome size" />
|
|
|
159 </when>
|
|
|
160 </conditional>
|
|
27
|
161 </xml>
|
|
10
|
162
|
|
27
|
163 <xml name="image_file_format">
|
|
22
|
164 <param name="outFileFormat" type="select" label="Image file format">
|
|
|
165 <option value="png" selected="true">png</option>
|
|
|
166 <option value="pdf">pdf</option>
|
|
|
167 <option value="svg">svg</option>
|
|
|
168 <option value="eps">eps</option>
|
|
|
169 <option value="emf">emf</option>
|
|
|
170 </param>
|
|
27
|
171 </xml>
|
|
22
|
172
|
|
27
|
173 <xml name="input_save_matrix_values">
|
|
|
174 <param name="saveMatrix" type="boolean" label="Save the matrix of values underlying the heatmap"/>
|
|
|
175 </xml>
|
|
|
176
|
|
|
177 <xml name="input_graphic_output_settings">
|
|
10
|
178 <conditional name="output" >
|
|
|
179 <param name="showOutputSettings" type="select" label="Show advanced output settings" >
|
|
|
180 <option value="no" selected="true">no</option>
|
|
|
181 <option value="yes">yes</option>
|
|
|
182 </param>
|
|
|
183 <when value="no" />
|
|
|
184 <when value="yes">
|
|
25
|
185 <yield />
|
|
27
|
186 <param name="saveData" type="boolean" label="Save the data underlying the average profile"/>
|
|
10
|
187 <param name="saveSortedRegions" type="boolean" label="Save the regions after skipping zeros or min/max threshold values" help="The order of the regions in the file follows the sorting order selected. This is useful, for example, to generate other heatmaps keeping the sorting of the first heatmap."/>
|
|
|
188 </when>
|
|
|
189 </conditional>
|
|
27
|
190 </xml>
|
|
10
|
191
|
|
27
|
192 <xml name="input_image_file_format">
|
|
25
|
193 <param name="outFileFormat" type="select" label="Image file format">
|
|
|
194 <option value="png" selected="true">png</option>
|
|
|
195 <option value="pdf">pdf</option>
|
|
|
196 <option value="svg">svg</option>
|
|
|
197 <option value="eps">eps</option>
|
|
|
198 <option value="emf">emf</option>
|
|
|
199 </param>
|
|
27
|
200 </xml>
|
|
25
|
201
|
|
27
|
202 <xml name="output_image_file_format">
|
|
24
|
203 <data format="png" name="outFileName" label="${tool.name} image">
|
|
|
204 <change_format>
|
|
|
205 <when input="output.outFileFormat" value="pdf" format="pdf" />
|
|
|
206 <when input="output.outFileFormat" value="svg" format="svg" />
|
|
|
207 <when input="output.outFileFormat" value="eps" format="eps" />
|
|
|
208 <when input="output.outFileFormat" value="emf" format="emf" />
|
|
|
209 </change_format>
|
|
|
210 </data>
|
|
27
|
211 </xml>
|
|
25
|
212
|
|
27
|
213 <xml name="output_save_matrix_values">
|
|
|
214 <data format="tabular" name="outFileNameMatrix" label="${tool.name} on ${on_string}: Heatmap values">
|
|
|
215 <filter>
|
|
|
216 ((
|
|
|
217 output['showOutputSettings'] == 'yes' and
|
|
|
218 output['saveMatrix'] is True
|
|
|
219 ))
|
|
|
220 </filter>
|
|
|
221 </data>
|
|
|
222 </xml>
|
|
|
223
|
|
|
224 <xml name="output_graphic_outputs">
|
|
22
|
225 <data format="tabular" name="outFileNameData" label="${tool.name} on ${on_string}: averages per matrix column">
|
|
|
226 <filter>
|
|
|
227 ((
|
|
|
228 output['showOutputSettings'] == 'yes' and
|
|
|
229 output['saveData'] is True
|
|
|
230 ))
|
|
|
231 </filter>
|
|
|
232 </data>
|
|
|
233 <data format="bed" name="outFileSortedRegions" label="${tool.name} on ${on_string}: sorted/filtered regions">
|
|
|
234 <filter>
|
|
|
235 ((
|
|
|
236 output['showOutputSettings'] == 'yes' and
|
|
|
237 output['saveSortedRegions'] is True
|
|
|
238 ))
|
|
|
239 </filter>
|
|
|
240 </data>
|
|
27
|
241 </xml>
|
|
22
|
242
|
|
27
|
243 <xml name="colormap">
|
|
10
|
244 <param name="colorMap" type="select" label="Color map to use for the heatmap" help=" Available color map names can be found here: http://www.astro.lsa.umich.edu/~msshin/science/code/matplotlib_cm/">
|
|
|
245 <option value="RdYlBu" selected="true">RdYlBu</option>
|
|
|
246 <option value="Accent">Accent</option>
|
|
|
247 <option value="Spectral">Spectral</option>
|
|
|
248 <option value="Set1">Set1</option>
|
|
|
249 <option value="Set2">Set2</option>
|
|
|
250 <option value="Set3">Set3</option>
|
|
|
251 <option value="Dark2">Dark2</option>
|
|
|
252 <option value="Reds">Reds</option>
|
|
|
253 <option value="Oranges">Oranges</option>
|
|
|
254 <option value="Greens">Greens</option>
|
|
|
255 <option value="Blues">Blues</option>
|
|
|
256 <option value="Greys">Greys</option>
|
|
|
257 <option value="Purples">Purples</option>
|
|
|
258 <option value="Paired">Paired</option>
|
|
|
259 <option value="Pastel1">Pastel1</option>
|
|
|
260 <option value="Pastel2">Pastel2</option>
|
|
|
261 <option value="spring">spring</option>
|
|
|
262 <option value="summer">summer</option>
|
|
|
263 <option value="autumn">autumn</option>
|
|
|
264 <option value="winter">winter</option>
|
|
|
265 <option value="hot">hot</option>
|
|
|
266 <option value="coolwarm">coolwarm</option>
|
|
|
267 <option value="cool">cool</option>
|
|
|
268 <option value="seismic">seismic</option>
|
|
|
269 <option value="terrain">terrain</option>
|
|
|
270 <option value="ocean">ocean</option>
|
|
|
271 <option value="rainbow">rainbow</option>
|
|
|
272 <option value="bone">bone</option>
|
|
|
273 <option value="flag">flag</option>
|
|
|
274 <option value="prism">prism</option>
|
|
|
275 <option value="cubehelix">cubehelix</option>
|
|
|
276 <option value="binary">binary</option>
|
|
|
277 <option value="pink">pink</option>
|
|
|
278 <option value="gray">gray</option>
|
|
|
279 <option value="copper">copper</option>
|
|
|
280 <option value="BrBG">BrBG</option>
|
|
|
281 <option value="BuGn">BuGn</option>
|
|
|
282 <option value="BuPu">BuPu</option>
|
|
|
283 <option value="GnBu">GnBu</option>
|
|
|
284 <option value="OrRd">OrRd</option>
|
|
|
285 <option value="PiYG">PiYG</option>
|
|
|
286 <option value="PRGn">PRGn</option>
|
|
|
287 <option value="PuOr">PuOr</option>
|
|
|
288 <option value="PuRd">PuRd</option>
|
|
|
289 <option value="PuBu">PuBu</option>
|
|
|
290 <option value="RdBu">RdBu</option>
|
|
|
291 <option value="RdGy">RdGy</option>
|
|
|
292 <option value="RdPu">RdPu</option>
|
|
|
293 <option value="YlGn">YlGn</option>
|
|
|
294 <option value="PuBuGn">PuBuGn</option>
|
|
|
295 <option value="RdYlGn">RdYlGn</option>
|
|
|
296 <option value="YlGnBu">YlGnBu</option>
|
|
|
297 <option value="YlOrBr">YlOrBr</option>
|
|
|
298 <option value="YlOrRd">YlOrRd</option>
|
|
|
299 <option value="gist_gray">gist_gray</option>
|
|
|
300 <option value="gist_stern">gist_stern</option>
|
|
|
301 <option value="gist_earth">gist_earth</option>
|
|
|
302 <option value="gist_yarg">gist_yarg</option>
|
|
|
303 <option value="gist_ncar">gist_ncar</option>
|
|
|
304 <option value="gist_rainbow">gist_rainbow</option>
|
|
|
305 <option value="gist_heat">gist_heat</option>
|
|
|
306 <option value="gnuplot">gnuplot</option>
|
|
|
307 <option value="gnuplot2">gnuplot2</option>
|
|
|
308 <option value="CMRmap">CMRmap</option>
|
|
|
309 <option value="bwr">bwr</option>
|
|
|
310 <option value="hsv">hsv</option>
|
|
|
311 <option value="brg">brg</option>
|
|
|
312 <option value="jet">jet</option>
|
|
|
313 <option value="afmhot">afmhot</option>
|
|
|
314 <option value="Accent_r">Accent reversed</option>
|
|
|
315 <option value="Spectral_r">Spectral reversed</option>
|
|
|
316 <option value="Set1_r">Set1 reversed</option>
|
|
|
317 <option value="Set2_r">Set2 reversed</option>
|
|
|
318 <option value="Set3_r">Set3 reversed</option>
|
|
|
319 <option value="Dark2_r">Dark2 reversed</option>
|
|
|
320 <option value="Reds_r">Reds reversed</option>
|
|
|
321 <option value="Oranges_r">Oranges reversed</option>
|
|
|
322 <option value="Greens_r">Greens reversed</option>
|
|
|
323 <option value="Blues_r">Blues reversed</option>
|
|
|
324 <option value="Greys_r">Greys reversed</option>
|
|
|
325 <option value="Purples_r">Purples reversed</option>
|
|
|
326 <option value="Paired_r">Paired reversed</option>
|
|
|
327 <option value="Pastel1_r">Pastel1 reversed</option>
|
|
|
328 <option value="Pastel2_r">Pastel2 reversed</option>
|
|
|
329 <option value="spring_r">spring reversed</option>
|
|
|
330 <option value="summer_r">summer reversed</option>
|
|
|
331 <option value="autumn_r">autumn reversed</option>
|
|
|
332 <option value="winter_r">winter reversed</option>
|
|
|
333 <option value="hot_r">hot reversed</option>
|
|
|
334 <option value="coolwarm_r">coolwarm reversed</option>
|
|
|
335 <option value="cool_r">cool reversed</option>
|
|
|
336 <option value="seismic_r">seismic reversed</option>
|
|
|
337 <option value="terrain_r">terrain reversed</option>
|
|
|
338 <option value="ocean_r">ocean reversed</option>
|
|
|
339 <option value="rainbow_r">rainbow reversed</option>
|
|
|
340 <option value="bone_r">bone reversed</option>
|
|
|
341 <option value="flag_r">flag reversed</option>
|
|
|
342 <option value="prism_r">prism reversed</option>
|
|
|
343 <option value="cubehelix_r">cubehelix reversed</option>
|
|
|
344 <option value="binary_r">binary reversed</option>
|
|
|
345 <option value="pink_r">pink reversed</option>
|
|
|
346 <option value="gray_r">gray reversed</option>
|
|
|
347 <option value="copper_r">copper reversed</option>
|
|
|
348 <option value="BrBG_r">BrBG reversed</option>
|
|
|
349 <option value="BuGn_r">BuGn reversed</option>
|
|
|
350 <option value="BuPu_r">BuPu reversed</option>
|
|
|
351 <option value="GnBu_r">GnBu reversed</option>
|
|
|
352 <option value="OrRd_r">OrRd reversed</option>
|
|
|
353 <option value="PiYG_r">PiYG reversed</option>
|
|
|
354 <option value="PRGn_r">PRGn reversed</option>
|
|
|
355 <option value="PuOr_r">PuOr reversed</option>
|
|
|
356 <option value="PuRd_r">PuRd reversed</option>
|
|
|
357 <option value="PuBu_r">PuBu reversed</option>
|
|
|
358 <option value="RdBu_r">RdBu reversed</option>
|
|
|
359 <option value="RdGy_r">RdGy reversed</option>
|
|
|
360 <option value="RdPu_r">RdPu reversed</option>
|
|
|
361 <option value="YlGn_r">YlGn reversed</option>
|
|
|
362 <option value="PuBuGn_r">PuBuGn reversed</option>
|
|
|
363 <option value="RdYlBu_r">RdYlBu reversed</option>
|
|
|
364 <option value="RdYlGn_r">RdYlGn reversed</option>
|
|
|
365 <option value="YlGnBu_r">YlGnBu reversed</option>
|
|
|
366 <option value="YlOrBr_r">YlOrBr reversed</option>
|
|
|
367 <option value="YlOrRd_r">YlOrRd reversed</option>
|
|
|
368 <option value="gist_gray_r">gist_gray reversed</option>
|
|
|
369 <option value="gist_stern_r">gist_stern reversed</option>
|
|
|
370 <option value="gist_earth_r">gist_earth reversed</option>
|
|
|
371 <option value="gist_yarg_r">gist_yarg reversed</option>
|
|
|
372 <option value="gist_ncar_r">gist_ncar reversed</option>
|
|
|
373 <option value="gist_rainbow_r">gist_rainbow reversed</option>
|
|
|
374 <option value="gist_heat_r">gist_heat reversed</option>
|
|
|
375 <option value="gnuplot_r">gnuplot reversed</option>
|
|
|
376 <option value="gnuplot2_r">gnuplot2 reversed</option>
|
|
|
377 <option value="CMRmap_r">CMRmap reversed</option>
|
|
|
378 <option value="bwr_r">bwr reversed</option>
|
|
|
379 <option value="hsv_r">hsv reversed</option>
|
|
|
380 <option value="brg_r">brg reversed</option>
|
|
|
381 <option value="jet_r">jet reversed</option>
|
|
|
382 <option value="afmhot_r">afmhot reversed</option>
|
|
|
383 </param>
|
|
|
384
|
|
27
|
385 </xml>
|
|
10
|
386
|
|
|
387 </macros>
|
