Mercurial > repos > pjbriggs > macs21
annotate macs21_wrapper.xml @ 7:0c6b14f3fefc draft
Attempt to fix installation of tool dependencies for UCSC/bigwig generation, and trap for undefined genome build.
| author | pjbriggs |
|---|---|
| date | Tue, 21 Apr 2015 08:29:16 -0400 |
| parents | 344dd37d1704 |
| children | 78c15c0a96ae |
| rev | line source |
|---|---|
| 2 | 1 <tool id="macs2_1_peakcalling" name="MACS2.1.0" version="2.1.0-2"> |
| 0 | 2 <requirements> |
| 3 <requirement type="package" version="2.7">python</requirement> | |
| 2 | 4 <requirement type="package" version="1.9">numpy</requirement> |
| 0 | 5 <requirement type="package" version="2.1.0.20140616">macs2</requirement> |
| 6 | 6 <requirement type="package" version="3.1.2">R</requirement> |
|
5
9f6d5cae9181
Uploaded yet another attempt to fix tool dependencies!
pjbriggs
parents:
4
diff
changeset
|
7 <requirement type="package" version="1.0">ucsc_tools_for_macs21</requirement> |
| 0 | 8 </requirements> |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
9 <description>Model-based Analysis of ChIP-Seq</description> |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
10 <command interpreter="python"> |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
11 macs21_wrapper.py |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
12 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
13 ## Major command |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
14 $major_command.major_command_selector |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
15 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
16 ## ChIP-seq input |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
17 $major_command.input_chipseq_file1 |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
18 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
19 ## ChIP-seq control |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
20 #if str($major_command.input_control_file1) != 'None' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
21 -c $major_command.input_control_file1 |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
22 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
23 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
24 ## Call peaks |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
25 #if str($major_command.major_command_selector) == 'callpeak' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
26 --format=$major_command.input_chipseq_file1.extension |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
27 --name="$experiment_name" |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
28 --bw=$major_command.bw |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
29 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
30 ## Genome size |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
31 #if str($major_command.genome_size.gsize) == '' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
32 --gsize=$major_command.genome_size.user_defined_gsize |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
33 #else: |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
34 --gsize=$major_command.genome_size.gsize |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
35 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
36 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
37 ## Broad peaks |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
38 #if str($major_command.broad_options.broad_regions) == 'broad' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
39 --broad --broad-cutoff=$major_command.broad_options.broad_cutoff |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
40 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
41 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
42 ## (no)model options |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
43 #if str($major_command.nomodel_type.nomodel_type_selector) == 'nomodel' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
44 --nomodel --extsize=$major_command.nomodel_type.extsize |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
45 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
46 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
47 ## pq value select options |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
48 #if str($major_command.pq_options.pq_options_selector) == 'qvalue' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
49 --qvalue=$major_command.pq_options.qvalue |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
50 #else |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
51 --pvalue=$major_command.pq_options.pvalue |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
52 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
53 ## |
| 2 | 54 ## Bedgraph options |
| 55 #if $major_command.bdg_options.bdg == True | |
| 56 -B $major_command.bdg_options.spmr | |
| 57 #end if | |
| 58 ## | |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
59 ## Advanced options |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
60 #if str($major_command.advanced_options.advanced_options_selector) == 'on' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
61 --mfold $major_command.advanced_options.mfoldlo $major_command.advanced_options.mfoldhi |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
62 $major_command.advanced_options.nolambda |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
63 $major_command.advanced_options.call_summits |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
64 #if str($major_command.advanced_options.keep_duplicates.keep_dup) == '' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
65 --keep-dup $major_command.advanced_options.keep_duplicates.maximum_tags |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
66 #else |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
67 --keep-dup $major_command.advanced_options.keep_duplicates.keep_dup |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
68 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
69 #else |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
70 ## Defaults if advanced options not set |
| 2 | 71 --mfold 10 30 --keep-dup 1 |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
72 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
73 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
74 ## Output files |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
75 --output-summits=$output_summits_bed_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
76 --output-extra-files=$output_extra_files |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
77 --output-extra-files-path=$output_extra_files.files_path |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
78 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
79 ## Narrow/broad peak outputs |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
80 #if str($major_command.broad_options.broad_regions) == 'broad' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
81 --output-broadpeaks=$output_broadpeaks_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
82 --output-gappedpeaks=$output_gappedpeaks_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
83 #else |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
84 --output-narrowpeaks=$output_narrowpeaks_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
85 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
86 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
87 ## Bedgraph outputs |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
88 #if str($major_command.bdg_options.bdg) == 'True' |
| 2 | 89 --output-pileup=$output_treat_pileup_file |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
90 --output-lambda-bedgraph=$output_lambda_bedgraph_file |
| 2 | 91 #if str($major_command.bdg_options.make_bigwig) == 'True' |
| 92 --output-bigwig=$output_bigwig_file | |
| 93 --length=$GALAXY_DATA_INDEX_DIR/shared/ucsc/chrom/${major_command.input_chipseq_file1.dbkey}.len | |
| 94 #end if | |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
95 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
96 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
97 ## XLS/interval output |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
98 #if str($major_command.xls_to_interval) == 'True' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
99 --output-xls-to-interval=$output_xls_to_interval_peaks_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
100 #else |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
101 --output-peaks=$output_peaks_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
102 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
103 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
104 ## |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
105 ## Compare .bdg files |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
106 #if str($major_command.major_command_selector) == 'bdgcmp' |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
107 -m $major_command.bdgcmp_options.bdgcmp_options_selector |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
108 -p $major_command.pseudocount |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
109 --output-bdgcmp $output_bdgcmp_file |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
110 #end if |
|
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
111 </command> |
| 0 | 112 <inputs> |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
113 <!--experiment name used as base for output file names --> |
| 0 | 114 <param name="experiment_name" type="text" value="MACS2.1.0 in Galaxy" size="50" |
| 115 label="Experiment Name"/> | |
|
1
d0986d2be693
Substantial reimplementation of internals, also renamed id and version.
pjbriggs
parents:
0
diff
changeset
|
116 <!--select a major MACS2 command--> |
| 0 | 117 <conditional name="major_command"> |
| 118 <param name="major_command_selector" type="select" label="Select action to be performed"> | |
| 119 <option value="callpeak">Peak Calling</option> | |
| 120 <option value="bdgcmp">Compare .bdg Files</option> | |
| 121 </param> | |
| 122 <!--callpeak option of macs2--> | |
| 123 <when value="callpeak"> | |
| 124 <!--choose 'broad' or 'narrow' regions--> | |
| 125 <conditional name="broad_options"> | |
| 126 <param name="broad_regions" type="select" label="Type of region to call" | |
| 127 help="Broad regions are formed by linking nearby enriched regions"> | |
| 128 <option value="" selected="true">Narrow regions</option> | |
| 129 <option value="broad">Broad regions</option> | |
| 130 </param> | |
| 131 <when value="broad"> | |
| 132 <param name="broad_cutoff" type="float" | |
| 133 label="Cutoff for broad regions" | |
| 134 value="0.1" help="default: 0.1 (--broad-cutoff)"/> | |
| 135 </when> | |
| 136 </conditional> | |
| 137 <param name="input_chipseq_file1" type="data" format="bed,sam,bam" | |
| 138 label="ChIP-seq read file" /> | |
| 139 <param name="input_control_file1" type="data" format="bed,sam,bam" optional="True" | |
| 140 label="ChIP-seq control read file" /> | |
| 141 <conditional name="genome_size"> | |
| 142 <param name="gsize" type="select" label="Effective genome size" | |
| 143 help="Either pre-defined (for common organisms), or user-defined (--gsize)"> | |
| 144 <option value="hs" selected="true">Human (2.7e9)</option> | |
| 145 <option value="mm">Mouse (1.87e9)</option> | |
| 146 <option value="ce">C. elegans (9e7)</option> | |
| 147 <option value="dm">Fruitfly (1.2e8)</option> | |
| 148 <option value="">User-defined</option> | |
| 149 </param> | |
| 150 <when value=""> | |
| 151 <!-- User-defined effective genome size --> | |
| 152 <param name="user_defined_gsize" type="float" value="" | |
| 153 label="Enter effective genome size (number of bases)" | |
| 154 help="e.g. '1.0e+9' or '1000000000'" /> | |
| 155 </when> | |
| 156 </conditional> | |
| 157 <param name="bw" type="integer" label="Band width" value="300" help="(--bw)"/> | |
| 158 <param name="xls_to_interval" label="Include XLS file from MACS" | |
| 159 type="boolean" truevalue="True" falsevalue="False" checked="True" | |
| 160 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."/> | |
| 161 | |
| 162 <conditional name="bdg_options"> | |
| 163 <param name="bdg" | |
| 2 | 164 label="Save treatment and control lambda pileups in bedGraph" |
| 0 | 165 type="boolean" truevalue="-B" falsevalue="" checked="False" /> |
| 166 <when value="-B"> | |
| 167 <param name="spmr" | |
| 168 type="boolean" truevalue="--SPMR" falsevalue="" checked="False" | |
| 169 label="Save signal per million reads for fragment pileup profiles" | |
| 170 help="(--SPMR)" /> | |
| 2 | 171 <param name="make_bigwig" type="boolean" checked="True" |
| 172 truevalue="True" falsevalue="" | |
| 173 label="Also generate bigWig file from bedGraph" | |
| 174 help="bigWig file can used in subsequent analyses e.g. CEAS" /> | |
| 0 | 175 </when> |
| 176 <when value=""> | |
| 177 <!-- Display nothing --> | |
| 178 </when> | |
| 179 </conditional> | |
| 180 | |
| 181 <conditional name="pq_options"> | |
| 182 <param name="pq_options_selector" type="select" | |
| 183 label="Select p-value or q-value" help="default uses q-value"> | |
| 184 <option value="qvalue">q-value</option> | |
| 185 <option value="pvalue">p-value</option> | |
| 186 </param> | |
| 187 <when value="pvalue"> | |
| 188 <param name="pvalue" type="float" | |
| 189 label="p-value cutoff for binding region detection" | |
| 190 value="1e-2" help="default: 1e-2 (--pvalue)"/> | |
| 191 </when> | |
| 192 <when value="qvalue"> | |
| 193 <param name="qvalue" type="float" | |
| 194 label="q-value cutoff for binding region detection" | |
| 195 value="0.01" help="default: 0.01 (--qvalue)"/> | |
| 196 </when> | |
| 197 </conditional> | |
| 198 <conditional name="advanced_options"> | |
| 199 <param name="advanced_options_selector" type="select" | |
| 200 label="Display advanced options"> | |
| 201 <option value="off">Hide</option> | |
| 202 <option value="on">Display</option> | |
| 203 </param> | |
| 204 <when value="on"> | |
| 205 <param name="mfoldlo" type="integer" | |
| 206 label="Select the regions with MFOLD high-confidence enrichment ratio against background to build model (lower-limit)" | |
| 207 value="10" help="(--mfold)"/> | |
| 208 <param name="mfoldhi" type="integer" | |
| 209 label="Select the regions with MFOLD high-confidence enrichment ratio against background to build model (upper-limit)" | |
| 210 value="30" help="(--mfold)"/> | |
| 211 <param name="nolambda" | |
| 212 label="Use fixed background lambda as local lambda for every binding region" | |
| 213 type="boolean" truevalue="--nolambda" falsevalue="" checked="False" | |
| 214 help="(--nolambda)"/> | |
| 215 <param name="call_summits" | |
| 216 label="Detect subpeaks within binding region" | |
| 217 type="boolean" truevalue="--call-summits" falsevalue="" checked="False" | |
| 218 help="(--call-summits)"/> | |
| 219 <conditional name="keep_duplicates"> | |
| 220 <param name="keep_dup" type="select" | |
| 221 label="Use of duplicate reads"> | |
| 222 <option value="auto">Automatically calculate maximum number of duplicates to keep (auto)</option> | |
| 223 <option value="all">Use all duplicates (all)</option> | |
| 224 <option value="" selected="true">Manually specify maxium number of duplicates</option> | |
| 225 </param> | |
| 226 <when value=""> | |
| 227 <param name="maximum_tags" type="integer" value="1" | |
| 228 label="Maxium number of duplicated tags to keep at each location"/> | |
| 229 </when> | |
| 230 </conditional> | |
| 231 </when> | |
| 232 <when value="off"> | |
| 233 <!--display nothing--> | |
| 234 </when> | |
| 235 </conditional> | |
| 236 <conditional name="nomodel_type"> | |
| 237 <param name="nomodel_type_selector" type="select" label="Build Model"> | |
| 238 <option value="nomodel">Do not build the shifting model (--nomodel enabled)</option> | |
| 239 <option value="create_model" selected="true">Build the shifting model (--nomodel disabled)</option> | |
| 240 </param> | |
| 241 <when value="nomodel"> | |
| 242 <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)"/> | |
| 243 </when> | |
| 244 </conditional> | |
| 245 </when> | |
| 246 | |
| 247 <!--callpeak option of macs2--> | |
| 248 <when value="bdgcmp"> | |
| 249 <param name="input_chipseq_file1" type="data" format="bed,sam,bam" | |
| 250 label="ChIP-seq read file" /> | |
| 251 <param name="input_control_file1" type="data" format="bed,sam,bam" optional="True" | |
| 252 label="ChIP-seq control read file" /> | |
| 253 <param name="pseudocount" type="float" label="Set pseudocount" value="0.00001" | |
| 254 help="default: 0.00001 (-p)"/> | |
| 255 <conditional name="bdgcmp_options"> | |
| 256 <param name="bdgcmp_options_selector" type="select" | |
| 257 label="Select action to be performed"> | |
| 258 <option value="ppois">ppois</option> | |
| 259 <option value="qpois">qpois</option> | |
| 260 <option value="subtract">subtract</option> | |
| 261 <option value="logFE">logFE</option> | |
| 262 <option value="FE">FE</option> | |
| 263 <option value="logLR">logLR</option> | |
| 264 </param> | |
| 265 </conditional> | |
| 266 </when> | |
| 267 </conditional> | |
| 268 </inputs> | |
| 269 | |
| 270 <outputs> | |
| 271 <!--callpeaks output--> | |
| 272 <data name="output_extra_files" format="html" | |
| 273 label="${tool.name}: callpeak on ${on_string} (html report)"> | |
| 274 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 275 </data> | |
| 276 <data name="output_summits_bed_file" format="bed" | |
| 277 label="${tool.name}: callpeak on ${on_string} (summits: bed)"> | |
| 278 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 279 </data> | |
| 280 <data name="output_peaks_file" format="xls" | |
| 281 label="${tool.name}: callpeak on ${on_string} (peaks: xls)"> | |
| 282 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 283 <filter>major_command['xls_to_interval'] is False</filter> | |
| 284 </data> | |
| 285 <data name="output_narrowpeaks_file" format="interval" | |
| 286 label="${tool.name}: callpeak on ${on_string} (peaks: narrowPeak)"> | |
| 287 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 288 <filter>major_command['broad_options']['broad_regions'] == ''</filter> | |
| 289 </data> | |
| 290 <data name="output_broadpeaks_file" format="interval" | |
| 291 label="${tool.name}: callpeak on ${on_string} (peaks: broadPeak)"> | |
| 292 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 293 <filter>major_command['broad_options']['broad_regions'] == 'broad'</filter> | |
| 294 </data> | |
| 295 <data name="output_gappedpeaks_file" format="interval" | |
| 296 label="${tool.name}: callpeak on ${on_string} (peaks: gappedPeak)"> | |
| 297 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 298 <filter>major_command['broad_options']['broad_regions'] == 'broad'</filter> | |
| 299 </data> | |
| 300 <data name="output_xls_to_interval_peaks_file" format="interval" | |
| 301 label="${tool.name}: callpeak on ${on_string} (peaks: interval)"> | |
| 302 <filter>major_command['xls_to_interval'] is True</filter> | |
| 303 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 304 </data> | |
| 305 <data name="output_treat_pileup_file" format="bedgraph" | |
| 306 label="${tool.name}: callpeak on ${on_string} (treat pileup: bedGraph)"> | |
| 307 <filter>major_command['bdg_options']['bdg'] is True</filter> | |
| 308 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 309 </data> | |
| 310 <data name="output_lambda_bedgraph_file" format="bedgraph" | |
| 311 label="${tool.name}: callpeak on ${on_string} (control lambda: bedGraph)"> | |
| 312 <filter>major_command['bdg_options']['bdg'] is True</filter> | |
| 313 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 314 </data> | |
| 2 | 315 <data name="output_bigwig_file" format="bigwig" |
| 316 label="${tool.name}: callpeak on ${on_string} (treat pileup: bigWig)"> | |
| 317 <filter>major_command['major_command_selector'] == 'callpeak'</filter> | |
| 318 <filter>major_command['bdg_options']['bdg'] is True</filter> | |
| 319 <filter>major_command['bdg_options']['make_bigwig'] is True</filter> | |
| 320 </data> | |
| 0 | 321 <!--bdgcmp output--> |
| 322 <data name="output_bdgcmp_file" format="bdg" | |
| 323 label="${tool.name}: bdgcmp on ${on_string} (bdg)"> | |
| 324 <filter>major_command['major_command_selector'] == 'bdgcmp'</filter> | |
| 325 </data> | |
| 326 </outputs> | |
| 327 <tests> | |
| 328 <!--none yet for macs2--> | |
| 329 </tests> | |
| 330 <help> | |
| 331 | |
| 332 .. class:: warningmark | |
| 333 | |
| 334 **This is a modified version of the standard Galaxy toolshed "MACS2" tool, | |
| 335 which has been customised for users at the University of Manchester to work | |
| 336 with MACS 2.1.0.** | |
| 337 | |
| 338 It is based on the 16:14f378e35191 revision of the tool at | |
| 339 | |
| 340 * http://toolshed.g2.bx.psu.edu/view/modencode-dcc/macs2 | |
| 341 | |
| 342 ------ | |
| 343 | |
| 344 **What it does** | |
| 345 | |
| 346 With the improvement of sequencing techniques, chromatin immunoprecipitation | |
| 347 followed by high throughput sequencing (ChIP-Seq) is getting popular to study | |
| 348 genome-wide protein-DNA interactions. To address the lack of powerful ChIP-Seq | |
| 349 analysis method, we present a novel algorithm, named Model-based Analysis of | |
| 350 ChIP-Seq (MACS), for identifying transcript factor binding sites. MACS captures | |
| 351 the influence of genome complexity to evaluate the significance of enriched | |
| 352 ChIP regions, and MACS improves the spatial resolution of binding sites through | |
| 353 combining the information of both sequencing tag position and orientation. MACS | |
| 354 can be easily used for ChIP-Seq data alone, or with control sample with the | |
| 355 increase of specificity. | |
| 356 | |
| 357 View the original MACS2 documentation: | |
| 358 https://github.com/taoliu/MACS/blob/master/README.rst | |
| 359 | |
| 360 ------ | |
| 361 | |
| 362 **Usage** | |
| 363 | |
| 364 **Peak Calling**: Main MACS2 Function to Call peaks from alignment results. | |
| 365 | |
| 366 **Compare .bdg files**: Deduct noise by comparing two signal tracks in bedGraph. | |
| 367 | |
| 368 | |
| 369 ------ | |
| 370 | |
| 371 **Citation** | |
| 372 | |
| 373 For the underlying tool, please cite Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137. | |
| 374 | |
| 375 Integration of MACS2 with Galaxy performed by Ziru Zhou ( ziruzhou@gmail.com ). Please send your comments/questions to modENCODE DCC at help@modencode.org. | |
| 376 </help> | |
| 377 </tool> |