用途：BWA-MEM2 基于 BWT-FM 索引，将 Illumina 短读长 DNA 序列（约 50–250 bp）比对到参考基因组，是 WGS、WES、ChIP-seq、ATAC-seq 等 DNA 比对的标准工具，产出可供 GATK HaplotypeCaller、Picard 使用的 SAM/BAM。
主要能力：产出带 @RG read group 头的 SAM；对双端读长标记 proper pair，并将嵌合/拆分读长解析为补充比对；-M 可将拆分读长标记为 secondary（兼容 Picard），-Y 对补充比对使用软剪切（推荐用于 GATK）。它比原版 BWA-MEM 快约 2 倍且结果相同。
使用方式：先确认环境是否已有该工具，可运行 command -v bwa-mem2；在 pixi 项目中用 pixi run bwa-mem2 调用。安装：conda install -c bioconda bwa-mem2 samtools。建立索引：bwa-mem2 index GRCh38.fa。比对并排序：bwa-mem2 mem -t 16 -R "@RG\tID:sample1\tSM:sample1\tPL:ILLUMINA" GRCh38.fa sample1_R1.fastq.gz sample1_R2.fastq.gz | samtools sort -@ 8 -o sample1.sorted.bam。随后可用 samtools index、fixmate、markdup 去重，再进入 GATK。
适用场景：WGS/WES 变异检测前比对、ChIP-seq/ATAC-seq 产生 BAM 供 MACS3 峰值调用、生成 GATK 兼容 BAM、将旧 FASTQ 重新比对新装配。读长 ≥ 50 bp；更短读长可考虑 BWA-backtrack；跨剪接位点的 RNA-seq 应改用 STAR；需本地比对或最小化索引时可用 Bowtie2。
限制与注意：人类基因组索引约需 28 GB 内存，小鼠约 6–8 GB，索引人类约 25–35 分钟；D2 中的样本名决定下游联合基因分型分组，缺失或重复会静默合并或拆分样本。典型问题包括映射率低、GATK 报 @RG 缺失、索引内存不足、双端计数不平衡（交错 FASTQ 可加 -p）、索引文件缺失、速度慢、补充比对影响下游工具、GATK BQSR 缺已知变异 VCF 等。参数默认与取值范围见资料：-t 默认 1、-k 默认 19、-c 默认 500、-T 默认 30 等。收费与价格暂未提供。