opensource/galaxy
3
0
Fork 0
mirror of https://github.com/galaxyproject/galaxy.git synced 2026-09-24 16:30:27 +08:00
Code Issues Packages Projects Releases Wiki Activity

Allow BAM's set_meta() to use samtools 1 to generate the index, if the samtools found on $PATH is samtools 1.

Browse Source
This commit is contained in:
Nate Coraor
2015-01-28 10:38:22 -05:00
parent b05ffd2d52
commit 13e8b3e7de
1 changed files with 11 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
lib/galaxy/datatypes/binary.py
+11 -2
View File
@@ -264,10 +264,19 @@ class Bam( Binary ):
##$ samtools index
##Usage: samtools index <in.bam> [<out.index>]
stderr_name = tempfile.NamedTemporaryFile( prefix = "bam_index_stderr" ).name
command = 'samtools index %s %s' % ( dataset.file_name, index_file.file_name )
proc = subprocess.Popen( args=command, shell=True, stderr=open( stderr_name, 'wb' ) )
command = [ 'samtools', 'index', dataset.file_name, index_file.file_name ]
proc = subprocess.Popen( args=command, stderr=open( stderr_name, 'wb' ) )
exit_code = proc.wait()
#Did index succeed?
if exit_code == -6:
# SIGABRT, most likely samtools 1.0+ which does not accept the index name parameter.
command = [ 'samtools', 'index', dataset.file_name ]
proc = subprocess.Popen( args=command, stderr=open( stderr_name, 'wb' ) )
exit_code = proc.wait()
if os.path.exists( os.path.join( dataset.file_name, '.bai' ) ):
shutil.move( os.path.join( dataset.file_name, '.bai' ), index_file.file_name )
else:
open( stderr_name, 'ab+' ).write( 'Galaxy attempted to build the BAM index with samtools 1.0+ but failed\n')
stderr = open( stderr_name ).read().strip()
if stderr:
if exit_code != 0: