Skip to content
Snippets Groups Projects
Commit 00f77e8a authored by Studer Gabriel's avatar Studer Gabriel
Browse files

use default argument for db size when creating new blast db

Newer version have a hard limit for maximum file size that has been
exceeded by the arguments provided by the blast binding.
We simply use the default instead. Thanks to Andrius Merkys for the proposed fix.
parent 15f67e43
Branches
Tags
No related merge requests found
...@@ -149,21 +149,20 @@ def CreateDB(infasta, dbout, mkdb_cmd=None): ...@@ -149,21 +149,20 @@ def CreateDB(infasta, dbout, mkdb_cmd=None):
if mkdb_cmd==None: if mkdb_cmd==None:
try: try:
exe=settings.Locate('formatdb') exe=settings.Locate('formatdb')
args=[exe, '-i', infasta, '-v',str(10000),'-n',dbout] args=[exe, '-i', infasta, '-n', dbout]
except: except:
try: try:
exe=settings.Locate('makeblastdb') exe=settings.Locate('makeblastdb')
args=[exe, '-in', infasta, ' -max_file_sz 10GB -out', dbout,'-dbtype','prot'] args=[exe, '-in', infasta, '-out', dbout, '-dbtype', 'prot']
except: except:
raise RuntimeError('could not find makeblastdb/formatdb executable') raise RuntimeError('could not find makeblastdb/formatdb executable')
else: else:
if os.path.basename(mkdb_cmd)=='makeblastdb': if os.path.basename(mkdb_cmd)=='makeblastdb':
exe=settings.Locate('makeblastdb',explicit_file_name=mkdb_cmd) exe=settings.Locate('makeblastdb',explicit_file_name=mkdb_cmd)
args=[exe, '-in', infasta, ' -max_file_sz 10GB -out', dbout,'-dbtype','prot'] args=[exe, '-in', infasta, '-out', dbout, '-dbtype', 'prot']
elif os.path.basename(mkdb_cmd)=='formatdb': elif os.path.basename(mkdb_cmd)=='formatdb':
exe=settings.Locate('formatdb',explicit_filename=mkdb_cmd) exe=settings.Locate('formatdb',explicit_filename=mkdb_cmd)
args=[exe, '-i', infasta, '-v',str(10000),'-n',dbout] args=[exe, '-i', infasta, '-n', dbout]
else: else:
raise IOError('mkdb command must either be the path to formatdb or makeblastdb!') raise IOError('mkdb command must either be the path to formatdb or makeblastdb!')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment