2

我修改了如下的一段代码来解析来自 BLAST XML 输出的所需信息。

import csv
from Bio.Blast import NCBIXML
blast_records = NCBIXML.parse(open('PGblast.xml', 'rU'))

output = csv.writer(open('PGhit.csv','w'), delimiter =',',
                    quoting=csv.QUOTE_NONNUMERIC)
output.writerow(["Query","Hit ID", "Hit Def", "E-Value"])

E_VALUE_THRESH = 0.00000000000000001

for blast_record in blast_records:
    for alignment in blast_record.alignments:
        for hsp in alignment.hsps:
            if hsp.expect < E_VALUE_THRESH:
                output.writerow([blast_record.query[:8],
                                 alignment.hit_id, alignment.hit_def,hsp.expect])

blast_records.close()

该代码允许我解析 E 值被截断的命中。但是，我希望仅解析 BLAST XML 输出中的最佳命中或前 3 个命中，因为 BLAST 输出文件的大小很大。

解析每个命中结果将花费大量时间来处理，实际上我不希望所有命中结果。

有人可以帮我吗？

4

1 回答 1

1

仅解析每个 Hit 的前 3 个 HSP，而不解析整个文件，需要您编写自己的自定义 XML 解析器。Biopython 的 NCBIXML 不这样做。

但是，如果您正在寻找速度改进，您可以尝试新的 SearchIO 子模块 (http://biopython.org/wiki/SearchIO)。它有一个新的 BLAST XML 解析器，它应该比旧的 NCBIXML 解析器更快。旧的解析器依赖于纯 Python XML 解析器，而 SearchIO 中的新解析器则cElementTree尽可能使用。

该子模块仍然是新的和实验性的，因此在正式发布之前可能还会有一些变化。如果您有兴趣，这里还有一个草稿教程：http: //bow.web.id/biopython/Tutorial.html#htoc96。

于 2012-12-15T03:59:58.910 回答