# HG changeset patch # User iuc # Date 1504619516 14400 # Node ID b0d5e752c0c5c5daefb7c32dda39743c712a0ac4 # Parent cb98f29aa1dc011b37a9c629ad7bfe36d079caf0 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit a78ca56350b1fcf9d7de2a1112155c9bd0844348 diff -r cb98f29aa1dc -r b0d5e752c0c5 bedToBam.xml --- a/bedToBam.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/bedToBam.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + converter macros.xml @@ -10,7 +10,7 @@ bedtools bedtobam $bed12 -mapq $mapq - -g $genome + -g @GENOME_FILE@ -i '$input' > '$output' ]]> @@ -20,7 +20,7 @@ - + @@ -30,6 +30,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 bedpeToBam.xml --- a/bedpeToBam.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/bedpeToBam.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + converter macros.xml @@ -10,13 +10,13 @@ bedtools bedpetobam -mapq $mapq -i '$input' - -g $genome + -g @GENOME_FILE@ > '$output' ]]> - + @@ -27,6 +27,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 complementBed.xml --- a/complementBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/complementBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + Extract intervals not represented by an interval file macros.xml @@ -9,13 +9,13 @@ "$output" ]]> - + @@ -23,6 +23,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 fisherBed.xml --- a/fisherBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/fisherBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + calculate Fisher statistic between two feature files macros.xml @@ -13,7 +13,7 @@ -a '$inputA' -b '$inputB' -f $overlap - -g $genome + -g @GENOME_FILE@ $reciprocal $m > '$output' @@ -22,7 +22,7 @@ - + @@ -37,6 +37,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 flankBed.xml --- a/flankBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/flankBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + create new intervals from the flanks of existing intervals macros.xml @@ -10,7 +10,7 @@ flankBed $pct $strand - -g $genome + -g @GENOME_FILE@ -i '$input' #if $addition.addition_select == 'b': @@ -24,7 +24,7 @@ - + @@ -39,6 +39,7 @@ + @@ -46,6 +47,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 genomeCoverageBed.xml --- a/genomeCoverageBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/genomeCoverageBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + compute the coverage over an entire genome macros.xml @@ -8,12 +8,7 @@ - + @@ -65,11 +60,11 @@ label="Report regions with zero coverage" help="If set, regions without any coverage will also be reported. (-bga)" /> + help="Each bedGraph coverage value is multiplied by this factor before being reported. Useful for normalizing coverage by, e.g., reads per million (RPM). (-scale)" /> + help="Combine all positions with a depth >= max into a single bin in the histogram. (-max)" /> @@ -99,6 +94,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 macros.xml --- a/macros.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/macros.xml Tue Sep 05 09:51:56 2017 -0400 @@ -49,10 +49,68 @@ 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. (-split)" /> - - - + + + + + + + + + + + + + + + + +#if $genome_file_opts.genome_file_opts_selector == "loc": + '$genome_file_opts.genome.fields.len_path' +#elif $genome_file_opts.genome_file_opts_selector == "hist": + '$genome_file_opts.genome' +#end if + + +#if $genome.genome_choose == "-g": + #if $genome.genome_file_opts.genome_file_opts_selector == "loc": + -g '$genome.genome_file_opts.genome.fields.len_path' + #elif $genome.genome_file_opts.genome_file_opts_selector == "hist": + -g '$genome.genome_file_opts.genome' + #end if +#end if + + +#if $type.type_select == "genome": + #if $type.genome_file_opts.genome_file_opts_selector == "loc": + -g '$type.genome_file_opts.genome.fields.len_path' + #elif $type.genome_file_opts.genome_file_opts_selector == "hist": + -g '$type.genome_file_opts.genome' + #end if +#end if + + +#if $empty.empty_selector == "-empty": + #if $empty.genome_file_opts.genome_file_opts_selector == "loc": + -g '$empty.genome_file_opts.genome.fields.len_path' + #elif $empty.genome_file_opts.genome_file_opts_selector == "hist": + -g '$empty.genome_file_opts.genome' + #end if +#end if + + +#if $input_type.input_type_select == "bed": + -i '$input_type.input' + #if $input_type.genome_file_opts.genome_file_opts_selector == "loc": + -g '$input_type.genome_file_opts.genome.fields.len_path' + #elif $input_type.genome_file_opts.genome_file_opts_selector == "hist": + -g '$input_type.genome_file_opts.genome' + #end if +#elif $input_type.input_type_select == "bam": + -ibam '$input_type.input' +#end if + + make interval windows across a genome macros.xml @@ -9,7 +9,7 @@ - + @@ -81,6 +81,7 @@ + @@ -88,6 +89,7 @@ + @@ -97,6 +99,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 mapBed.xml --- a/mapBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/mapBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + apply a function to a column for each overlapping interval macros.xml @@ -8,17 +8,15 @@ "${output}" ]]> @@ -36,16 +34,19 @@ - + + + + - + - + @@ -86,6 +87,19 @@ + + + + + + + + + + + + + + identifies common intervals among multiple interval files macros.xml @@ -11,10 +11,8 @@ $header $cluster -filler "${filler}" - #if $zero.value == True: - -empty - -g $genome - #end if + $empty.empty_selector + @GENOME_FILE_UNION@ #if str($tag.tag_select) == "tag": #set files = '" "'.join( [ str( $file ) for $file in $tag.inputs ] ) @@ -45,17 +43,23 @@ - + - - + + + + + + + + + + @@ -67,23 +71,24 @@ - + - + - + - - + + - + - - + + + diff -r cb98f29aa1dc -r b0d5e752c0c5 randomBed.xml --- a/randomBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/randomBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + generate random intervals in a genome macros.xml @@ -8,7 +8,7 @@ - + @@ -28,6 +28,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 shuffleBed.xml --- a/shuffleBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/shuffleBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + randomly redistrubute intervals in a genome macros.xml @@ -8,7 +8,7 @@ - + @@ -67,6 +67,7 @@ + @@ -75,6 +76,7 @@ + @@ -83,6 +85,7 @@ + @@ -91,6 +94,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 slopBed.xml --- a/slopBed.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/slopBed.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + adjust the size of intervals macros.xml @@ -10,7 +10,7 @@ bedtools slop $pct $strand - -g '$genome' + -g @GENOME_FILE@ -i '$inputA' #if $addition.addition_select == 'b': -b $addition.b @@ -24,7 +24,7 @@ - + @@ -40,6 +40,7 @@ + @@ -47,6 +48,7 @@ + diff -r cb98f29aa1dc -r b0d5e752c0c5 test-data/mapBed3.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBed3.bed Tue Sep 05 09:51:56 2017 -0400 @@ -0,0 +1,6 @@ +chr1 10 20 a1 1 + +chr1 50 60 a2 2 - +chr1 80 90 a3 3 - +chr9 10 20 a1 1 + +chr9 50 60 a2 2 - +chr9 80 90 a3 3 - diff -r cb98f29aa1dc -r b0d5e752c0c5 test-data/mapBed4.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBed4.bed Tue Sep 05 09:51:56 2017 -0400 @@ -0,0 +1,10 @@ +chr1 12 14 b1 2 + +chr1 13 15 b2 5 - +chr1 16 18 b3 5 + +chr1 82 85 b4 2 - +chr1 85 87 b5 3 + +chr9 12 14 b1 2 + +chr9 13 15 b2 5 - +chr9 16 18 b3 5 + +chr9 82 85 b4 2 - +chr9 85 87 b5 3 + diff -r cb98f29aa1dc -r b0d5e752c0c5 test-data/mapBed_result5.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBed_result5.bed Tue Sep 05 09:51:56 2017 -0400 @@ -0,0 +1,6 @@ +chr1 10 20 a1 1 + 2,5 +chr1 50 60 a2 2 - . +chr1 80 90 a3 3 - 2 +chr9 10 20 a1 1 + 2,5 +chr9 50 60 a2 2 - . +chr9 80 90 a3 3 - 2 diff -r cb98f29aa1dc -r b0d5e752c0c5 tool-data/dbkeys.loc --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/dbkeys.loc Tue Sep 05 09:51:56 2017 -0400 @@ -0,0 +1,13 @@ +#This file lists the locations and dbkeys of all the genome files +#See here for details: http://bedtools.readthedocs.io/en/latest/content/general-usage.html#genome-file-format +#You can add elements to this data table using the data_manager_fetch_genome_dbkeys_all_fasta data manager. +#Alternatively, index files created using samtools faidx (http://www.htslib.org/doc/faidx.html) contain 3 extra columns, but work fine. +#This file has the format (white space characters are TAB characters): +# +# +# +#So, dbkeys.loc could look something like this: +# +#apiMel3 Honeybee (Apis mellifera): apiMel3 /path/to/genome/apiMel3/apiMel3.len +#hg19 Human (Homo sapiens): hg19 Canonical /path/to/genome/hg19/hg19canon.len +# diff -r cb98f29aa1dc -r b0d5e752c0c5 tool_data_table_conf.xml.sample --- a/tool_data_table_conf.xml.sample Fri Aug 18 15:32:52 2017 -0400 +++ b/tool_data_table_conf.xml.sample Tue Sep 05 09:51:56 2017 -0400 @@ -8,5 +8,10 @@ value, dbkey, name, path -
+ + + + value, name, len_path + +
diff -r cb98f29aa1dc -r b0d5e752c0c5 unionBedGraphs.xml --- a/unionBedGraphs.xml Fri Aug 18 15:32:52 2017 -0400 +++ b/unionBedGraphs.xml Tue Sep 05 09:51:56 2017 -0400 @@ -1,4 +1,4 @@ - + combines coverage intervals from multiple BEDGRAPH files macros.xml @@ -10,10 +10,8 @@ unionBedGraphs $header -filler "${filler}" - #if $zero.value == True: - -empty - -g $genome - #end if + $empty.empty_selector + @GENOME_FILE_UNION@ #if str($tag.tag_select) == "tag": #set files = '" "'.join( [ str( $file ) for $file in $tag.inputs ] ) @@ -43,14 +41,20 @@ - + - - + + + + + + + + + + @@ -61,27 +65,28 @@ - + - + - + - - + + - + - - + + + - + @@ -94,8 +99,9 @@ - - + + +