我正在编写一个处理 BLAST 结果的脚本,每次都必须重新请求搜索,这在时间上相当麻烦,因为我仍在学习该库。 目前我正在使用:
blastResult = NCBIWWW.qblast("blastp", "swissprot", AAsequence)
from Bio.Blast import NCBIWWW
from Bio.Blast import NCBIXML
# Perform the BLAST search
result_handle = NCBIWWW.qblast("blastp", "swissprot", "KAF8065798.1")
# Parse the BLAST result and retrieve the query_id
blast_record = NCBIXML.read(result_handle)
query_id = blast_record.query
# You can use the query_id to retrieve the BLAST request later
# For example, you can use the query_id to check the status of the BLAST request
# Close the result handle
result_handle.close()
# Example query_id
print(query_id)
# Output: "Query_1"
Lipase [Scenedesmus sp. PABB004]