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