# HG changeset patch # User lparsons # Date 1397243128 14400 # Node ID 227f9d3f0e323f259f67e828e023c3c4e0bde5b4 # Parent 3ffe4e2572a742b9573f5cec44ac3b412e1de8a6 Updated HTSeq package to version 0.6.1, fixed input format string, updated dependency definitions diff -r 3ffe4e2572a7 -r 227f9d3f0e32 .hgignore --- a/.hgignore Mon Jan 27 14:12:46 2014 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -^package$ -syntax: glob -*.pyc diff -r 3ffe4e2572a7 -r 227f9d3f0e32 .hgtags --- a/.hgtags Mon Jan 27 14:12:46 2014 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -142ec9462fa644cd29642ea0ad10e9889ce75619 0.3 -3e7b1e8a69f9d79b8de603f24e466deb7b192100 0.3.release2 -496a7dc0786f8367264520dfff6e8e8e744de3b2 0.3.1 diff -r 3ffe4e2572a7 -r 227f9d3f0e32 README.md --- a/README.md Mon Jan 27 14:12:46 2014 -0500 +++ b/README.md Fri Apr 11 15:05:28 2014 -0400 @@ -9,7 +9,7 @@ Installation ------------ -Installtion directly from the [Galaxy +Installation directly from the [Galaxy Toolshed](http://toolshed.g2.bx.psu.edu/view/lparsons/htseq_count) is recommended. @@ -24,7 +24,7 @@ License ------- -Copyright (c) 2012-2013, Lance R. Parsons +Copyright (c) 2012-2014, Lance R. Parsons All rights reserved. Licensed under the BSD 2-Clause License: diff -r 3ffe4e2572a7 -r 227f9d3f0e32 fabfile.py --- a/fabfile.py Mon Jan 27 14:12:46 2014 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,26 +0,0 @@ -from fabric.api import local -from fabric.operations import prompt -import os - - -def package(): - ''' - Pacakge for upload to toolshed - packaging 'test' version (default) uses local directory - otherwise, specify a mercurial tag to package - ''' - package_dir = 'package' - base_filename = os.path.join(package_dir, 'htseq-count') - version = prompt("Enter version number for package [test]:") - revision_option = '' - if version != '': - revision_option = '-r "%s"' % version - else: - version = 'test' - version_filename = '%s_%s.tar.gz' % (base_filename, version) - local('mkdir -p %s' % package_dir) - local('rm -f %s' % version_filename) - if version == 'test': - local('tar czvf %s --exclude "fabfile.*" --exclude "%s" --exclude ".hg*" *' % (version_filename, package_dir)) - else: - local('hg archive -t tgz %s -X "fabfile.*" -X "package" -X ".hg*" -p . "%s"' % (revision_option, version_filename)) diff -r 3ffe4e2572a7 -r 227f9d3f0e32 htseq-count.xml --- a/htseq-count.xml Mon Jan 27 14:12:46 2014 -0500 +++ b/htseq-count.xml Fri Apr 11 15:05:28 2014 -0400 @@ -1,11 +1,10 @@ - + - Count aligned reads in a BAM file that overlap features in a GFF file htseq-count -h | grep version | sed 's/^\(.*\)*\(version .*\)\./\2/' numpy - htseq + htseq samtools - picard ##set up input files @@ -18,18 +17,9 @@ #set $reference_fasta_filename = str( $samout_conditional.reference_source.ref_file.fields.path ) #end if #end if - #if str($singlepaired) == "paired": - ln -s $samfile local_input.sam && - java -Xmx2G -jar "\$JAVA_JAR_PATH/SortSam.jar" VALIDATION_STRINGENCY=LENIENT SORT_ORDER=queryname O=prepared_input.sam I=local_input.sam TMP_DIR="${__new_file_path__}" - || echo "Error running Picard MergeSamFiles" >&2 && - #else: - #if $samfile.extension == "bam": - samtools view $samfile | - #else - ln -s $samfile prepared_input.sam && - #end if - #end if htseq-count + --format=$samfile.extension + --order=pos --mode=$mode --stranded=$stranded --minaqual=$minaqual @@ -38,27 +28,14 @@ #if $samout_conditional.samout: --samout=$__new_file_path__/${samoutfile.id}_tmp #end if - #if str($singlepaired) == "paired": - prepared_input.sam - #else: - #if $samfile.extension == "bam": - - - #else: - prepared_input.sam - #end if - #end if + $samfile $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: && samtools view -Su -t ${reference_fasta_filename}.fai $__new_file_path__/${samoutfile.id}_tmp | samtools sort -o - sorted > $samoutfile #end if - - - Paired libraries will be sorted by read name prior to counting. - - - + Mode to handle reads overlapping more than one feature. @@ -72,7 +49,7 @@ - + Skip all reads with alignment quality lower than the given minimum value @@ -188,8 +165,7 @@ The following figure illustrates the effect of these three modes: -.. image:: $PATH_TO_IMAGES/count_modes.png - :width: 500 +.. image:: count_modes.png Strandedness diff -r 3ffe4e2572a7 -r 227f9d3f0e32 tool_dependencies.xml --- a/tool_dependencies.xml Mon Jan 27 14:12:46 2014 -0500 +++ b/tool_dependencies.xml Fri Apr 11 15:05:28 2014 -0400 @@ -2,24 +2,20 @@ - + - - - - - + - + - https://pypi.python.org/packages/source/H/HTSeq/HTSeq-0.5.4p5.tar.gz + https://pypi.python.org/packages/source/H/HTSeq/HTSeq-0.6.1.tar.gz - +