njaved/samtools_bedtools

samtools 1.11 + bedtools 镜像文档

镜像概述

本Docker镜像基于miniconda3基础镜像构建，集成了samtools 1.11和bedtools两款生物信息学常用工具。samtools用于处理SAM/BAM格式的高通量测序比对数据，bedtools则专注于基因组区域的交并集分析及相关操作。镜像提供了便捷的生物信息学数据分析环境，适用于基因组学研究中的数据处理流程。

核心功能与特性

samtools 1.11 功能

• SAM/BAM文件格式转换与处理（排序、索引、合并、拆分）
• 比对质量统计与过滤（碱基质量、映射质量筛选）
• 特定区域数据提取（按染色体、位置范围提取reads）
• BAM文件索引创建（.bai格式索引）
• 计算测序深度与覆盖度

bedtools 功能

• 基因组区域交、并、差集运算（intersect、union、subtract）
• 覆盖度分析（计算基因组区域被reads覆盖的深度与长度）
• 格式转换（BED与其他基因组格式互转）
• 窗口分析（滑动窗口统计基因组特征）
• 配对区域分析（寻找配对的基因组区域关系）

镜像特性

• 基于miniconda3，提供Python环境支持
• 工具版本固定（samtools 1.11），确保分析结果可重复
• 轻量级设计，减少冗余依赖
• 支持直接挂载数据目录，方便本地文件操作

使用场景

• 基因组测序数据（WGS、WES、RNA-seq）比对结果处理
• 变异检测流程中的BAM文件预处理（排序、索引、去重）
• 基因区域注释与功能分析
• 测序覆盖度评估与质量控制
• 比较基因组学中的区域差异分析

使用方法

基本使用

拉取镜像

docker
docker pull [镜像仓库地址]/samtools1.11-bedtools:latest

运行容器

docker
docker run -it --rm -v /本地数据目录:/data samtools1.11-bedtools:latest /bin/bash

• -v /本地数据目录:/data：将本地数据目录挂载到容器内的/data目录，实现文件共享
• --rm：容器退出后自动删除
• -it：交互式终端模式

工具使用示例

samtools示例：BAM文件排序与索引

docker
# 进入容器后执行
samtools sort /data/input.bam -o /data/sorted.bam
samtools index /data/sorted.bam

bedtools示例：计算基因区域覆盖度

docker
# 进入容器后执行
bedtools coverage -a /data/genes.bed -b /data/sorted.bam > /data/coverage.txt

Docker Compose配置（可选）

创建docker-compose.yml文件：

yaml
version: '3'
services:
  samtools-bedtools:
    image: [镜像仓库地址]/samtools1.11-bedtools:latest
    volumes:
      - ./local_data:/data
    tty: true
    stdin_open: true

启动服务：

bash
docker-compose up -d

进入容器：

bash
docker-compose exec samtools-bedtools /bin/bash

注意事项

• 确保本地数据目录权限正确，避免容器内无法读取/写入文件
• 大型BAM文件处理建议增加容器内存限制（--memory参数）
• 工具具体参数可通过samtools --help或bedtools --help查看详细说明