BioPython的Bio.align与Bio.pairwise2对齐实现基于什么算法？

两个模块的核心对齐逻辑均采用生物信息学领域通用的标准动态规划比对算法，无未公开的私有算法，也没有为这两个模块的实现单独发表专属算法论文，具体溯源如下：

Bio.pairwise2 模块 #

这是BioPython早期提供的纯Python实现的成对序列比对模块：

• 全局序列对齐功能基于Needleman-Wunsch算法实现，是1970年提出的经典全局双序列比对动态规划算法
• 局部序列对齐功能基于Smith-Waterman算法实现，是1981年提出的经典局部双序列比对动态规划算法
• 支持线性、仿射两类空位罚分模式，其中仿射空位罚分逻辑基于Gotoh在1982年提出的经典改进方案实现，可分别设置空位开放罚分与空位延伸罚分
• 模块同时支持自定义替换矩阵、自定义匹配/错配得分，所有逻辑均严格遵循上述经典算法的标准流程，无算法层面的自定义修改

Bio.Align 模块 #

这是BioPython 1.75版本后推出的新版对齐模块，用于逐步替代旧的Bio.pairwise2模块：

• 核心成对全局/局部对齐的算法逻辑与Bio.pairwise2完全一致，同样基于Needleman-Wunsch、Smith-Waterman经典算法，仿射空位罚分沿用Gotoh 1982年的改进方案
• 区别仅为工程层面的优化：核心计算逻辑用C语言重写，相比纯Python实现的Bio.pairwise2运行速度有数量级提升，同时优化了大序列比对的内存占用，算法本身没有改动
• 模块额外封装的多序列比对读写、比对结果操作等工具类功能，不涉及核心对齐算法的修改

可引用学术文献 #

科研写作中可根据使用的功能直接引用对应领域通用文献即可：

• 全局成对对齐引用：Needleman S B, Wunsch C D. A general method applicable to the search for similarities in the amino acid sequence of two proteins. Journal of Molecular Biology, 1970, 48(3): 443-453.
• 局部成对对齐引用：Smith T F, Waterman M S. Identification of common molecular subsequences. Journal of Molecular Biology, 1981, 147(1): 195-197.
• 使用仿射空位罚分模式可补充引用：Gotoh O. An improved algorithm for matching biological sequences. Journal of Molecular Biology, 1982, 162(3): 705-708.
• 引用BioPython工具本身可使用官方正式发表文献：Cock P J A, Antao T, Chang J T, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics, 2009, 25(11): 1422-1423.

注：上述算法溯源在BioPython官方维护的文档与源码注释中均有明确标注，不存在未公开的算法调整，所有实现均遵循领域通用标准。

内容的提问来源于stack exchange，提问作者luca sala