# HG changeset patch # User iuc # Date 1484663185 18000 # Node ID 88cea300225a137fd1e4fba7c698f5f326a6cf37 # Parent 0c01ed8d0e60d09a8ee7b0145e9f7b7daf05893d planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/htseq_count commit da411ad3afa5b8fb7d28fdb1b13a777854e87212 diff -r 0c01ed8d0e60 -r 88cea300225a htseq-count.xml --- a/htseq-count.xml Fri Sep 23 07:48:44 2016 -0400 +++ b/htseq-count.xml Tue Jan 17 09:26:25 2017 -0500 @@ -1,8 +1,8 @@ - + - Count aligned reads in a BAM file that overlap features in a GFF file htseq - samtools + samtools @@ -23,26 +23,26 @@ #if $samout_conditional.samout == "Yes": #if str( $samout_conditional.reference_source.reference_source_selector ) == "history": ln -s "${samout_conditional.reference_source.ref_file}" "${reference_fasta_filename}" && - samtools faidx "${reference_fasta_filename}" 2>&1 || echo "Error running samtools faidx for htseq-count" >&2 && + samtools faidx '${reference_fasta_filename}' 2>&1 || echo "Error running samtools faidx for htseq-count" >&2 && #else: #set $reference_fasta_filename = str( $samout_conditional.reference_source.ref_file.fields.path ) #end if #end if - + #if $force_sort == "True": #if $samfile.extension == 'bam': - samtools sort -n "$samfile" "name_sorted_alignment" && + samtools sort -n -o 'name_sorted_alignment.bam' '$samfile' && #else - samtools view -Su -t "${reference_fasta_filename}.fai" "$samfile" | samtools sort -n - "name_sorted_alignment" && + samtools view -Su -t '${reference_fasta_filename}.fai' '$samfile' | samtools sort -n -o 'name_sorted_alignment.bam' - && #end if #end if - + htseq-count --mode=$mode --stranded=$stranded --minaqual=$minaqual - --type="$featuretype" - --idattr="$idattr" + --type='$featuretype' + --idattr='$idattr' #if $samout_conditional.samout == "Yes": --samout='$__new_file_path__/${samoutfile.id}_tmp' #end if @@ -55,16 +55,16 @@ --format=$samfile.extension '$samfile' #end if - + "$gfffile" | awk '{if ($1 ~ "no_feature|ambiguous|too_low_aQual|not_aligned|alignment_not_unique") print $0 | "cat 1>&2"; else print $0}' > '$counts' 2> '$othercounts' - + #if $samout_conditional.samout == "Yes": && samtools view -Su - -t "${reference_fasta_filename}.fai" - "$__new_file_path__/${samoutfile.id}_tmp" - | samtools sort -o - name_sorted_alignment > "$samoutfile" + -t '${reference_fasta_filename}.fai' + '$__new_file_path__/${samoutfile.id}_tmp' + | samtools sort -o '$samoutfile' - #end if ]]> @@ -179,7 +179,6 @@ - @@ -285,7 +284,6 @@ Public License v3. Part of the 'HTSeq' framework. ]]> - 10.1093/bioinformatics/btu638 diff -r 0c01ed8d0e60 -r 88cea300225a test-data/htseq-test_samout.bam Binary file test-data/htseq-test_samout.bam has changed