|
11
|
1 <tool id="bedtools_map" name="MapBed" version="0.2.0">
|
|
|
2 <description></description>
|
|
|
3 <macros>
|
|
|
4 <import>macros.xml</import>
|
|
|
5 </macros>
|
|
|
6 <expand macro="requirements" />
|
|
|
7 <expand macro="stdio" />
|
|
|
8 <command>
|
|
|
9 bedtools map
|
|
|
10 -a $inputA
|
|
|
11 -b $inputB
|
|
|
12 $strand
|
|
|
13 -o $operation
|
|
|
14 -c $col
|
|
|
15 -f $overlap
|
|
|
16 $reciprocal
|
|
|
17 $split
|
|
|
18 $header
|
|
|
19 #if $genome.genome_choose == "-g" :
|
|
|
20 -g $genome.genome
|
|
|
21 #end if
|
|
|
22 > $output
|
|
|
23 </command>
|
|
|
24 <inputs>
|
|
|
25 <param format="bed,vcf,gff,gff3" name="inputA" type="data" label="BED/VCF/GFF file A"/>
|
|
|
26 <param format="bed,gff,vcf,gff3" name="inputB" type="data" label="BED/VCF/GFF file B"/>
|
|
|
27 <param name="col" type="integer" value="5" label="Specify the column from the B file to map onto intervals in A" />
|
|
|
28 <expand macro="overlap" />
|
|
|
29 <param name="reciprocal" type="boolean" checked="false" truevalue="-r" falsevalue="" label="Require reciprocal overlap." help="If set, the overlap between the BAM alignment and the BED interval must affect the above fraction of both the alignment and the BED interval." />
|
|
|
30 <expand macro="strand2" />
|
|
|
31 <param name="operation" type="select" label="Specify the operation">
|
|
|
32 <option value="sum">Sum - numeric only</option>
|
|
|
33 <option value="absmin">AbsMin - numeric only</option>
|
|
|
34 <option value="absmax">AbsMax - numeric only</option>
|
|
|
35 <expand macro="math_options" />
|
|
|
36 <expand macro="additional_math_options" />
|
|
|
37 </param>
|
|
|
38 <param name="split" type="boolean" checked="true" 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." />
|
|
|
39 <param name="header" type="boolean" checked="false" truevalue="-header" falsevalue="" label="Print the header from the A file prior to results." />
|
|
|
40 <conditional name="genome">
|
|
|
41 <param name="genome_choose" type="boolean" checked="false" truevalue="-g" falsevalue="" label="Treat split/spliced BAM or BED12 entries as distinct BED intervals when computing coverage." help="" />
|
|
|
42 <when value="-g">
|
|
|
43 <expand macro="genome" />
|
|
|
44 </when>
|
|
|
45 </conditional>
|
|
|
46 </inputs>
|
|
|
47 <outputs>
|
|
|
48 <data format_source="inputA" name="output" metadata_source="inputA" label="Mapping of ${inputB.name} into ${inputA.name}"/>
|
|
|
49 </outputs>
|
|
|
50 <help>
|
|
|
51
|
|
|
52 **What it does**
|
|
|
53
|
|
|
54 bedtools map allows one to map overlapping features in a B file onto features in an A file and apply statistics and/or summary operations on those features.
|
|
|
55
|
|
|
56 .. image:: $PATH_TO_IMAGES/map-glyph.png
|
|
|
57
|
|
|
58 .. class:: infomark
|
|
|
59
|
|
|
60 bedtools map requires each input file to be sorted by genome coordinate. For BED files, this can be done with sort -k1,1 -k2,2n. Other sorting criteria are allowed if a genome file (-g) is provides that specifies the expected chromosome order.
|
|
|
61
|
|
|
62 .. class:: infomark
|
|
|
63
|
|
|
64 The map tool is substantially faster in versions 2.19.0 and later. The plot below demonstrates the increased speed when, for example, counting the number of exome alignments that align to each exon. The bedtools times are compared to the bedops bedmap utility as a point of reference.
|
|
|
65
|
|
|
66 @REFERENCES@
|
|
|
67
|
|
|
68 </help>
|
|
|
69 </tool>
|
