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