Mercurial > repos > iuc > bedtools

comparison mapBed.xml @ 0:b8348686a0b9 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Imported from capsule None
author iuc
date Tue, 04 Nov 2014 01:45:04 -0500
parents
children 82aac94b06c3
comparison
equal deleted inserted replaced
-1:000000000000 0:b8348686a0b9
1 <tool id="bedtools_map" name="MapBed" version="@WRAPPER_VERSION@.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 &gt; $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 <expand macro="citations" />
70 </tool>