Biopython 是一个用 Python 编写的生物信息学工具包，其中包括了许多序列读取与处理的模块。以下是 Biopython 序列读取与处理常用函数的示例及解释。

1. 读取 FASTA 格式文件中的序列

from Bio import SeqIO

# 读取单个序列
record = SeqIO.read("example.fasta", "fasta")

# 读取多个序列
records = list(SeqIO.parse("example.fasta", "fasta"))

SeqIO.read 函数用于读取单个序列，返回一个 SeqRecord 类型的变量。SeqIO.parse 函数用于读取多个序列，返回一个可迭代的 SeqRecord 类型变量。

2. 读取 GenBank 格式文件中的序列

from Bio import SeqIO

# 读取单个序列
record = SeqIO.read("example.gb", "genbank")

# 读取多个序列
records = list(SeqIO.parse("example.gb", "genbank"))

与读取 FASTA 格式文件类似，SeqIO.read 函数用于读取单个序列，返回一个 SeqRecord 类型的变量。SeqIO.parse 函数用于读取多个序列，返回一个可迭代的 SeqRecord 类型变量。

3. 写入 FASTA 格式文件

from Bio import SeqIO

seq = "ATCGATCG"

with open("example.fasta", "w") as output:
    SeqIO.write(SeqRecord(Seq(seq)), output, "fasta")

SeqIO.write 函数用于将序列写入文件，并指定文件格式。此处将字符串 seq 包装为一个 SeqRecord 变量写入到文件中。

4. 写入 GenBank 格式文件

from Bio import SeqIO

seq = "ATCGATCG"

with open("example.gb", "w") as output:
    SeqIO.write(SeqRecord(Seq(seq)), output, "genbank")

与写入 FASTA 格式文件类似，SeqIO.write 函数用于将序列写入文件，并指定文件格式。此处将字符串 seq 包装为一个 SeqRecord 变量写入到文件中。 免责声明：本文内容通过AI工具匹配关键字智能整合而成，仅供参考，火山引擎不对内容的真实、准确或完整作任何形式的承诺。如有任何问题或意见，您可以通过联系service@volcengine.com进行反馈，火山引擎收到您的反馈后将及时答复和处理。