# HG changeset patch # User iuc # Date 1484945173 18000 # Node ID 620d5603d1a871bab83090eba13be2311da249a1 # Parent 6e5c95760ab15afb48b543f4aa4ad11de66cd291 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/htseq_count commit 89e3a62dd6cbc8bcec84d08c1710bfb4e7f5938f-dirty diff -r 6e5c95760ab1 -r 620d5603d1a8 htseq-count.xml --- a/htseq-count.xml Wed Sep 21 10:59:41 2016 -0400 +++ b/htseq-count.xml Fri Jan 20 15:46:13 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 6e5c95760ab1 -r 620d5603d1a8 test-data/htseq-test_samout.bam Binary file test-data/htseq-test_samout.bam has changed