AW: [BioPython] Problem with standalone BLAST
Andreas Schlicker
andreasschlicker at web.de
Wed Jun 2 17:26:43 EDT 2004
Hi Brad; hi all;
Here is the code I use:
from Bio.Blast import NCBIStandalone
results, errors = NCBIStandalone.blastall(os.path.join(self.dir, blastall'),
program, database, query, expectation = self.e)
save = file(output, 'w')
while 1:
line = results.readline()
if not line:
break
save.write(line)
results.close()
errors.close()
save.close()
Basically, just run Blast and save it's output to a file.
Most times this works fine but sometimes Blast does not finish. The process
stays alive but is sleeping after a while and never continues.
I tested it on these two different platforms:
SunOS 5.9
Python 2.3.3
BioPython 1.24
2 x 3.06 GHz Xeon / 4 GByte RAM
RedHat 7.3
Python-2.3.4
Numerical Python 23.1
mxTextTools 2.0.5
reportlab 1.19
Biopython 1.30
To me it looks like a problem with Python because running Blast from the
command line with the same command -- as generated from Biopython -- works
perfect.
Has anyone had similar problems? Or can think of a solution?
Thanks for your help.
Andreas
More information about the BioPython
mailing list