Thanks Youri, That's the bit of code I was worrying about (BAM grooming). I've filed an issue for this: https://github.com/galaxyproject/galaxy/issues/1336 Peter On Wed, Dec 16, 2015 at 9:47 AM, Y. Hoogstrate <y.hoogstrate@erasmusmc.nl> wrote:

Great news Peter,

This is how Galaxy seems to do it:

command = "samtools sort %s %s" % ( file_name, tmp_sorted_dataset_file_name_prefix ) proc = subprocess.Popen( args=command, shell=True, cwd=tmp_dir, stderr=open( stderr_name, 'wb' ) )

https://github.com/galaxyproject/galaxy/blob/dev/lib/galaxy/datatypes/binary...

Since galaxy ships with a builtin pysam (0.8.3) wheel, I think it would be handier to use pysam for sorting (once updated to the latest HTSlib). Pysam now also has all samtools c functions directly builtin:

pysam.index( "arg1", "arg2" ) is the equivalent of the command line argument *samtools index "arg1" "arg2"*

...

Youri