[Biopython] help with NCBIXML.parse

Peter Cock p.j.a.cock at googlemail.com
Wed Mar 28 13:08:28 UTC 2012


On Wed, Mar 28, 2012 at 1:54 PM,  <ferreirafm at usp.br> wrote:
> Hi there,
> What I'm doing wrong in the following piece of code?
> Thanks in advance,
> Fred

You seem to be calling BLAST multiple times in a loop and
trying to give it SeqRecord objects. It wants FASTA files,
and you can call BLAST once with a single FASTA query
file (containing multiple records) and a single database or
FASTA subject file (also containing multiple records).

As to the specific error, did you look at your blast_out.xml
file and what it said on line 88?

Peter



More information about the Biopython mailing list