[BioPython] Psi blast
Julie Bernauer
julie.bernauer at ibbmc.u-psud.fr
Mon Apr 25 07:47:50 EDT 2005
Hello !
I hope this question hasn't been answered yet. Please excuse me in case
it has.
I'd like to get the result of phi-blast after 5 runs (html page):
Here is a piece of my code :
while 1:
cur_record=s_iterator.next()
if cur_record is None:
break
f_record=Fasta.Record()
f_record.title=cur_record.accessions[0]
f_record.sequence=cur_record.sequence
print f_record.title
b_results=NCBIWWW.qblast('blastp','swissprot',f_record)
save_file=open('resultats/blast/'+id+'.html','w')
blast_results=b_results.read()
save_file.write(blast_results)
save_file.close()
I's like to have sthg like :
b_results=NCBIWWW.qblast('blastp','swissprot',f_record,run_psiblast=5)
Would you please help me ?
Thanks in advance
--
Julie BERNAUER
Equipe de Génomique Structurale http://www.genomics.eu.org
IBBMC - UMR 8619 - U.P.S. Bât.430 Tel. : +33 1 69 15 31 57
91405 Orsay - FRANCE Fax. : +33 1 69 85 37 15
More information about the BioPython
mailing list