# HG changeset patch # User iuc # Date 1504640414 14400 # Node ID 44867b59dbf2a0c8883ebd6cef64bccb23483708 # Parent e0cec48a4695b0baab45f8bcf9753efd72be7940 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit cf7a3674bc833bbd5fdd2ad02e724935ffc7a174 diff -r e0cec48a4695 -r 44867b59dbf2 bedToBam.xml --- a/bedToBam.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/bedToBam.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 bedpeToBam.xml --- a/bedpeToBam.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/bedpeToBam.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 complementBed.xml --- a/complementBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/complementBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 fisherBed.xml --- a/fisherBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/fisherBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 flankBed.xml --- a/flankBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/flankBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 genomeCoverageBed.xml --- a/genomeCoverageBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/genomeCoverageBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 macros.xml --- a/macros.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/macros.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 mapBed.xml --- a/mapBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/mapBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 randomBed.xml --- a/randomBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/randomBed.xml Tue Sep 05 15:40:14 2017 -0400 @@ -1,4 +1,4 @@ - + generate random intervals in a genome macros.xml @@ -8,7 +8,7 @@ - + @@ -28,6 +28,7 @@ + diff -r e0cec48a4695 -r 44867b59dbf2 shuffleBed.xml --- a/shuffleBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/shuffleBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 slopBed.xml --- a/slopBed.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/slopBed.xml Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 test-data/mapBed3.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBed3.bed Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 test-data/mapBed4.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBed4.bed Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 test-data/mapBed_result5.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBed_result5.bed Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 tool-data/dbkeys.loc.sample --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/dbkeys.loc.sample Tue Sep 05 15:40:14 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 e0cec48a4695 -r 44867b59dbf2 tool_data_table_conf.xml.sample --- a/tool_data_table_conf.xml.sample Fri Aug 18 15:33:08 2017 -0400 +++ b/tool_data_table_conf.xml.sample Tue Sep 05 15:40:14 2017 -0400 @@ -8,5 +8,10 @@ value, dbkey, name, path -
+ + + + value, name, len_path + +
diff -r e0cec48a4695 -r 44867b59dbf2 unionBedGraphs.xml --- a/unionBedGraphs.xml Fri Aug 18 15:33:08 2017 -0400 +++ b/unionBedGraphs.xml Tue Sep 05 15:40:14 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 @@ - - + + +