Mercurial > repos > iuc > bedtools
comparison mapBed.xml @ 1:82aac94b06c3 draft
Uploaded
| author | iuc |
|---|---|
| date | Thu, 08 Jan 2015 14:25:51 -0500 |
| parents | b8348686a0b9 |
| children | 457b09031d57 |
comparison
equal
deleted
inserted
replaced
| 0:b8348686a0b9 | 1:82aac94b06c3 |
|---|---|
| 4 <import>macros.xml</import> | 4 <import>macros.xml</import> |
| 5 </macros> | 5 </macros> |
| 6 <expand macro="requirements" /> | 6 <expand macro="requirements" /> |
| 7 <expand macro="stdio" /> | 7 <expand macro="stdio" /> |
| 8 <command> | 8 <command> |
| 9 <![CDATA[ | |
| 9 bedtools map | 10 bedtools map |
| 10 -a $inputA | 11 -a "${inputA}" |
| 11 -b $inputB | 12 -b "${inputB}" |
| 12 $strand | 13 $strand |
| 13 -o $operation | 14 -o $operation |
| 14 -c $col | 15 -c "${cols}" |
| 15 -f $overlap | 16 -f $overlap |
| 16 $reciprocal | 17 $reciprocal |
| 17 $split | 18 $split |
| 18 $header | 19 $header |
| 19 #if $genome.genome_choose == "-g" : | 20 #if $genome.genome_choose == "-g" : |
| 20 -g $genome.genome | 21 -g $genome.genome |
| 21 #end if | 22 #end if |
| 22 > $output | 23 > "${output}" |
| 24 ]]> | |
| 23 </command> | 25 </command> |
| 24 <inputs> | 26 <inputs> |
| 25 <param format="bed,vcf,gff,gff3" name="inputA" type="data" label="BED/VCF/GFF file A"/> | 27 <param format="bam,bed,vcf,gff,gff3" name="inputA" type="data" label="File A (BAM/BED/VCF/GFF)" /> |
| 26 <param format="bed,gff,vcf,gff3" name="inputB" type="data" label="BED/VCF/GFF file B"/> | 28 <param format="bam,bed,gff,vcf,gff3" name="inputB" type="data" label="File B (BAM/BED/VCF/GFF)" /> |
| 27 <param name="col" type="integer" value="5" label="Specify the column from the B file to map onto intervals in A" /> | 29 <expand macro="choose_columns" /> |
| 28 <expand macro="overlap" /> | 30 <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." /> | 31 <param name="reciprocal" type="boolean" checked="false" truevalue="-r" falsevalue="" |
| 32 label="Require reciprocal overlap" | |
| 33 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. (-r)" /> | |
| 30 <expand macro="strand2" /> | 34 <expand macro="strand2" /> |
| 31 <param name="operation" type="select" label="Specify the operation"> | 35 <expand macro="choose_operations"> |
| 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="math_options" /> |
| 36 <expand macro="additional_math_options" /> | 37 <expand macro="additional_math_options" /> |
| 37 </param> | 38 </expand> |
| 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 <expand macro="split" /> |
| 39 <param name="header" type="boolean" checked="false" truevalue="-header" falsevalue="" label="Print the header from the A file prior to results." /> | 40 <expand macro="print_header" /> |
| 40 <conditional name="genome"> | 41 <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 <param name="genome_choose" type="boolean" checked="false" truevalue="-g" falsevalue="" |
| 43 label="Treat split/spliced BAM or BED12 entries as distinct BED intervals when computing coverage." help="(-g)" /> | |
| 42 <when value="-g"> | 44 <when value="-g"> |
| 43 <expand macro="genome" /> | 45 <expand macro="genome" /> |
| 44 </when> | 46 </when> |
| 45 </conditional> | 47 </conditional> |
| 46 </inputs> | 48 </inputs> |
| 47 <outputs> | 49 <outputs> |
| 48 <data format_source="inputA" name="output" metadata_source="inputA" label="Mapping of ${inputB.name} into ${inputA.name}"/> | 50 <data format_source="inputA" name="output" metadata_source="inputA" label="Mapping of ${inputB.name} into ${inputA.name}"/> |
| 49 </outputs> | 51 </outputs> |
| 52 <tests> | |
| 53 <test> | |
| 54 <param name="inputA" value="mapBed1.bed" ftype="bed" /> | |
| 55 <param name="inputB" value="mapBed2.bed" ftype="bed" /> | |
| 56 <param name="cols" value="5" /> | |
| 57 <param name="operation" value="mean" /> | |
| 58 <output name="output" file="mapBed_result1.bed" ftype="bed" /> | |
| 59 </test> | |
| 60 <test> | |
| 61 <param name="inputA" value="mapBed1.bed" ftype="bed" /> | |
| 62 <param name="inputB" value="mapBed2.bed" ftype="bed" /> | |
| 63 <param name="cols" value="5" /> | |
| 64 <param name="operation" value="collapse" /> | |
| 65 <output name="output" file="mapBed_result2.bed" ftype="bed" /> | |
| 66 </test> | |
| 67 <test> | |
| 68 <param name="inputA" value="mapBed1.bed" ftype="bed" /> | |
| 69 <param name="inputB" value="mapBed2.bed" ftype="bed" /> | |
| 70 <param name="cols" value="5" /> | |
| 71 <param name="operation" value="collapse" /> | |
| 72 <param name="strand" value="-S" /> | |
| 73 <output name="output" file="mapBed_result3.bed" ftype="bed" /> | |
| 74 </test> | |
| 75 <test> | |
| 76 <param name="inputA" value="mapBed1.bed" ftype="bed" /> | |
| 77 <param name="inputB" value="mapBed2.bed" ftype="bed" /> | |
| 78 <param name="cols" value="5" /> | |
| 79 <param name="operation" value="collapse" /> | |
| 80 <param name="strand" value="-s" /> | |
| 81 <output name="output" file="mapBed_result4.bed" ftype="bed" /> | |
| 82 </test> | |
| 83 </tests> | |
| 50 <help> | 84 <help> |
| 51 | 85 <![CDATA[ |
| 52 **What it does** | 86 **What it does** |
| 53 | 87 |
| 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. | 88 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 | 89 |
| 56 .. image:: $PATH_TO_IMAGES/map-glyph.png | 90 .. image:: $PATH_TO_IMAGES/map-glyph.png |
| 62 .. class:: infomark | 96 .. class:: infomark |
| 63 | 97 |
| 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. | 98 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 | 99 |
| 66 @REFERENCES@ | 100 @REFERENCES@ |
| 67 | 101 ]]> |
| 68 </help> | 102 </help> |
| 69 <expand macro="citations" /> | 103 <expand macro="citations" /> |
| 70 </tool> | 104 </tool> |