1990 年前,序列比对只能通过 Smith-Waterman DP(O(mn))完成。随数据库增长,这已不可接受。BLAST 引入了一套启发式加速策略,在可接受的敏感性损失下将速度提升数十倍。BLAST 的 seed-and-extend 范式直接影响了后来的 BLAT、BWA-SW 和 BWA-MEM。
| 特性 | BLAST X-dropoff | BWA-MEM Z-dropoff |
|---|---|---|
| 公式 | bestScore − curr > X | bestScore − curr > Z + |x−y|×pgapExt |
| 单边 gap | 惩罚 | 不惩罚 |
| 设计目标 | 蛋白质搜索 | DNA 比对(含结构变异) |
Z-dropoff 是 X-dropoff 的改进:加入 |x−y|×pgapExt 后,坐标差增大时阈值自动放宽,允许一方有长插入。
| 算法 | 种子 | 种子长度 | 延伸 | 截断 |
|---|---|---|---|---|
| BLAST | 精确 k-mer | 11 (DNA) | 无空位 → SW | X-dropoff |
| BLAT | 精确 k-mer | 11 | 无空位 → SW | 无 |
| BWA-SW | DP 间 FM-index | 变长,允许差异 | SW | Z-best |
| BWA-MEM | SMEM | 变长,精确 | Banded DP | Z-dropoff |