python - 将 SeqIO 字典写入 Fasta 文件

Question

我最初使用 Bio.SeqIO.to_dict 语句将我的 fasta 序列转换为字典。我想将子集字典写回 fasta 文件。

Test 是一个 Python 字典，以 fasta header 作为键，以序列作为索引。

我的代码尝试利用 SeqIO.write:

with open("example.fasta", "w") as handle:
SeqIO.write(test, handle, "fasta")

AttributeError: 'str' object has no attribute 'id'

恐怕通过将 SeqIO 生成器对象转换为字典，我无法轻松返回其他 SeqIO 函数期望的输入。

Answer 1

我解决了这个问题。即使使用 SeqIO.to_dict 转换，字典的值也是原始生成器类。要将这个字典写回 Fasta，只需调用字典的值即可。

with open("example.fasta", "w") as handle:
    SeqIO.write(test.values(), handle, "fasta")