[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