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