annotate tools/new_operations/coverage.xml @ 0:9071e359b9a3

Uploaded
author xuebing
date Fri, 09 Mar 2012 19:37:19 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
1 <tool id="gops_coverage_1" name="Coverage">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
2 <description>of a set of intervals on second set of intervals</description>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
3 <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>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
4 <inputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
5 <param format="interval" name="input1" type="data" help="First dataset">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
6 <label>What portion of</label>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
7 </param>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
8 <param format="interval" name="input2" type="data" help="Second dataset">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
9 <label>is covered by</label>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
10 </param>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
11 </inputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
12 <outputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
13 <data format="interval" name="output" metadata_source="input1" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
14 </outputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
15 <code file="operation_filter.py"/>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
16 <tests>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
17 <test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
18 <param name="input1" value="1.bed" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
19 <param name="input2" value="2.bed" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
20 <output name="output" file="gops_coverage_out.interval" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
21 </test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
22 <test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
23 <param name="input1" value="1.bed" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
24 <param name="input2" value="2_mod.bed" ftype="interval"/>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
25 <output name="output" file="gops_coverage_out_diffCols.interval" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
26 </test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
27 <test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
28 <param name="input1" value="gops_bigint.interval" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
29 <param name="input2" value="gops_bigint2.interval" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
30 <output name="output" file="gops_coverage_out2.interval" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
31 </test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
32 </tests>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
33 <help>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
34
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
35 .. class:: infomark
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
36
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
37 **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.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
38
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
39 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.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
40
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
41 -----
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
42
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
43 **Screencasts!**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
44
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
45 See Galaxy Interval Operation Screencasts_ (right click to open this link in another window).
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
46
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
47 .. _Screencasts: http://wiki.g2.bx.psu.edu/Learn/Interval%20Operations
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
48
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
49 -----
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
50
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
51 **Example**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
52
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
53
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
54 if **First dataset** are genes ::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
55
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
56 chr11 5203271 5204877 NM_000518 0 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
57 chr11 5210634 5212434 NM_000519 0 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
58 chr11 5226077 5227663 NM_000559 0 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
59 chr11 5226079 5232587 BC020719 0 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
60 chr11 5230996 5232587 NM_000184 0 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
61
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
62 and **Second dataset** are repeats::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
63
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
64 chr11 5203895 5203991 L1MA6 500 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
65 chr11 5204163 5204239 A-rich 219 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
66 chr11 5211034 5211167 (CATATA)n 245 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
67 chr11 5211642 5211673 AT_rich 24 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
68 chr11 5226551 5226606 (CA)n 303 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
69 chr11 5228782 5228825 (TTTTTG)n 208 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
70 chr11 5229045 5229121 L1PA11 440 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
71 chr11 5229133 5229319 MER41A 1106 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
72 chr11 5229374 5229485 L2 244 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
73 chr11 5229751 5230083 MLT1A 913 -
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
74 chr11 5231469 5231526 (CA)n 330 +
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
75
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
76 the Result is the coverage density of repeats in the genes::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
77
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
78 chr11 5203271 5204877 NM_000518 0 - 172 0.107098
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
79 chr11 5210634 5212434 NM_000519 0 - 164 0.091111
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
80 chr11 5226077 5227663 NM_000559 0 - 55 0.034678
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
81 chr11 5226079 5232587 BC020719 0 - 860 0.132145
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
82 chr11 5230996 5232587 NM_000184 0 - 57 0.035827
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
83
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
84 For example, the following line of output::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
85
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
86 chr11 5203271 5204877 NM_000518 0 - 172 0.107098
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
87
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
88 implies that 172 nucleotides accounting for 10.7% of the this interval (chr11:5203271-5204877) overlap with repetitive elements.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
89
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
90 </help>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
91 </tool>