BWA:短读长序列比对工具

生命科学 2026-10-27 6 阅读

BWA 是基于 Burrows-Wheeler 变换的 DNA 序列比对软件,用于将短读长序列比对到大型参考基因组,如人类基因组。

这个工具是什么

BWA 是一个将 DNA 序列比对到大型参考基因组的软件包,由 Li H. 和 Durbin R. 开发,采用 Burrows-Wheeler 变换构建 FM-index 实现高效比对。它包含三种算法:BWA-backtrack 针对 100bp 以下的 Illumina 读长,BWA-SW 和 BWA-MEM 支持 70bp 至数兆碱基的较长序列。BWA-MEM 是最新算法,支持长读长和嵌合比对,通常推荐使用。BWA 以 GPLv3 许可证发布,主要用 C 语言编写,仅依赖 zlib 即可编译。

核心能力

  • 基于 FM-index 实现参考基因组的高效索引与比对
  • 提供 BWA-backtrack、BWA-SW 和 BWA-MEM 三种算法
  • BWA-MEM 支持长读长和嵌合比对,速度与准确性较好
  • 输出标准 SAM 格式,便于下游分析

科研中的典型用法

  • 将 Illumina 短读长序列比对到人类参考基因组,用于变异检测前的比对步骤
  • 对 70bp 至数兆碱基的较长序列进行比对,支持嵌合比对分析
  • 使用 BWA-MEM 算法处理克隆序列或组装 contig 的比对任务

快速上手

从 GitHub 克隆源码后执行 make 编译,生成单个可执行文件 bwa。使用前需用 bwa index ref.fa 构建参考基因组索引,然后通过 bwa mem ref.fa read1.fq read2.fq 进行双端比对,输出 SAM 格式。详细用法可查看 man ./bwa.1。

生态与相近工具

同领域有 minimap2(适用于长读长)、BWA-MEM2 和 BWA-MEM3(提供 BWA-MEM 兼容模式)、minibwa(将替代 BWA 用于短读长)。若需 BWA-MEM 的原始输出,可选用 BWA-MEM2 或 BWA-MEM3 兼容模式;新项目可考虑 minibwa。

项目信息

  • 仓库:lh3/bwa
  • 星标:1,767 · 派生:581
  • 主要语言:C
  • 许可证:GPL-3.0
  • 最近提交:2026-08-07

本文属于《学科研究工具知识库 · 生命科学》第 8 篇(按 GitHub 星标排序)。收录标准:该学科公开可获取的开源研究工具,星标与活跃度为主要参考,不代表对其性能或适用性的背书;选型前请结合自身场景评估。

评论 共 0 条

暂无评论

微信小程序

微信扫一扫体验

立即
投稿

微信公众账号

微信扫一扫加关注

发表
评论
返回
顶部