0
|
1 <tool id="macs2_bdgpeakcall" name="MACS2 bdgpeakcall" version="@VERSION_STRING@.0">
|
|
2 <description>Call peaks from bedGraph output</description>
|
|
3 <expand macro="requirements">
|
|
4 <requirement type="package" version="4.1.0">gnu_awk</requirement>
|
|
5 </expand>
|
|
6 <expand macro="version_command" />
|
|
7 <macros>
|
|
8 <import>macs2_macros.xml</import>
|
|
9 </macros>
|
|
10 <command>
|
|
11 macs2 bdgpeakcall
|
|
12 --ifile "${ infile }"
|
|
13 --cutoff "${ cutoff }"
|
|
14 --min-length "${ minlen }"
|
|
15 --max-gap "${ maxgap }"
|
|
16 ${ callsummits }
|
|
17 ${ notrackline }
|
|
18 --ofile "macs2_bdgpeakcall.bdg"
|
|
19 &&
|
|
20 awk '!/^track name/' "macs2_bdgpeakcall.bdg" > "${ outfile }"
|
|
21 </command>
|
|
22 <expand macro="stdio" />
|
|
23 <inputs>
|
|
24 <param name="infile" type="data" format="bedgraph" label="MACS score in bedGraph" />
|
|
25
|
|
26 <param name="cutoff" type="float" label="Cutoff for peaks" value="5.0"
|
|
27 help="Cutoff depends on which method you used for score track. If the file contains pvalue scores from MACS2, score 5 means pvalue 1e-5. DEFAULT: 5.0 (--cutoff)"/>
|
|
28 <param name="minlen" type="integer" value="200" label="Minimum length of peak" help="better to set it as d value (--min-length)"/>
|
|
29 <param name="maxgap" type="integer" value="30" label="Maximum gap between significant points in a peak" help="better to set it as tag size DEFAULT: 30 (--max-gap)"/>
|
|
30
|
|
31 <param name="callsummits" type="boolean" truevalue="--call-summits" falsevalue="" checked="False"
|
|
32 label="If set, MACS will use a more sophisticated approach to find all summits in each enriched peak region."
|
|
33 help="(-call-summits)"/>
|
|
34 <param name="notrackline" type="boolean" truevalue="" falsevalue="--no-trackline" checked="False"
|
|
35 label="Include trackline into bedGraph output"
|
|
36 help="required by USCS (--no-trackline)"/>
|
|
37 </inputs>
|
|
38 <outputs>
|
|
39 <data name="outfile" format="tabular" label="${tool.name} on ${on_string}" />
|
|
40 </outputs>
|
|
41 <tests>
|
|
42 <test>
|
|
43 <param name="infile" value="callpeak_treatment_part.bdg" ftype="bedgraph" />
|
|
44 <param name="cutoff" value="5.0"/>
|
|
45 <param name="minlen" value="200"/>
|
|
46 <param name="maxgap" value="30"/>
|
|
47 <output name="outfile" file="bdgpeakcall_on_callpeak_treatment.tabular"/>
|
|
48 </test>
|
|
49 </tests>
|
|
50 <help>
|
|
51 **What it does**
|
|
52
|
|
53 Call peaks from bedGraph output.
|
|
54 bdgpeakcall from macs2
|
|
55
|
|
56
|
|
57 @citation@
|
|
58 </help>
|
|
59 <expand macro="citations" />
|
|
60 </tool>
|