python怎么输出序列_BioPython：从Blast输出fi中提取序列id

我用这个代码提取了所有的结果from Bio.Blast import NCBIXML

for record in NCBIXML.parse(open("rpoD.xml")) :

print "QUERY: %s" % record.query

for align in record.alignments :

print " MATCH: %s..." % align.title[:60]

for hsp in align.hsps :

print " HSP, e=%f, from position %i to %i" \

% (hsp.expect, hsp.query_start, hsp.query_end)

if hsp.align_length < 60 :

print " Query: %s" % hsp.query

print " Match: %s" % hsp.match

print " Sbjct: %s" % hsp.sbjct

else :

print " Query: %s..." % hsp.query[:57]

print " Match: %s..." % hsp.match[:57]

print " Sbjct: %s..." % hsp.sbjct[:57]

print "Done"

或是为了少一些细节from Bio.Blast import NCBIXML

for record in NCBIXML.parse(open("NC_003197.xml")) :

#We want to ignore any queries with no search results:

if record.alignments :

print "QUERY: %s..." % record.query[:60]

for align in record.alignments :

for hsp in align.hsps :

print " %s HSP, e=%f, from position %i to %i" \

% (align.hit_id, hsp.expect, hsp.query_start, hsp.query_end)

print "Done"

我用过这个网站