annotate macs2_callpeak.xml @ 3:6d4babad010f draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
author iuc
date Wed, 27 Jan 2016 15:23:25 -0500
parents bfe57d6e0c4c
children 56e104999978
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 #else:
fe62ba547975 Uploaded
iuc
parents:
diff changeset
70 --foldenrichment "${ cutoff_options.foldenrichment }"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
71 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
72
fe62ba547975 Uploaded
iuc
parents:
diff changeset
73 ## model options
fe62ba547975 Uploaded
iuc
parents:
diff changeset
74 #if $nomodel_type.nomodel_type_selector == 'nomodel':
fe62ba547975 Uploaded
iuc
parents:
diff changeset
75 --nomodel
fe62ba547975 Uploaded
iuc
parents:
diff changeset
76 ##--shiftsize '$nomodel_type.shiftsize'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
77 --extsize '${ nomodel_type.extsize }'
fe62ba547975 Uploaded
iuc
parents:
diff changeset
78 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
79
fe62ba547975 Uploaded
iuc
parents:
diff changeset
80 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
81 #if 'peaks_tabular' in str($outputs).split(','):
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
82 &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
83 cp MACS2_peaks.xls "${ output_tabular }"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
84 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
85
fe62ba547975 Uploaded
iuc
parents:
diff changeset
86 ## 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
87 #if $nomodel_type.nomodel_type_selector == 'create_model' and 'pdf' in str($outputs).split(','):
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
88 &amp;&amp;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
89 Rscript MACS2_model.r > MACS2_model.r.log
fe62ba547975 Uploaded
iuc
parents:
diff changeset
90 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
91
fe62ba547975 Uploaded
iuc
parents:
diff changeset
92 #if 'html' in str($outputs).split(','):
fe62ba547975 Uploaded
iuc
parents:
diff changeset
93 ## 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
94 &amp;&amp;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
95 (
fe62ba547975 Uploaded
iuc
parents:
diff changeset
96 count=`ls -1 MACS2* 2>/dev/null | wc -l`;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
97 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
98 then
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
99 mkdir "${ output_extra_files.files_path }";
fe62ba547975 Uploaded
iuc
parents:
diff changeset
100 cp MACS2* "${ output_extra_files.files_path }";
fe62ba547975 Uploaded
iuc
parents:
diff changeset
101 python "\$MACS2_SCRIPT_PATH/dir2html.py" "${ output_extra_files.files_path }" $temp_stderr > "${ output_extra_files }";
fe62ba547975 Uploaded
iuc
parents:
diff changeset
102 fi;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
103 )
fe62ba547975 Uploaded
iuc
parents:
diff changeset
104 #end if
fe62ba547975 Uploaded
iuc
parents:
diff changeset
105 ;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
106 exit_code_for_galaxy=\$?;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
107 cat $temp_stderr 2&gt;&amp;1;
fe62ba547975 Uploaded
iuc
parents:
diff changeset
108 (exit \$exit_code_for_galaxy)
fe62ba547975 Uploaded
iuc
parents:
diff changeset
109 </command>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
110 <inputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
111 <param name="input_treatment_file" type="data" format="bam,sam,bed" multiple="True" label="ChIP-Seq Treatment File" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
112 <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
113
fe62ba547975 Uploaded
iuc
parents:
diff changeset
114 <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
115 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
116
fe62ba547975 Uploaded
iuc
parents:
diff changeset
117 <expand macro="conditional_effective_genome_size" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
118 <expand macro="band_width" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
119
fe62ba547975 Uploaded
iuc
parents:
diff changeset
120 <conditional name="cutoff_options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
121 <param name="cutoff_options_selector" type="select" label="Peak detection based on" help="default uses q-value">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
122 <option value="qvalue" selected="true">q-value</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
123 <option value="pvalue">p-value</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
124 <option value="foldenrichment">foldenrichment</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
125 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
126 <when value="pvalue">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
127 <param name="pvalue" type="float" value="" label="p-value cutoff for peak detection"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
128 help="default: not set (--pvalue)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
129 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
130 <when value="qvalue">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
131 <param name="qvalue" type="float" value="0.05" label="Minimum FDR (q-value) cutoff for peak detection"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
132 help="default: 0.05 (--qvalue)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
133 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
134 <when value="foldenrichment">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
135 <param name="foldenrichment" value="" type="integer" label="Foldenrichment cutoff for peak detection"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
136 help="(--foldenrichment)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
137 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
138 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
139
fe62ba547975 Uploaded
iuc
parents:
diff changeset
140 <conditional name="nomodel_type">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
141 <param name="nomodel_type_selector" type="select" label="Build Model">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
142 <option value="nomodel">Do not build the shifting model (--nomodel)</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
143 <option value="create_model" selected="true">Build the shifting model</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
144 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
145 <when value="create_model"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
146 <when value="nomodel">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
147 <!--<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
148 <param name="extsize" type="integer" value="100" label="The arbitrary extension size in bp"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
149 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
150 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
151 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
152
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
153 <param name="outputs" type="select" display="checkboxes" multiple="True" 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
154 <option value="peaks_tabular" selected="True">Peaks as tabular file</option>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
155 <!--<option value="narrow">narrow Peaks</option>-->
fe62ba547975 Uploaded
iuc
parents:
diff changeset
156 <option value="summits" selected="true">summits</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
157 <option value="bdg" selected="true">Scores in bedGraph files (--bdg)</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
158 <option value="html">Summary page (html)</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
159 <option value="pdf">Plot in PDF</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
160 <validator type="no_options" message="Please select at least one output file." />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
161 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
162
fe62ba547975 Uploaded
iuc
parents:
diff changeset
163 <conditional name="advanced_options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
164 <param name="advanced_options_selector" type="select" label="Advanced options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
165 <option value="off" selected="true">Hide advanced options</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
166 <option value="on">Display advanced options</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
167 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
168 <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
169 <param name="to_large" type="boolean" truevalue="--to-large" falsevalue="" checked="False"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
170 label="When set, scale the small sample up to the bigger sample"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
171 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
172 <param name="nolambda" type="boolean" truevalue="--nolambda" falsevalue="" checked="False"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
173 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
174 <param name="ratio" type="float" value="1.0"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
175 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
176 help="(--ratio)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
177 <param name="slocal" value="1000" type="integer" label="The small nearby region in basepairs to calculate dynamic lambda"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
178 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
179 <param name="llocal" value="10000" type="integer" label="The large nearby region in basepairs to calculate dynamic lambda"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
180 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
181 <conditional name="broad_options">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
182 <param name="broad_options_selector" type="select"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
183 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
184 <option value="nobroad" selected="true">No broad regions</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
185 <option value="broad">broad regions</option>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
186 </param>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
187 <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
188 <param name="broad_cutoff" type="float" label="Cutoff for broad region" value="0.1"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
189 help="value is either p-value or q-value as specified above (--broad-cutoff)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
190 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
191 <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
192 <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
193 label="Use a more sophisticated signal processing approach to find subpeak summits in each enriched peak region"
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
194 help="(--call-summits)"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
195 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
196 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
197 <expand macro="keep_duplicates" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
198 </when>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
199 <when value="off" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
200 </conditional>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
201 </inputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
202 <outputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
203 <!--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
204 <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
205 <filter>'peaks_tabular' in outputs</filter>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
206 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
207 <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
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_gappedpeaks" format="bed" from_work_dir="MACS2_peaks.gappedPeak" label="${tool.name} on ${on_string} (gapped 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'] == "on" and
fe62ba547975 Uploaded
iuc
parents:
diff changeset
219 advanced_options['broad_options']['broad_options_selector'] == "broad"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
220 ))
fe62ba547975 Uploaded
iuc
parents:
diff changeset
221 </filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
222 </data>
3
6d4babad010f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
223 <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
224 <filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
225 (
fe62ba547975 Uploaded
iuc
parents:
diff changeset
226 advanced_options['advanced_options_selector'] == "off" or
fe62ba547975 Uploaded
iuc
parents:
diff changeset
227 (
fe62ba547975 Uploaded
iuc
parents:
diff changeset
228 advanced_options['advanced_options_selector'] == "on" and
fe62ba547975 Uploaded
iuc
parents:
diff changeset
229 advanced_options['broad_options']['broad_options_selector'] == "nobroad"
fe62ba547975 Uploaded
iuc
parents:
diff changeset
230 ))
fe62ba547975 Uploaded
iuc
parents:
diff changeset
231 </filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
232 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
233 <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
234 <filter>'summits' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
235 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
236 <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
237 <filter>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
238 ((
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
239 '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
240 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
241 ))
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
242 </filter>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
243 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
244 <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
245 <filter>'bdg' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
246 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
247 <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
248 <filter>'bdg' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
249 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
250 <data name="output_extra_files" format="html" label="${tool.name} on ${on_string} (html report)">
fe62ba547975 Uploaded
iuc
parents:
diff changeset
251 <filter>'html' in outputs</filter>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
252 </data>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
253 </outputs>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
254 <tests>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
255 <test>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
256 <param name="input_control_file" value="Control_200K.bed" ftype="bed"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
257 <param name="input_treatment_file" value="ChIP_200K.bed" ftype="bed"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
258 <param name="cutoff_options_selector" value="qvalue"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
259 <param name="qvalue" value="0.05"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
260 <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
261 <param name="outputs" value="peaks_tabular,bdg"/>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
262 <param name="effective_genome_size_options_selector" value="user_defined" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
263 <param name="gsize" value="3300000000" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
264 <output name="output_control_lambda" compare="contains" file="callpeak_control_part.bdg"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
265 <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
266 <output name="output_tabular" compare="contains" file="callpeak_part.tabular"/>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
267 <output name="output_summits" compare="contains" file="callpeak_summits_part.bed"/>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
268 </test>
2
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
269 <test>
bfe57d6e0c4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/macs2 commit 63955994af5090ec444c03c221df0012d5ae4e74
iuc
parents: 0
diff changeset
270 <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
271 <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
272 <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
273 <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
274 <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
275 <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
276 <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
277 <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
278 <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
279 </test>
0
fe62ba547975 Uploaded
iuc
parents:
diff changeset
280 </tests>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
281 <help>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
282 **What it does**
fe62ba547975 Uploaded
iuc
parents:
diff changeset
283
fe62ba547975 Uploaded
iuc
parents:
diff changeset
284 With the improvement of sequencing techniques, chromatin immunoprecipitation followed by high throughput sequencing (ChIP-Seq)
fe62ba547975 Uploaded
iuc
parents:
diff changeset
285 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
286 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
287 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
288
fe62ba547975 Uploaded
iuc
parents:
diff changeset
289 View the original MACS2 documentation: https://github.com/taoliu/MACS/blob/master/README
fe62ba547975 Uploaded
iuc
parents:
diff changeset
290
fe62ba547975 Uploaded
iuc
parents:
diff changeset
291 ------
fe62ba547975 Uploaded
iuc
parents:
diff changeset
292
fe62ba547975 Uploaded
iuc
parents:
diff changeset
293 **Usage**
fe62ba547975 Uploaded
iuc
parents:
diff changeset
294
fe62ba547975 Uploaded
iuc
parents:
diff changeset
295 **Peak Calling**: Main MACS2 Function to Call peaks from alignment results.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
296
fe62ba547975 Uploaded
iuc
parents:
diff changeset
297 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
298 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
299
fe62ba547975 Uploaded
iuc
parents:
diff changeset
300 **Compare .bdg files**: Deduct noise by comparing two signal tracks in bedGraph.
fe62ba547975 Uploaded
iuc
parents:
diff changeset
301
fe62ba547975 Uploaded
iuc
parents:
diff changeset
302
fe62ba547975 Uploaded
iuc
parents:
diff changeset
303 @citation@
fe62ba547975 Uploaded
iuc
parents:
diff changeset
304 </help>
fe62ba547975 Uploaded
iuc
parents:
diff changeset
305 <expand macro="citations" />
fe62ba547975 Uploaded
iuc
parents:
diff changeset
306 </tool>