|
10
|
1 <tool id="bedtools_coveragebed_counts" name="Count intervals in one file overlapping intervals in another file." version="0.1.0">
|
|
|
2
|
|
|
3 <requirements>
|
|
|
4 <requirement type="binary">coverageBed</requirement>
|
|
|
5 </requirements>
|
|
|
6
|
|
|
7 <command>
|
|
|
8 coverageBed
|
|
|
9 #if $inputA.ext == "bam"
|
|
|
10 -abam '$inputA'
|
|
|
11 #else
|
|
|
12 -a '$inputA'
|
|
|
13 #end if
|
|
|
14 -b '$inputB'
|
|
|
15 -counts
|
|
|
16 $split
|
|
|
17 $strand
|
|
|
18 | sort -k1,1 -k2,2n
|
|
|
19 > '$output'
|
|
|
20 </command>
|
|
|
21
|
|
|
22 <inputs>
|
|
|
23 <param format="bed,bam" name="inputA" type="data" label="Count how many intervals in this BED or BAM file (source)">
|
|
|
24 <validator type="unspecified_build" />
|
|
|
25 </param>
|
|
|
26 <param format="bed" name="inputB" type="data" label="overlap the intervals in this BED file (target)">
|
|
|
27 <validator type="unspecified_build" />
|
|
|
28 </param>
|
|
|
29 <param name="split" type="boolean" checked="false" truevalue="-split" falsevalue="" label="Treat split/spliced BAM or BED12 entries as distinct BED intervals when computing coverage." help="If set, the coverage will be calculated based the spliced intervals only. For BAM files, this inspects the CIGAR N operation to infer the blocks for computing coverage. For BED12 files, this inspects the BlockCount, BlockStarts, and BlockEnds fields (i.e., columns 10,11,12). If this option is not set, coverage will be calculated based on the interval's START/END coordinates, and would include introns in the case of RNAseq data." />
|
|
|
30
|
|
|
31 <param name="strand" type="select" label="Count">
|
|
|
32 <option value="">overlaps on either strand</option>
|
|
|
33 <option value="-s">only overlaps occurring on the **same** strand.</option>
|
|
|
34 <option value="-S">only overlaps occurring on the **opposite** strand.</option>
|
|
|
35 </param>
|
|
|
36 </inputs>
|
|
|
37
|
|
|
38 <outputs>
|
|
|
39 <data format="bed" name="output" metadata_source="inputB" label="count of overlaps in ${inputA.name} on ${inputB.name}"/>
|
|
|
40 </outputs>
|
|
|
41
|
|
|
42 <help>
|
|
|
43
|
|
|
44 **What it does**
|
|
|
45
|
|
|
46 This tool converts counts the number of intervals in a BAM or BED file (the source) that overlap another BED file (the target).
|
|
|
47
|
|
|
48 .. class:: infomark
|
|
|
49
|
|
|
50 The output file will be comprised of each interval from your original target BED file, plus an additional column indicating the number of intervals in your source file that overlapped that target interval.
|
|
|
51
|
|
|
52 .. class:: warningmark
|
|
|
53
|
|
|
54 This tool requires that `bedtools`__ has been installed on your system.
|
|
|
55
|
|
|
56 ------
|
|
|
57
|
|
|
58 This tool is part of the `bedtools package`__ from the `Quinlan laboratory`__. If you use this tool, please cite `Quinlan AR, and Hall I.M. BEDTools: A flexible framework for comparing genomic features. Bioinformatics, 2010, 26, 6.`__
|
|
|
59
|
|
|
60 .. __: http://code.google.com/p/bedtools/
|
|
|
61 .. __: http://code.google.com/p/bedtools/
|
|
|
62 .. __: http://cphg.virginia.edu/quinlan/
|
|
|
63 .. __: http://bioinformatics.oxfordjournals.org/content/26/6/841.short
|
|
|
64
|
|
|
65 </help>
|
|
|
66 </tool>
|
