生信技能27 - BAM文件的批量合并

最新推荐文章于 2025-05-05 19:09:27 发布

生信与基因组学

最新推荐文章于 2025-05-05 19:09:27 发布

阅读量1.2k

点赞数

分类专栏：生信分析项目实战技能合集文章标签：数据分析数据挖掘 python

本文链接：https://blog.csdn.net/LittleComputerRobot/article/details/133343641

版权

生信分析项目实战技能合集专栏收录该内容

75 篇文章 ¥29.90 ¥99.00

订阅专栏

超级会员免费看

本文介绍了如何利用Samtools的merge命令批量合并BAM文件，包括将当前目录下所有BAM文件合并以及按特定条件进行合并，重点讲解了如何合并指定区域的BAM文件。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Samtools merge命令

当有多个样本的bam文件时，可以使用samtools的merge命令将这些bam文件进行合并为一个bam文件。合并如果为sort后的bam文件将保持现有排序顺序。

merge命令格式：

 samtools merge [-nur1f] [-h inh.sam] [-R reg] [-b <list>] <out.bam> <in1.bam> [<in2.bam><in3.bam>…<inN.bam]

 参数：

    -l 指定压缩等级；

    -b FILE 输入文件列表，一个文件一行；

    -f 强制覆盖同名输出文件；

    -h FILE 指定FILE内的’@’头复制到输出bam文件中并替换输出文件的文件头。否则，输出文件的文件头将从第一个输入文件复制过来；

    -n 设定输入比对文件是以read名进行排序的而不是以染色体坐标排序的；

    -R STRING 合并输入文件的指定区域；

    -r 使RG标签添加到每一个比对文件上，标签值来自文件名；

    -u 输出的bam文件不压缩；

    -c 当多个输入文件包含相同的@RG头ID时，只保留第一个到合并后输出的文件。当合并多个相同样本的不同文件时，非常有用。

    -p 与-c参数类似，对于要