Mercurial > repos > devteam > coverage

comparison coverage.xml @ 0:1e864693a1c0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Imported from capsule None
author devteam
date Tue, 01 Apr 2014 10:53:30 -0400
parents
children 16b62a897aac
comparison
equal deleted inserted replaced
-1:000000000000 0:1e864693a1c0
1 <tool id="gops_coverage_1" name="Coverage" version="0.0.1">
2 <description>of a set of intervals on second set of intervals</description>
3 <requirements>
4 <requirement type="package" version="0.7.1">bx-python</requirement>
5 <requirement type="package" version="1.0.0">galaxy-ops</requirement>
6 </requirements>
7 <command interpreter="python">gops_coverage.py $input1 $input2 $output -1 ${input1.metadata.chromCol},${input1.metadata.startCol},${input1.metadata.endCol},${input1.metadata.strandCol} -2 ${input2.metadata.chromCol},${input2.metadata.startCol},${input2.metadata.endCol},${input2.metadata.strandCol}</command>
8 <inputs>
9 <param format="interval" name="input1" type="data" help="First dataset">
10 <label>What portion of</label>
11 </param>
12 <param format="interval" name="input2" type="data" help="Second dataset">
13 <label>is covered by</label>
14 </param>
15 </inputs>
16 <outputs>
17 <data format="interval" name="output" metadata_source="input1" />
18 </outputs>
19 <code file="operation_filter.py"/>
20 <tests>
21 <test>
22 <param name="input1" value="1.bed" />
23 <param name="input2" value="2.bed" />
24 <output name="output" file="gops_coverage_out.interval" />
25 </test>
26 <test>
27 <param name="input1" value="1.bed" />
28 <param name="input2" value="2_mod.bed" ftype="interval"/>
29 <output name="output" file="gops_coverage_out_diffCols.interval" />
30 </test>
31 <test>
32 <param name="input1" value="gops_bigint.interval" />
33 <param name="input2" value="gops_bigint2.interval" />
34 <output name="output" file="gops_coverage_out2.interval" />
35 </test>
36 </tests>
37 <help>
38
39 .. class:: infomark
40
41 **TIP:** If your dataset does not appear in the pulldown menu -> it is not in interval format. Use "edit attributes" to set chromosome, start, end, and strand columns.
42
43 Find the coverage of intervals in the first dataset on intervals in the second dataset. The coverage is added as two columns, the first being bases covered, and the second being the fraction of bases covered by that interval.
44
45 -----
46
47 **Screencasts!**
48
49 See Galaxy Interval Operation Screencasts_ (right click to open this link in another window).
50
51 .. _Screencasts: http://wiki.g2.bx.psu.edu/Learn/Interval%20Operations
52
53 -----
54
55 **Example**
56
57
58 if **First dataset** are genes ::
59
60 chr11 5203271 5204877 NM_000518 0 -
61 chr11 5210634 5212434 NM_000519 0 -
62 chr11 5226077 5227663 NM_000559 0 -
63 chr11 5226079 5232587 BC020719 0 -
64 chr11 5230996 5232587 NM_000184 0 -
65
66 and **Second dataset** are repeats::
67
68 chr11 5203895 5203991 L1MA6 500 +
69 chr11 5204163 5204239 A-rich 219 +
70 chr11 5211034 5211167 (CATATA)n 245 +
71 chr11 5211642 5211673 AT_rich 24 +
72 chr11 5226551 5226606 (CA)n 303 +
73 chr11 5228782 5228825 (TTTTTG)n 208 +
74 chr11 5229045 5229121 L1PA11 440 +
75 chr11 5229133 5229319 MER41A 1106 +
76 chr11 5229374 5229485 L2 244 -
77 chr11 5229751 5230083 MLT1A 913 -
78 chr11 5231469 5231526 (CA)n 330 +
79
80 the Result is the coverage density of repeats in the genes::
81
82 chr11 5203271 5204877 NM_000518 0 - 172 0.107098
83 chr11 5210634 5212434 NM_000519 0 - 164 0.091111
84 chr11 5226077 5227663 NM_000559 0 - 55 0.034678
85 chr11 5226079 5232587 BC020719 0 - 860 0.132145
86 chr11 5230996 5232587 NM_000184 0 - 57 0.035827
87
88 For example, the following line of output::
89
90 chr11 5203271 5204877 NM_000518 0 - 172 0.107098
91
92 implies that 172 nucleotides accounting for 10.7% of the this interval (chr11:5203271-5204877) overlap with repetitive elements.
93
94 </help>
95 </tool>