Mercurial > repos > pjbriggs > macs21
annotate macs21_wrapper.xml @ 3:4124781932db draft
Updated to MACS 2.1.1 and use conda for dependency resolution.
| author | pjbriggs |
|---|---|
| date | Tue, 20 Mar 2018 11:25:04 -0400 |
| parents | 00d73c812399 |
| children | 11cf21ee4242 |
| rev | line source |
|---|---|
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
1 <tool id="macs2_1_peakcalling" name="MACS2.1.1" version="2.1.1.0"> |
| 0 | 2 <description>Model-based Analysis of ChIP-Seq: peak calling</description> |
| 3 <requirements> | |
| 4 <requirement type="package" version="2.7">python</requirement> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
5 <requirement type="package" version="2.1.1">macs2</requirement> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
6 <requirement type="package" version="3.4">R</requirement> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
7 <requirement type="package" version="357">ucsc-fetchchromsizes</requirement> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
8 <requirement type="package" version="357">ucsc-bedclip</requirement> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
9 <requirement type="package" version="357">ucsc-bedsort</requirement> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
10 <requirement type="package" version="357">ucsc-bedgraphtobigwig</requirement> |
| 0 | 11 </requirements> |
| 12 <version_command>macs2 --version</version_command> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
13 <command><![CDATA[ |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
14 python $__tool_directory__/macs21_wrapper.py callpeak |
| 0 | 15 ## |
| 16 ## ChIP-seq input | |
| 17 $input_chipseq_file1 | |
| 18 ## | |
| 19 ## ChIP-seq control | |
| 20 #if str($input_control_file1) != 'None' | |
| 21 -c $input_control_file1 | |
| 22 #end if | |
| 23 ## | |
| 1 | 24 --format=$format |
| 0 | 25 --name="$experiment_name" |
| 26 --bw=$bw | |
| 27 ## | |
| 28 ## Genome size | |
| 29 #if str($genome_size.gsize) == '' | |
| 30 --gsize=$genome_size.user_defined_gsize | |
| 31 #else: | |
| 32 --gsize=$genome_size.gsize | |
| 33 #end if | |
| 34 ## | |
| 35 ## Broad peaks | |
| 36 #if str($broad_options.broad_regions) == 'broad' | |
| 37 --broad --broad-cutoff=$broad_options.broad_cutoff | |
| 38 #end if | |
| 39 ## | |
| 40 ## (no)model options | |
| 41 #if str($nomodel_type.nomodel_type_selector) == 'nomodel' | |
| 42 --nomodel --extsize=$nomodel_type.extsize | |
| 43 #end if | |
| 44 ## | |
| 45 ## pq value select options | |
| 46 #if str($pq_options.pq_options_selector) == 'qvalue' | |
| 47 --qvalue=$pq_options.qvalue | |
| 48 #else | |
| 49 --pvalue=$pq_options.pvalue | |
| 50 #end if | |
| 51 ## | |
| 52 ## Bedgraph options | |
| 53 #if $bdg_options.bdg | |
| 54 -B $bdg_options.spmr | |
| 55 #end if | |
| 56 ## | |
| 57 ## Advanced options | |
| 58 #if $advanced_options.advanced_options_selector | |
| 59 --mfold $advanced_options.mfoldlo $advanced_options.mfoldhi | |
| 60 $advanced_options.nolambda | |
| 61 $advanced_options.call_summits | |
| 62 #if str($advanced_options.keep_duplicates.keep_dup) == '' | |
| 63 --keep-dup $advanced_options.keep_duplicates.maximum_tags | |
| 64 #else | |
| 65 --keep-dup $advanced_options.keep_duplicates.keep_dup | |
| 66 #end if | |
| 67 #else | |
| 68 ## Defaults if advanced options not set | |
| 69 --mfold 10 30 --keep-dup 1 | |
| 70 #end if | |
| 71 ## | |
| 72 ## Output files | |
| 73 --output-summits=$output_summits_bed_file | |
| 74 --output-extra-files=$output_extra_files | |
| 75 --output-extra-files-path=$output_extra_files.files_path | |
| 76 ## | |
| 77 ## Narrow/broad peak outputs | |
| 78 #if str($broad_options.broad_regions) == 'broad' | |
| 79 --output-broadpeaks=$output_broadpeaks_file | |
| 80 --output-gappedpeaks=$output_gappedpeaks_file | |
| 81 #else | |
| 82 --output-narrowpeaks=$output_narrowpeaks_file | |
| 83 #end if | |
| 84 ## | |
| 85 ## Bedgraph outputs | |
| 86 #if $bdg_options.bdg | |
| 87 --output-pileup=$output_treat_pileup_file | |
| 88 --output-lambda-bedgraph=$output_lambda_bedgraph_file | |
| 89 #if $bdg_options.make_bigwig | |
| 90 --output-bigwig=$output_bigwig_file | |
| 91 --length=$GALAXY_DATA_INDEX_DIR/shared/ucsc/chrom/${input_chipseq_file1.dbkey}.len | |
| 92 #end if | |
| 93 #end if | |
| 94 ## | |
| 95 ## XLS/interval output | |
| 96 #if str($xls_to_interval) == 'True' | |
| 97 --output-xls-to-interval=$output_xls_to_interval_peaks_file | |
| 98 #else | |
| 99 --output-peaks=$output_peaks_file | |
| 100 #end if | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
101 ]]></command> |
| 0 | 102 <inputs> |
| 103 <!--experiment name used as base for output file names --> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
104 <param name="experiment_name" type="text" value="MACS2.1 in Galaxy" size="50" |
| 0 | 105 label="Experiment Name"/> |
| 106 <!--choose 'broad' or 'narrow' regions--> | |
| 107 <conditional name="broad_options"> | |
| 108 <param name="broad_regions" type="select" label="Type of region to call" | |
| 109 help="Broad regions are formed by linking nearby enriched regions"> | |
| 110 <option value="" selected="true">Narrow regions</option> | |
| 111 <option value="broad">Broad regions</option> | |
| 112 </param> | |
| 113 <when value="broad"> | |
| 114 <param name="broad_cutoff" type="float" | |
| 115 label="Cutoff for broad regions" | |
| 116 value="0.1" help="default: 0.1 (--broad-cutoff)"/> | |
| 117 </when> | |
| 118 </conditional> | |
| 1 | 119 <param name="format" type="select" label="Format of input read data" |
| 120 help="Specify the format of the input data and whether or not it is paired end (--format)"> | |
| 121 <option value="BAMPE" selected="true">BAM (paired-end)</option> | |
| 122 <option value="BAM">BAM (single-end)</option> | |
| 123 <option value="BEDPE">BED (paired-end)</option> | |
| 124 <option value="BED">BED (single-end)</option> | |
| 125 <option value="SAMPE">SAM (paired-end)</option> | |
| 126 <option value="SAM">SAM (single-end)</option> | |
| 127 </param> | |
| 0 | 128 <param name="input_chipseq_file1" type="data" format="bed,sam,bam" |
| 129 label="ChIP-seq read file" /> | |
| 130 <param name="input_control_file1" type="data" format="bed,sam,bam" optional="True" | |
| 131 label="ChIP-seq control read file" /> | |
| 132 <conditional name="genome_size"> | |
| 133 <param name="gsize" type="select" label="Effective genome size" | |
| 134 help="Either pre-defined (for common organisms), or user-defined (--gsize)"> | |
| 135 <option value="hs" selected="true">Human (2.7e9)</option> | |
| 136 <option value="mm">Mouse (1.87e9)</option> | |
| 137 <option value="ce">C. elegans (9e7)</option> | |
| 138 <option value="dm">Fruitfly (1.2e8)</option> | |
| 139 <option value="">User-defined</option> | |
| 140 </param> | |
| 141 <when value=""> | |
| 142 <!-- User-defined effective genome size --> | |
| 143 <param name="user_defined_gsize" type="float" value="" | |
| 144 label="Enter effective genome size (number of bases)" | |
| 145 help="e.g. '1.0e+9' or '1000000000'" /> | |
| 146 </when> | |
| 147 </conditional> | |
| 148 <param name="bw" type="integer" label="Band width" value="300" help="(--bw)"/> | |
| 149 <param name="xls_to_interval" label="Include XLS file from MACS" | |
| 150 type="boolean" truevalue="True" falsevalue="False" checked="True" | |
| 151 help="MACS2 XLS file will be output to the history in 'interval' format (suitable for subsequent analysis in Galaxy). Note that start positions are 1-based."/> | |
| 152 | |
| 153 <conditional name="bdg_options"> | |
| 154 <param name="bdg" | |
| 155 label="Save treatment and control lambda pileups in bedGraph" | |
| 156 type="boolean" truevalue="-B" falsevalue="" checked="False" /> | |
| 157 <when value="-B"> | |
| 158 <param name="spmr" | |
| 159 type="boolean" truevalue="--SPMR" falsevalue="" checked="False" | |
| 160 label="Save signal per million reads for fragment pileup profiles" | |
| 161 help="(--SPMR)" /> | |
| 162 <param name="make_bigwig" type="boolean" checked="True" | |
| 163 truevalue="True" falsevalue="" | |
| 164 label="Also generate bigWig file from bedGraph" | |
| 165 help="bigWig file can used in subsequent analyses e.g. CEAS" /> | |
| 166 </when> | |
| 167 <when value=""> | |
| 168 <!-- Display nothing --> | |
| 169 </when> | |
| 170 </conditional> | |
| 171 | |
| 172 <conditional name="pq_options"> | |
| 173 <param name="pq_options_selector" type="select" | |
| 174 label="Select p-value or q-value" help="default uses q-value"> | |
| 175 <option value="qvalue">q-value</option> | |
| 176 <option value="pvalue">p-value</option> | |
| 177 </param> | |
| 178 <when value="pvalue"> | |
| 179 <param name="pvalue" type="float" | |
| 180 label="p-value cutoff for binding region detection" | |
| 181 value="1e-2" help="default: 1e-2 (--pvalue)"/> | |
| 182 </when> | |
| 183 <when value="qvalue"> | |
| 184 <param name="qvalue" type="float" | |
| 185 label="q-value cutoff for binding region detection" | |
| 186 value="0.01" help="default: 0.01 (--qvalue)"/> | |
| 187 </when> | |
| 188 </conditional> | |
| 189 <conditional name="advanced_options"> | |
| 190 <param name="advanced_options_selector" | |
| 191 type="boolean" truevalue="on" falsevalue="off" checked="False" | |
| 192 label="Use advanced options?" /> | |
| 193 <when value="on"> | |
| 194 <param name="mfoldlo" type="integer" | |
| 195 label="Select the regions with MFOLD high-confidence enrichment ratio against background to build model (lower-limit)" | |
| 196 value="10" help="(--mfold)"/> | |
| 197 <param name="mfoldhi" type="integer" | |
| 198 label="Select the regions with MFOLD high-confidence enrichment ratio against background to build model (upper-limit)" | |
| 199 value="30" help="(--mfold)"/> | |
| 200 <param name="nolambda" | |
| 201 label="Use fixed background lambda as local lambda for every binding region" | |
| 202 type="boolean" truevalue="--nolambda" falsevalue="" checked="False" | |
| 203 help="(--nolambda)"/> | |
| 204 <param name="call_summits" | |
| 205 label="Detect subpeaks within binding region" | |
| 206 type="boolean" truevalue="--call-summits" falsevalue="" checked="False" | |
| 207 help="(--call-summits)"/> | |
| 208 <conditional name="keep_duplicates"> | |
| 209 <param name="keep_dup" type="select" | |
| 210 label="Use of duplicate reads"> | |
| 211 <option value="auto">Automatically calculate maximum number of duplicates to keep (auto)</option> | |
| 212 <option value="all">Use all duplicates (all)</option> | |
| 213 <option value="" selected="true">Manually specify maxium number of duplicates</option> | |
| 214 </param> | |
| 215 <when value=""> | |
| 216 <param name="maximum_tags" type="integer" value="1" | |
| 217 label="Maxium number of duplicated tags to keep at each location"/> | |
| 218 </when> | |
| 219 </conditional> | |
| 220 </when> | |
| 221 <when value="off"> | |
| 222 <!--display nothing--> | |
| 223 </when> | |
| 224 </conditional> | |
| 225 <conditional name="nomodel_type"> | |
| 226 <param name="nomodel_type_selector" type="select" label="Build Model"> | |
| 227 <option value="nomodel">Do not build the shifting model (--nomodel enabled)</option> | |
| 228 <option value="create_model" selected="true">Build the shifting model (--nomodel disabled)</option> | |
| 229 </param> | |
| 230 <when value="nomodel"> | |
| 231 <param name="extsize" type="integer" label="Arbitrary extension size in bp" value="200" help="Used as fragment size to extend each read towards 3' end (--extsize)"/> | |
| 232 </when> | |
| 233 </conditional> | |
| 234 </inputs> | |
| 235 | |
| 236 <outputs> | |
| 237 <!--callpeaks output--> | |
| 238 <data name="output_extra_files" format="html" | |
| 239 label="${tool.name}: callpeak on ${on_string} (html report)"> | |
| 240 </data> | |
| 241 <data name="output_summits_bed_file" format="bed" | |
| 242 label="${tool.name}: callpeak on ${on_string} (summits: bed)"> | |
| 243 </data> | |
| 244 <data name="output_peaks_file" format="xls" | |
| 245 label="${tool.name}: callpeak on ${on_string} (peaks: xls)"> | |
| 246 <filter>xls_to_interval is False</filter> | |
| 247 </data> | |
| 248 <data name="output_narrowpeaks_file" format="interval" | |
| 249 label="${tool.name}: callpeak on ${on_string} (peaks: narrowPeak)"> | |
| 250 <filter>broad_options['broad_regions'] == ''</filter> | |
| 251 </data> | |
| 252 <data name="output_broadpeaks_file" format="interval" | |
| 253 label="${tool.name}: callpeak on ${on_string} (peaks: broadPeak)"> | |
| 254 <filter>broad_options['broad_regions'] == 'broad'</filter> | |
| 255 </data> | |
| 256 <data name="output_gappedpeaks_file" format="interval" | |
| 257 label="${tool.name}: callpeak on ${on_string} (peaks: gappedPeak)"> | |
| 258 <filter>broad_options['broad_regions'] == 'broad'</filter> | |
| 259 </data> | |
| 260 <data name="output_xls_to_interval_peaks_file" format="interval" | |
| 261 label="${tool.name}: callpeak on ${on_string} (peaks: interval)"> | |
| 262 <filter>xls_to_interval is True</filter> | |
| 263 </data> | |
| 264 <data name="output_treat_pileup_file" format="bedgraph" | |
| 265 label="${tool.name}: callpeak on ${on_string} (treat pileup: bedGraph)"> | |
| 266 <filter>bdg_options['bdg'] is True</filter> | |
| 267 </data> | |
| 268 <data name="output_lambda_bedgraph_file" format="bedgraph" | |
| 269 label="${tool.name}: callpeak on ${on_string} (control lambda: bedGraph)"> | |
| 270 <filter>bdg_options['bdg'] is True</filter> | |
| 271 </data> | |
| 272 <data name="output_bigwig_file" format="bigwig" | |
| 273 label="${tool.name}: callpeak on ${on_string} (treat pileup: bigWig)"> | |
| 274 <filter>bdg_options['bdg'] is True</filter> | |
| 275 <filter>bdg_options['make_bigwig'] is True</filter> | |
| 276 </data> | |
| 277 </outputs> | |
| 278 <tests> | |
| 279 <!-- Peak calling without bigwig output --> | |
| 280 <test> | |
| 281 <!-- Inputs --> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
282 <param name="experiment_name" value="test_MACS2.1.1" /> |
| 0 | 283 <param name="broad_regions" value="" /> |
| 1 | 284 <param name="format" value="BED" /> |
| 0 | 285 <param name="input_chipseq_file1" value="test_region_IP.bed" dbkey="galGal3" |
| 286 ftype="bed" /> | |
| 287 <param name="input_control_file1" value="test_region_Input.bed" | |
| 288 ftype="bed" /> | |
| 289 <param name="gsize" value="" /> | |
| 290 <param name="user_defined_gsize" value="775000000.0" /> | |
| 291 <param name="bw" value="300" /> | |
| 292 <param name="xls_to_interval" value="true" /> | |
| 293 <param name="bdg_options|bdg" value="-B" /> | |
| 294 <param name="bdg_options|spmr" value="--SPMR" /> | |
| 295 <param name="bdg_options|make_bigwig" value="false" /> | |
| 296 <param name="pq_options_selector" value="qvalue" /> | |
| 297 <param name="qvalue" value="0.05" /> | |
| 298 <param name="advanced_options_selector" value="true" /> | |
| 299 <param name="advanced_options|mfoldlo" value="5" /> | |
| 300 <param name="advanced_options|mfoldhi" value="50" /> | |
| 301 <param name="advanced_options|nolambda" value="" /> | |
| 302 <param name="advanced_options|call_summits" value="" /> | |
| 303 <param name="advanced_options|keep_duplicates" value="" /> | |
| 304 <param name="advanced_options|maximum_tags" value="1" /> | |
| 305 <param name="nomodel_type_selector" value="nomodel" /> | |
| 306 <param name="nomodel_type|extsize" value="243" /> | |
| 307 <!-- Outputs --> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
308 <output name="output_extra_files" file="test_MACS2.1.1_html_report.zip" |
| 0 | 309 compare="sim_size" delta="1500" /> |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
310 <output name="output_summits_bed_file" file="test_MACS2.1.1_summits.bed" /> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
311 <output name="output_narrowpeaks_file" file="test_MACS2.1.1_peaks_narrowPeak.interval" /> |
| 0 | 312 <output name="output_xls_to_interval_peaks_file" |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
313 file="test_MACS2.1.1_peaks.xls.re_match" |
| 0 | 314 compare="re_match" lines_diff="1" /> |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
315 <output name="output_treat_pileup_file" file="test_MACS2.1.1_treat_pileup.bdg" /> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
316 <output name="output_lambda_bedgraph_file" file="test_MACS2.1.1_control_lambda.bdg" /> |
| 0 | 317 </test> |
| 318 <!-- Peak calling with bigwig output --> | |
| 319 <test> | |
| 320 <!-- Inputs --> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
321 <param name="experiment_name" value="test_MACS2.1.1" /> |
| 0 | 322 <param name="broad_regions" value="" /> |
| 1 | 323 <param name="format" value="BED" /> |
| 0 | 324 <param name="input_chipseq_file1" value="test_region_IP.bed" dbkey="galGal3" |
| 325 ftype="bed" /> | |
| 326 <param name="input_control_file1" value="test_region_Input.bed" | |
| 327 ftype="bed" /> | |
| 328 <param name="gsize" value="" /> | |
| 329 <param name="user_defined_gsize" value="775000000.0" /> | |
| 330 <param name="bw" value="300" /> | |
| 331 <param name="xls_to_interval" value="true" /> | |
| 332 <param name="bdg_options|bdg" value="-B" /> | |
| 333 <param name="bdg_options|spmr" value="--SPMR" /> | |
| 334 <param name="bdg_options|make_bigwig" value="true" /> | |
| 335 <param name="pq_options_selector" value="qvalue" /> | |
| 336 <param name="qvalue" value="0.05" /> | |
| 337 <param name="advanced_options_selector" value="true" /> | |
| 338 <param name="advanced_options|mfoldlo" value="5" /> | |
| 339 <param name="advanced_options|mfoldhi" value="50" /> | |
| 340 <param name="advanced_options|nolambda" value="" /> | |
| 341 <param name="advanced_options|call_summits" value="" /> | |
| 342 <param name="advanced_options|keep_duplicates" value="" /> | |
| 343 <param name="advanced_options|maximum_tags" value="1" /> | |
| 344 <param name="nomodel_type_selector" value="nomodel" /> | |
| 345 <param name="nomodel_type|extsize" value="243" /> | |
| 346 <!-- Outputs --> | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
347 <output name="output_extra_files" file="test_MACS2.1.1_bw_html_report.zip" |
| 0 | 348 compare="sim_size" delta="2500" /> |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
349 <output name="output_summits_bed_file" file="test_MACS2.1.1_summits.bed" /> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
350 <output name="output_narrowpeaks_file" file="test_MACS2.1.1_peaks_narrowPeak.interval" /> |
| 0 | 351 <output name="output_xls_to_interval_peaks_file" |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
352 file="test_MACS2.1.1_peaks.xls.re_match" |
| 0 | 353 compare="re_match" lines_diff="1" /> |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
354 <output name="output_treat_pileup_file" file="test_MACS2.1.1_treat_pileup.bdg" /> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
355 <output name="output_lambda_bedgraph_file" file="test_MACS2.1.1_control_lambda.bdg" /> |
|
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
356 <output name="output_bigwig_file" file="test_MACS2.1.1_treat_pileup.bw" |
| 0 | 357 compare="sim_size" /> |
| 358 </test> | |
| 359 </tests> | |
| 360 <help> | |
| 361 **What it does** | |
| 362 | |
|
3
4124781932db
Updated to MACS 2.1.1 and use conda for dependency resolution.
pjbriggs
parents:
2
diff
changeset
|
363 MACS (Model-based Analysis of ChIP-seq) 2.1.1 provides algorithms for identifying |
| 0 | 364 transcript factor binding sites. The program can be used either for ChIP-Seq data alone, |
| 365 or with control sample data to improve specificity. | |
| 366 | |
| 367 View the MACS2 documentation at: | |
| 368 https://github.com/taoliu/MACS/blob/master/README.rst | |
| 369 | |
| 370 ------ | |
| 371 | |
| 372 **Usage** | |
| 373 | |
| 374 The tool interfaces with the **callpeak** function in MACS, which calls peaks from | |
| 375 alignment results. | |
| 376 | |
| 377 ------ | |
| 378 | |
| 379 **Credits** | |
| 380 | |
| 381 This Galaxy tool was based on the MACS2 tool hosted in the Galaxy toolshed at | |
| 382 | |
| 383 * http://toolshed.g2.bx.psu.edu/view/modencode-dcc/macs2 | |
| 384 | |
| 385 (specifically the 16:14f378e35191 revision of the tool) which is credited to Ziru | |
| 386 Zhou. This version is a reimplemented version developed within the Bioinformatics | |
| 387 Core Facility at the University of Manchester, which uses more up-to-date Galaxy | |
| 388 syntax and adds some extra features. | |
| 389 | |
| 390 The tool runs Tao Liu's MACS2 software: | |
| 391 | |
| 392 * https://github.com/taoliu/MACS | |
| 393 | |
| 394 The reference for MACS is: | |
| 395 | |
| 396 * Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, | |
| 397 Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS). | |
| 398 Genome Biol. 2008;9(9):R137. | |
| 399 | |
| 400 Please kindly acknowledge both this Galaxy tool and the MACS2 package if you | |
| 401 use it. | |
| 402 </help> | |
| 403 <citations> | |
| 404 <!-- | |
| 405 See https://wiki.galaxyproject.org/Admin/Tools/ToolConfigSyntax#A.3Ccitations.3E_tag_set | |
| 406 Can be either DOI or Bibtex | |
| 407 Use http://www.bioinformatics.org/texmed/ to convert PubMed to Bibtex | |
| 408 --> | |
| 409 <citation type="doi">10.1186/gb-2008-9-9-r137</citation> | |
| 410 </citations> | |
| 411 </tool> |