annotate genomeCoverageBed.xml @ 11:7308cc546a36 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 19967671ddd4e750f53d71274e5c12325a402811
author iuc
date Mon, 17 Oct 2016 10:13:35 -0400
parents 8e6b7c3597a8
children 44867b59dbf2
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
11
7308cc546a36 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 19967671ddd4e750f53d71274e5c12325a402811
iuc
parents: 7
diff changeset
1 <tool id="bedtools_genomecoveragebed" name="Genome Coverage" version="@WRAPPER_VERSION@.0">
4
607c0576c6ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
2 <description>compute the coverage over an entire genome</description>
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
3 <macros>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
4 <import>macros.xml</import>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
5 </macros>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
6 <expand macro="requirements" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
7 <expand macro="stdio" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
8 <command>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
9 <![CDATA[
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
10 bedtools genomecov
7
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
11 #if $input_type.input_type_select == 'bam'
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
12 -ibam '$input_type.input'
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
13 #else
7
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
14 -i '$input_type.input'
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
15 -g '$input_type.genome'
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
16 #end if
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
17
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
18 $split
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
19 $strand
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
20
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
21 #if str($report.report_select) == "bg":
2
457b09031d57 Uploaded
iuc
parents: 1
diff changeset
22 #if $report.zero_regions:
457b09031d57 Uploaded
iuc
parents: 1
diff changeset
23 $report.zero_regions
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
24 #else:
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
25 -bg
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
26 #end if
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
27
2
457b09031d57 Uploaded
iuc
parents: 1
diff changeset
28 #if str($report.scale):
457b09031d57 Uploaded
iuc
parents: 1
diff changeset
29 -scale $report.scale
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
30 #end if
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
31 #else:
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
32 #if str($report.max):
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
33 -max $report.max
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
34 #end if
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
35 #end if
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
36 $d
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
37 $dz
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
38 $five
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
39 $three
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
40
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
41 > '$output'
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
42 ]]>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
43 </command>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
44 <inputs>
7
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
45 <conditional name="input_type">
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
46 <param name="input_type_select" type="select" label="Input type">
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
47 <option value="bed">BED/VCF/GFF</option>
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
48 <option value="bam">BAM</option>
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
49 </param>
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
50 <when value="bed">
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
51 <param format="bed,vcf,gff,gff3" name="input" type="data" label="BED/VCF/GFF file" />
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
52 <expand macro="genome" />
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
53 </when>
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
54 <when value="bam">
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
55 <param format="bam" name="input" type="data" label="BAM file" />
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
56 </when>
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
57 </conditional>
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
58 <conditional name="report">
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
59 <param name="report_select" type="select" label="Output type">
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
60 <option value="bg" selected="true">BedGraph coverage file</option>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
61 <option value="hist">Data suiteable for Histogram</option>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
62 </param>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
63 <when value="bg">
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
64 <param name="zero_regions" type="boolean" checked="False" truevalue="-bga" falsevalue=""
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
65 label="Report regions with zero coverage" help="If set, regions without any coverage will also be reported. (-bga)" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
66 <param name="scale" type="float" value="1.0"
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
67 label="Scale the coverage by a constant factor"
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
68 help="Each bedGraph coverage value is multiplied by this factor before being reported. Useful for normalizing coverage by, e.g., reads per million (RPM). (-scale)"/>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
69 </when>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
70 <when value="hist">
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
71 <param name="max" type="integer" label="Specify max depth" value="0"
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
72 help="Combine all positions with a depth >= max into a single bin in the histogram. (-max)"/>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
73 </when>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
74 </conditional>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
75 <expand macro="split" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
76 <param name="strand" type="select" label="Calculate coverage based on" help="(-strand)">
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
77 <option value="">both strands combined</option>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
78 <option value="-strand +">positive strand only</option>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
79 <option value="-strand -">negative strand only</option>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
80 </param>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
81
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
82 <param name="d" type="boolean" checked="False" truevalue="-d" falsevalue=""
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
83 label="Report the depth at each genome position with 1-based coordinates" help="(-d)" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
84 <param name="dz" type="boolean" checked="False" truevalue="-dz" falsevalue=""
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
85 label="Report the depth at each genome position with 0-based coordinatess" help="(-dz)" />
4
607c0576c6ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 6692e3a4fa1bf6e9a407735afdbb2454ed32b316
iuc
parents: 2
diff changeset
86 <param name="five" type="boolean" checked="False" truevalue="-5" falsevalue=""
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
87 label="Calculate coverage of 5’ positions" help="Instead of entire interval. (-5)" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
88 <param name="three" type="boolean" checked="False" truevalue="-3" falsevalue=""
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
89 label="Calculate coverage of 3’ positions" help="Instead of entire interval. (-3)" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
90 </inputs>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
91 <outputs>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
92 <data format="bedgraph" name="output">
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
93 <change_format>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
94 <when input="report.report_select" value="hist" format="tabular" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
95 </change_format>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
96 </data>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
97 </outputs>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
98 <tests>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
99 <test>
7
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
100 <param name="input_type_select" value="bed" />
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
101 <param name="input" value="genomeCoverageBed1.bed" ftype="bed" />
7
8e6b7c3597a8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 964b3159deeade1c90d20cef92dc15a14195edb7
iuc
parents: 4
diff changeset
102 <param name="genome" value="genomeCoverageBed1.len" ftype="tabular" />
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
103 <param name="report_select" value="hist" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
104 <output name="output" file="genomeCoverageBed_result1.bed" ftype="tabular" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
105 </test>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
106 </tests>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
107 <help>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
108 <![CDATA[
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
109 **What it does**
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
110
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
111 This tool calculates the genome-wide coverage of intervals defined in a BAM or BED file and reports them in BedGraph format.
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
112
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
113 .. image:: $PATH_TO_IMAGES/genomecov-glyph.png
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
114
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
115 .. class:: warningmark
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
116
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
117 The input BED or BAM file must be sorted by chromosome name (but doesn't necessarily have to be sorted by start position).
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
118
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
119 -----
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
120
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
121 **Example 1**
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
122
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
123 Input (BED format)-
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
124 Overlapping, un-sorted intervals::
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
125
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
126 chr1 140 176
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
127 chr1 100 130
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
128 chr1 120 147
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
129
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
130
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
131 Output (BedGraph format)-
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
132 Sorted, non-overlapping intervals, with coverage value on the 4th column::
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
133
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
134 chr1 100 120 1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
135 chr1 120 130 2
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
136 chr1 130 140 1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
137 chr1 140 147 2
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
138 chr1 147 176 1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
139
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
140 -----
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
141
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
142 **Example 2 - with ZERO-Regions selected (assuming hg19)**
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
143
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
144 Input (BED format)-
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
145 Overlapping, un-sorted intervals::
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
146
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
147 chr1 140 176
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
148 chr1 100 130
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
149 chr1 120 147
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
150
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
151
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
152 BedGraph output will contain five columns:
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
153
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
154 * 1. Chromosome name (or 'genome' for whole-genome coverage)
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
155 * 2. Coverage depth
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
156 * 3. The number of bases on chromosome (or genome) with depth equal to column 2.
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
157 * 4. The size of chromosome (or entire genome) in base pairs
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
158 * 5. The fraction of bases on chromosome (or entire genome) with depth equal to column 2.
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
159
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
160 **Example Output**:
11
7308cc546a36 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 19967671ddd4e750f53d71274e5c12325a402811
iuc
parents: 7
diff changeset
161
1
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
162 chr2L 0 1379895 23011544 0.0599653
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
163 chr2L 1 837250 23011544 0.0363839
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
164 chr2L 2 904442 23011544 0.0393038
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
165 chr2L 3 913723 23011544 0.0397072
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
166 chr2L 4 952166 23011544 0.0413778
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
167 chr2L 5 967763 23011544 0.0420555
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
168 chr2L 6 986331 23011544 0.0428624
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
169 chr2L 7 998244 23011544 0.0433801
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
170 chr2L 8 995791 23011544 0.0432735
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
171 chr2L 9 996398 23011544 0.0432999
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
172
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
173
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
174 @REFERENCES@
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
175 ]]>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
176 </help>
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
177 <expand macro="citations" />
82aac94b06c3 Uploaded
iuc
parents:
diff changeset
178 </tool>