使用Python的Bio.SeqIO库进行生物序列的合并和比对

Bio.SeqIO是Biopython库中的一个模块，用于处理生物序列文件的输入和输出。它支持常见的文件格式，如FASTA、GenBank、FASTQ等，并提供了一些常用的功能，如序列合并和比对等。

在下面的例子中，我们将演示如何使用Bio.SeqIO进行序列的合并和比对。

首先，我们需要安装Biopython库。可以使用以下命令在终端中安装：

pip install biopython

安装完成后，我们就可以开始使用Bio.SeqIO库了。

1. 序列合并

假设我们有两个FASTA文件，分别是seq1.fasta和seq2.fasta，文件内容如下：

seq1.fasta:

>seq1
ATCG

seq2.fasta:

>seq2
GATT

现在我们要将这两个序列合并成一个文件，代码如下：

from Bio import SeqIO

records = []
for file in ['seq1.fasta', 'seq2.fasta']:
    with open(file, 'r') as f:
        record = SeqIO.read(f, 'fasta')
        records.append(record)

merged_record = records[0] + records[1]

with open('merged.fasta', 'w') as f:
    SeqIO.write(merged_record, f, 'fasta')

代码中，我们首先创建一个空列表records用于存储待合并的序列对象。然后，我们使用SeqIO.read()函数读取FASTA文件，并将读取到的序列对象存入records列表中。最后，通过"+"运算符将两个序列对象合并成merged_record，并使用SeqIO.write()函数将合并后的序列写入到merged.fasta文件中。

2. 序列比对

假设我们有两个FASTA文件，分别是ref.fasta和query.fasta，文件内容如下：

ref.fasta:

>ref
ATCG

query.fasta:

>query
AGGT

现在我们要对ref和query序列进行比对，代码如下：

from Bio import SeqIO
from Bio import pairwise2

# 读取参考序列和查询序列
with open('ref.fasta', 'r') as f:
    ref_record = SeqIO.read(f, 'fasta')
with open('query.fasta', 'r') as f:
    query_record = SeqIO.read(f, 'fasta')

# 进行序列比对
alignments = pairwise2.align.globalxx(ref_record.seq, query_record.seq)

# 打印比对结果
for alignment in alignments:
    print(pairwise2.format_alignment(*alignment))

代码中，我们首先使用SeqIO.read()函数读取ref.fasta和query.fasta中的序列对象。然后，我们使用pairwise2.align.globalxx()函数对ref序列和query序列进行全局比对，得到一系列比对结果。最后，我们通过pairwise2.format_alignment()函数将比对结果格式化为字符串，并打印到屏幕上。

以上就是使用Bio.SeqIO进行生物序列的合并和比对的示例。Bio.SeqIO库提供了丰富的功能，可以用于处理各种生物序列相关的任务，如读取、写入、合并、比对等。通过灵活运用这些功能，我们可以更方便地处理和分析生物序列数据。