Mercurial > repos > iuc > bedtools
comparison intersectBed.xml @ 20:df56e1b12d0c draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 9f7b029951fa0f854c9c6a5f241ce9a20ae4f22a
author | iuc |
---|---|
date | Wed, 10 Jan 2018 19:21:33 -0500 |
parents | fadebae7e69b |
children | 95a3b2c25bd1 |
comparison
equal
deleted
inserted
replaced
19:a8eabd2838f6 | 20:df56e1b12d0c |
---|---|
1 <tool id="bedtools_intersectbed" name="Intersect intervals" version="@WRAPPER_VERSION@.0"> | 1 <tool id="bedtools_intersectbed" name="Intersect intervals" version="@WRAPPER_VERSION@.1"> |
2 <description>find overlapping intervals in various ways</description> | 2 <description>find overlapping intervals in various ways</description> |
3 <macros> | 3 <macros> |
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 <![CDATA[ |
10 #set inputBs = '" "'.join( [ str( $file ) for $file in $inputB ] ) | |
11 #set modes = ' '.join( str($overlap_mode).split(',') ) | 10 #set modes = ' '.join( str($overlap_mode).split(',') ) |
12 | 11 |
13 bedtools intersect | 12 bedtools intersect |
14 #if $inputA.is_of_type('bam'): | 13 #if $inputA.is_of_type('bam'): |
15 -abam "${inputA}" | 14 -abam '${inputA}' |
16 #else: | 15 #else: |
17 -a "${inputA}" | 16 -a '${inputA}' |
18 #end if | 17 #end if |
19 | 18 |
20 -b "${inputBs}" | 19 #if str($reduce_or_iterate.reduce_or_iterate_selector) == 'iterate': |
20 -b '$reduce_or_iterate.inputB' | |
21 #else: | |
22 #set inputBs = "' '".join([str($file) for $file in $reduce_or_iterate.inputB]) | |
23 -b '$inputBs' | |
24 #end if | |
21 $split | 25 $split |
22 $strand | 26 $strand |
23 #if str($fraction) != "None" and str($fraction): | 27 #if str($fraction) != "None" and str($fraction): |
24 -f "${fraction}" | 28 -f '${fraction}' |
25 #end if | 29 #end if |
26 $reciprocal | 30 $reciprocal |
27 $invert | 31 $invert |
28 $once | 32 $once |
29 $header | 33 $header |
30 $modes | 34 $modes |
31 $count | 35 $count |
32 > "${output}" | 36 > '${output}' |
33 ]]> | 37 ]]> |
34 </command> | 38 </command> |
35 <inputs> | 39 <inputs> |
36 <param format="bed,bam,vcf,gff,gff3" name="inputA" type="data" label="BED/VCF/GFF/BAM file"/> | 40 <param format="bed,bam,vcf,gff,gff3" name="inputA" type="data" label="File A to intersect with B" help="BAM/BED/GFF/VCF format" /> |
37 <param format="bed,bam,gff,vcf,gff3" name="inputB" type="data" multiple="True" label="One or more BAM/BED/GFF/VCF file(s)"/> | 41 <conditional name="reduce_or_iterate"> |
42 <param name='reduce_or_iterate_selector' type='select' label='Combined or separate output files'> | |
43 <option value='iterate' selected='true'>One output file per 'input B' file</option> | |
44 <option value='reduce'>Single output containing intersections of any 'input B' lines with A </option> | |
45 </param> | |
46 <when value='iterate'> | |
47 <param format="bam,bed,gff,gff3,vcf" name="inputB" type="data" | |
48 label="File(s) B to intersect with A" help="BAM/BED/GFF/VCF format"/> | |
49 </when> | |
50 <when value='reduce'> | |
51 <param format="bam,bed,gff,gff3,vcf" name="inputB" type="data" multiple="true" | |
52 label="File(s) B to intersect with A" help="BAM/BED/GFF/VCF format"/> | |
53 </when> | |
54 </conditional> | |
38 <expand macro="strand2" /> | 55 <expand macro="strand2" /> |
39 <param name="overlap_mode" type="select" multiple="True" label="What should be written to the output file?"> | 56 <param name="overlap_mode" type="select" multiple="True" label="What should be written to the output file?"> |
40 <option value="-wa" selected="True">Write the original entry in A for each overlap (-wa)</option> | 57 <option value="-wa" selected="True">Write the original entry in A for each overlap (-wa)</option> |
41 <option value="-wb">Write the original entry in B for each overlap. Useful for knowing what A overlaps. Restricted by the fraction- and reciprocal option (-wb)</option> | 58 <option value="-wb">Write the original entry in B for each overlap. Useful for knowing what A overlaps. Restricted by the fraction- and reciprocal option (-wb)</option> |
42 <option value="-wo">Write the original A and B entries plus the number of base pairs of overlap between the two features. Only A features with overlap are reported. Restricted by the fraction- and reciprocal option (-wo)</option> | 59 <option value="-wo">Write the original A and B entries plus the number of base pairs of overlap between the two features. Only A features with overlap are reported. Restricted by the fraction- and reciprocal option (-wo)</option> |
51 help="Alignments are only retained if the overlap with the an interval in the BED file comprises at least this fraction of the BAM alignment's length. For example, to require that the overlap affects 50% of the BAM alignment, use 0.50. (-f)"/> | 68 help="Alignments are only retained if the overlap with the an interval in the BED file comprises at least this fraction of the BAM alignment's length. For example, to require that the overlap affects 50% of the BAM alignment, use 0.50. (-f)"/> |
52 <!-- -r --> | 69 <!-- -r --> |
53 <expand macro="reciprocal" /> | 70 <expand macro="reciprocal" /> |
54 <!-- -v --> | 71 <!-- -v --> |
55 <param name="invert" type="boolean" checked="false" truevalue="-v" falsevalue="" | 72 <param name="invert" type="boolean" checked="false" truevalue="-v" falsevalue="" |
56 label="Report only those alignments that **do not** overlap the BED file" | 73 label="Report only those alignments that **do not** overlap with file(s) B" |
57 help="(-v)"/> | 74 help="(-v)"/> |
58 <!-- -u --> | 75 <!-- -u --> |
59 <param name="once" type="boolean" checked="false" truevalue="-u" falsevalue="" | 76 <param name="once" type="boolean" checked="false" truevalue="-u" falsevalue="" |
60 label="Write the original A entry _once_ if _any_ overlaps found in B." | 77 label="Write the original A entry _once_ if _any_ overlaps found in B." |
61 help="Just report the fact >=1 hit was found. (-u)" /> | 78 help="Just report the fact >=1 hit was found. (-u)" /> |