# HG changeset patch # User devteam # Date 1452290529 18000 # Node ID 6069ffa8b240444532d33726ef45fe0a94d717b0 # Parent 7bee165c3a446a489302fd5ce0efba05ed31dd4e planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd diff -r 7bee165c3a44 -r 6069ffa8b240 bwa-mem.xml --- a/bwa-mem.xml Thu Dec 03 17:38:14 2015 -0500 +++ b/bwa-mem.xml Fri Jan 08 17:02:09 2016 -0500 @@ -1,5 +1,5 @@ - + - map medium and long reads (> 100 bp) against reference genome read_group_macros.xml @@ -8,41 +8,26 @@ +/dev/null`; ## Linux + if [ $? -ne 0 ]; then + size=\$(stat -f %z "${reference_source.ref_file}"); ## OSX + fi && + if [ "\$size" -lt 2000000000 ]; then + echo "Reference genome size is \$size bytes, generating BWA index with is algorithm"; + bwa index -a is "${reference_fasta_filename}"; + else + echo "Reference genome size is \$size bytes, generating BWA index with bwtsw algorithm"; + bwa index -a bwtsw "${reference_fasta_filename}"; + fi + ) && #else: #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path ) #end if @@ -132,9 +117,10 @@ "${fastq_input.fastq_input1}" #end if - | samtools view -Sb - > temporary_bam_file.bam && + | samtools view -Sb - > temporary_bam_file.bam && samtools sort -f temporary_bam_file.bam ${bam_output} +]]> diff -r 7bee165c3a44 -r 6069ffa8b240 bwa.xml --- a/bwa.xml Thu Dec 03 17:38:14 2015 -0500 +++ b/bwa.xml Fri Jan 08 17:02:09 2016 -0500 @@ -1,5 +1,5 @@ - + - map short reads (< 100 bp) against reference genome read_group_macros.xml @@ -68,38 +68,26 @@ +/dev/null`; ## Linux + if [ $? -ne 0 ]; then + size=\$(stat -f %z "${reference_source.ref_file}"); ## OSX + fi && + if [ "\$size" -lt 2000000000 ]; then + echo "Reference genome size is \$size bytes, generating BWA index with is algorithm"; + bwa index -a is "${reference_fasta_filename}"; + else + echo "Reference genome size is \$size bytes, generating BWA index with bwtsw algorithm"; + bwa index -a bwtsw "${reference_fasta_filename}"; + fi + ) && #else: #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path ) #end if @@ -134,7 +122,7 @@ "${input_type.fastq_input1}" #end if - > first.sai && + > first.sai && bwa aln -t "\${GALAXY_SLOTS:-1}" @@ -149,7 +137,7 @@ "${input_type.fastq_input2}" #end if - > second.sai && + > second.sai && bwa sampe @@ -178,7 +166,7 @@ "${reference_fasta_filename}" "${input_type.fastq_input1}" - > first.sai && + > first.sai && bwa samse @@ -202,7 +190,7 @@ "${reference_fasta_filename}" "${input_type.bam_input}" - > first.sai && + > first.sai && bwa aln -t "\${GALAXY_SLOTS:-1}" @@ -211,7 +199,7 @@ @command_options@ "${reference_fasta_filename}" "${input_type.bam_input}" - > second.sai && + > second.sai && bwa sampe @@ -238,7 +226,7 @@ "${reference_fasta_filename}" "${input_type.bam_input}" - > first.sai && + > first.sai && bwa samse @@ -251,9 +239,10 @@ "${reference_fasta_filename}" first.sai "${input_type.bam_input}" #end if - | samtools view -Sb - > temporary_bam_file.bam && + | samtools view -Sb - > temporary_bam_file.bam && samtools sort -f temporary_bam_file.bam ${bam_output} +]]>