| Previous changeset 25:bf27106652f3 (2017-02-08) Next changeset 27:9f164587a92f (2017-06-06) |
|
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/freebayes commit 0c59a0b32e4875120e697951a187b81c662bc39c |
|
modified:
freebayes.xml |
| b |
| diff -r bf27106652f3 -r a028d13cd860 freebayes.xml --- a/freebayes.xml Wed Feb 08 12:45:05 2017 -0500 +++ b/freebayes.xml Thu May 04 11:34:13 2017 -0400 |
| b |
| @@ -1,4 +1,4 @@ -<tool id="freebayes" name="FreeBayes" version="@DEPENDENCY_VERSION@-2"> +<tool id="freebayes" name="FreeBayes" version="@DEPENDENCY_VERSION@-3"> <description>bayesian genetic variant detector</description> <macros> <import>macros.xml</import> @@ -35,14 +35,25 @@ ln -s -f '${Tabixized_input}' 'input_variant_vcf.vcf.gz.tbi' && #end if - #for $bam_count, $input_bam in enumerate( $reference_source.input_bams ): - samtools view -H b_${bam_count}.bam | - grep "^@SQ" | - cut -f 2- | - awk '{ gsub("^SN:","",$1); - gsub("^LN:","",$2); - print $1"\t0\t"$2; }' >> regions_all.bed && - #end for + ##if user has specified a region or target file, just use instead of calculating a set of unique regions + + #if str( $target_limit_type.target_limit_type_selector ) == "limit_by_target_file": + ln -s '${target_limit_type.input_target_bed}' regions_all.bed && + #elif str( $target_limit_type.target_limit_type_selector ) == "limit_by_region": + printf '${target_limit_type.region_chromosome}\t${target_limit_type.region_start}\t${target_limit_type.region_end}' > regions_all.bed && + #else + ##divide up the regions in the bam file for efficient processing + #for $bam_count, $input_bam in enumerate( $reference_source.input_bams ): + samtools view -H b_${bam_count}.bam | + grep "^@SQ" | + cut -f 2- | + awk '{ gsub("^SN:","",$1); + gsub("^LN:","",$2); + print $1"\t0\t"$2; }' >> regions_all.bed && + #end for + #end if + + sort -u regions_all.bed > regions_uniq.bed && ## split into even small chunks, this has some disatvantages and will not be used for the moment @@ -73,12 +84,6 @@ ## Outputs --vcf './vcf_output/part_\$i.vcf' - #if str( $target_limit_type.target_limit_type_selector ) == "limit_by_target_file": - --targets '${target_limit_type.input_target_bed}' - #elif str( $target_limit_type.target_limit_type_selector ) == "limit_by_region": - --region '${target_limit_type.region_chromosome}:${target_limit_type.region_start}..${target_limit_type.region_end}' - #end if - ##advanced options #if str( $options_type.options_type_selector ) == "simple": ##do nothing as command like build up to this point is sufficinet for simple diploid calling |