[BioPython] Is the bug still there??
Eirik Sønneland
eirik.sonneland at student.umb.no
Tue Feb 22 06:59:09 EST 2005
Thanks, Martina.
You made me see an error and it seems like it runs ok for now...even the
parsing:-)
regards,
Eirik
Martina wrote:
> Hallo Eric,
>
> I just tried the same. It takes a while and produces a warning, but
> this (from the docu) worked:
>
> INFILE = 'fastaTest.txt'
> from Bio import Fasta
> from Bio.Blast import NCBIWWW
>
> file_for_blast = open(INFILE, 'r')
> f_iterator = Fasta.Iterator(file_for_blast)
> f_record = f_iterator.next()
> b_results = NCBIWWW.qblast('blastn', 'month', f_record)
> save_file = open('my_blast.out', 'w')
> blast_results = b_results.read()
> save_file.write(blast_results)
> save_file.close()
>
> So the problem must be the parsing of the results, not the qblast.
More information about the BioPython
mailing list