生信软件3 - mapping比对bam文件质量评估工具 qualimap

最新推荐文章于 2025-08-18 09:48:01 发布
原创 最新推荐文章于 2025-08-18 09:48:01 发布 · 1.4k 阅读 · 4 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#数据挖掘 #python #pandas #numpy

Qualimap是一款强大的生物信息学工具，它提供了对BAM文件的质量控制（QC）功能，包括单个样本、RNA-seq样本和多样本的统计分析。用户可以通过简单的命令行接口进行操作，生成详细的报告和可视化图形。此外，文章还提到了其他相关生信软件如visNano和picard。

qualimap简介

qualimap可用于统计bam文件，输出结果包含可视化图形和详细统计信息，以及每条contig的mapping信息。

软件官网及案例

http://qualimap.conesalab.org/
案例
http://qualimap.conesalab.org/doc_html/samples.html#bam-samples

安装

wget https://bitbucket.org/kokonech/qualimap/downloads/qualimap_v2.2.1.zip

# 解压
unzip qualimap_v2.2.1.zip

# 进入软件目录
cd qualimap_v2.2.1
# 查看帮助信息
./qualimap -h

软件模块

bamqc: 用于单个bam文件的QC统计
rnaseq: 用于转录组RNA-Seq样本的bam文件的QC统计
multi-bamqc: 用于多样本的bam文件分组QC统计
counts: 用于转录组数据计数的统计，用于量化表达水平
clustering: 用于表观遗传特征的聚类
comp-couts: 输入bam文件和注释文件，计算映射到每个区域reads的数量

运行

1. 执行bamqc模块命令

bamqc模块用于单个NGS样本bam文件的统计。

qualimap bamqc -bam sample.bam \  # 指定bam文件路径
					-outformat PDF:HTML \  # 输出文件格式
					-outdir out \        #输出文件目录
					-nt 12 \             # 线程数
					--java-mem-size=10G  #设置最大内存

程序运行结束后，统计信息在report.pdf文件中查看。

report.pdf中部分图表如下所示：

Summary - Globals

2. 执行rnaseq模块命令

rnaseq模块用于RNA-seq数据的bam文件的统计。

qualimap rnaseq -bam sample.bam \  # 指定bam文件路径
					-outformat PDF:HTML \  # 输出文件格式
					-outdir out \        #输出文件目录
					-nt 12 \             # 线程数
					--java-mem-size=10G  #设置最大内存

3. 执行multi-bamqc模块命令

multi-bamqc模块用于多样本NGS的bam文件的统计和比较。

qualimap multi-bamqc -r \ # -r指定输入bam文件
					-d qualimap.list \     # 输入文件列表
					-outformat PDF:HTML \  # 输出文件格式
					-outdir out \        #输出文件目录
					-nt 12 \             # 线程数
					--java-mem-size=10G  #设置最大内存

其中输入文件列表qualimap.list有三列，每行一个样本，第一列样品名称，第二列 包含路径的bam文件/bamqc结果目录，第三列组名。

生信软件文章推荐

生信软件1 - 测序下机文件比对结果可视化工具 visNano

生信软件2 - 下游比对数据的统计工具 picard

RNA-seq分析：Step1（软件的安装及配置）
沉香best的博客
08-19 3488
好几个月没有跑RNA-seq分析了，为防止遗忘，特整理分析流程于转录组分析专栏。RNA-seq分析是生信分析流程比较入门的操作，常规的分析主要包括差异基因表达分析、GO和KEGG分析和WGCNA分析。一般来说，前两个分析是最为常见的，第三个主要存在于纯转录组分析文章中。本文主要讲述转录组分析的第一步，软件的安装和配置。
生信分析进阶3 - pysam操作bam文件统计unique reads和mapped reads高级技巧合辑
LittleComputerRobot的博客
05-31 1822
pysam统计unique reads和mapped reads高级技巧
bamcov:快速计算和可视化比对文件中的序列覆盖率
03-18
┌┐ ┌─┐┌┬┐┌─┐┌─┐┬ ┬ Quick calculation and visualization ├┴┐├─┤││││ │ │└┐┌┘ of sequence coverage on the terminal └─┘┴ ┴┴ ┴└─┘└─┘ └┘ v0.1 安装 git clone --recurse-submodules https://github.com/fbreitwieser/bamcov cd bamcov make make test 确保使用--recurse-submodules克隆存储库以获取htslib，否则必须在克隆的存储库中键入git submodule update --init --recursive 。 用法 显示带有标题的表格输出（默认）： ./bamcov -H test.ba
BamQC:映射的质量控制分析程序
05-18
BamQC BAM文件的质量控制应用程序 BamQC是一个应用程序，它使用包含映射数据的BAM文件并对其进行一系列测试以生成全面的QC报告。 这将帮助您理解数据，并告诉您所分析的文件是否存在异常。 每个测试都被标记为通过，警告或失败，具体取决于它与正常的大型数据集的期望值相差多少，且没有明显的偏差。 需要强调的是，警告甚至故障不一定表示您的数据存在问题，而仅是不寻常的问题。 样品的生物学性质可能意味着您期望结果中存在这种特殊偏差。 交互式图形或命令行 BamQC既可以作为交互式图形应用程序运行，也可以让您在单个应用程序中查看多个文件的结果。 或者，您可以以非交互方式运行该程序（例如，作为管道的一部分），该程序将为您处理的每个文件生成HTML报告。 跨平台Java BamQC是一个用Java编写的跨平台应用程序。 从理论上讲，它应该在具有合适的Java运行时环境的任何平台上运行。 已在运
[samtools] tview查看mapping比对信息
weixin_34233618的博客
03-17 1305
samtools view -bS test.sam > test.bam samtools sort test.bam -o test_sorted samtools index -b test_sorted.bam samtools tviewtest_sorted.bam ref.fa 转载于:https://www.cnblogs.com/xiaofeiIDO/p/65...
bam文件批量统计-qualimap multi-bamqc
qq_42680661的博客
03-01 942
Qualimap 2 是一个用 Java 和 R 编写的独立于平台的应用程序，它提供了图形用户界面 （GUI） 和命令行界面，以促进对齐排序数据及其衍生物（如特征计数）的质量控制。
QualiMap：一款强大的二代测序比对文件质控工具
公众号/简说基因，知乎/简宝玉
09-26 1560
在生物信息学中，数据质量的评估和可视化是很重要的一环。今天我们来聊聊一个常用的工具——Qualimap，它是一个用于评估高通量测序数据质量的开源软件，尤其是对RNA-seq和DNA测序数据的分析非常友好。无论你是本科生还是刚接触生物信息学的新人，Qualimap 都能帮助你快速上手了解你的测序数据质量。Qualimap的功能特点多样化的模块支持Qualimap 提供多个模块，主要包括：•BAM ...
生信软件49 - 全基因组亚硫酸氢盐测序（WGBS）比对与甲基化水平调用工具BSseeker2
LittleComputerRobot的博客
07-31 550
全基因组亚硫酸氢盐测序（WGBS）比对与甲基化水平调用工具BSseeker2
生信软件47 - 超低测序深度的全基因组测序cfDNA肿瘤分数估计工具ichorCNA
LittleComputerRobot的博客
06-24 921
超低测序深度的全基因组测序cfDNA肿瘤分数估计工具ichorCNA
生信软件46 - 三代测序低深度全基因组测序结构变异SV检测工具NanoVar
LittleComputerRobot的博客
05-06 852
三代测序低深度全基因组测序结构变异SV检测工具NanoVar
minion_multiQC:用于小兵数据的多功能QC
05-16
用于牛津纳米Kong的Minion Sequencer数据的多功能QC 什么？ 单个脚本适用于在牛津纳米Kong的Minion测序仪的单个流通池的输出上使用一系列流行的工具执行QC。 为什么？ 因为有很多好的工具，而且如果您真的想知道自己的数据是否有用，那么一次使用其中的许多工具将很有用。 依赖项/使用的工具 喝咖啡-首先要安装一些东西... 长鳍金枪鱼： : minion_QC.R： : NGMLR： ： qualimap： ://qualimap.bioinfo.cipf.es/ samtools： ： miniasm： : 配额： : //quast.sourceforge.net/download.html 安装所有这些工具，然后确保所有内容都在路径中的某个位置，以便您可以成功运行每个工具，如下所示： read_fast5_basecaller
生信软件51 - 生物信息高效命令行find,cut, xargs和seqtk技巧
LittleComputerRobot的博客
08-18 987
生物信息高效命令行find,cut, xargs和seqtk技巧
Qualimap 对比对好的排序后的bam文件的统计分析
weixin_41869644的博客
03-26 4568
1.Qualimap这个软件输入的bam必须是sort后的。 2.使用命令如下： qualimap --java-mem-size=35G bamqc -bam /export/rawdata/0.NGS/BJXWZ-201812001/DM18A2292-1/outdir/2.BwaMapOutDir/DM18A2292-1.sort.bam -outdir bamqc_result -o...
linux查看bam文件格式,bam格式文件处理大全（四）
weixin_39703551的博客
05-13 1256
sam文件是短序列比对生成的文件，是二代测序中最核心的文件。在RNAseq，变异检测等分析中，都需要首先生成sam文件格式。bam文件是sam格式的二进制格式，转换为二进制之后，可以减小文件的存储。掌握sam/bam文件的操作是处理二代测序数据的非常重要的内容，例如sam与bam的转换，排序，建立索引，reads计数等等操作。16 统计bam并绘图samtoolsstatsA1.sorted...
转录组比对番外篇：aligner 与 mapping 区别
冷冻工厂
12-06 696
导读 考虑到在转录组比对时，有许多软件可以使用，但很少有介绍它们之间的差别。因此，本文主要介绍 STAR, KALLISTO, SALMON 之间的区别。 1. 定义 STAR 是 alignerKallisto/salmon 是 mapping 2. aligner vs. mapping 那么aligner 与 mapping的区别[1]是什么呢？ 2.1. aligner 当我们比对一个 read 时，我们不仅要询问它可能在基因组中的位置，还要询问碱基与碱基的确切对应关系。例如，我们想要得到类似这样
mapping、bam结果提取、覆盖率、数据深度、序列拼接
Curry_chenhu的博客
07-25 4256
提取比对到参考序列结果 samtools view -bF 4 abc.bam > abc.F.bam 配对比对 samtools view -bF 12 abc.bam > abc.F.bam #bam排序 samtools sort *.bam -o *.bam bam转fq bedtools bamtofastq -i * -fq *_1.fq -fq2 *_2.fq samtools fastq NOf41_human-mapped.bam > NOf41_human_test.f
各种bam文件统计质量值的软件
Cassiel60的博客
04-04 8031
1.bamdst -- a BAM Depth Stat Tool 对bam文件进行统计 安装bamdst https://github.com/shiquan/bamdst git clone https://github.com/shiquan/bamdst.git cd bamdst make make后直接help后，查看用法，有时候可能会提示bamdst不存在，本来已经存在...
samtools常用命令详解
weixin_33883178的博客
03-05 1466
samtools的说明文档：http://samtools.sourceforge.net/samtools.shtml samtools是一个用于操作sam和bam文件的工具合集。包含有许多命令。以下是常用命令的介绍 1. view view命令的主要功能是：将sam文件转换成bam文件；然后对bam文件进行各种操作，比如数据的排序(不属于本命令的功能)和提取(这些操作 是对bam文...
第二章 RNA-seq一般分析流程全套
热门推荐
laurelgirl的博客
02-22 3万+
提前申明：该贴参照于https://www.jianshu.com/p/e8cd62ba14fe 标题1. 用conda安装RNA-seq所需要 #启动conda自设环境 conda activate python3 或者用 source activate python3 #安装所需软件(conda可以同时安装多个软件，但是建议初学者还是选择逐一安装，避免出现错误) conda inst...
生信基础概念之unique reads VS multi-mapping reads
最新发布
09-02
在生信领域，unique reads（唯一比对的读段）和 multi - mapping reads（多重比对的读段）是在将测序得到的短序列（reads）比对到参考基因组时出现的两种不同情况。 Unique reads 指的是那些在参考基因组中只有一个最佳比对位置的 reads。这些 reads 能够明确地定位到基因组的某一个特定区域，因此在后续的分析中，它们所携带的信息通常比较可靠，能够为基因表达定量、变异检测等分析提供准确的数据支持。例如在基因表达定量分析中，unique reads 可以准确地反映某个基因的转录本丰度。 Multi - mapping reads 则是那些可以比对到参考基因组中多个不同位置的 reads。这种情况通常是由于基因组中存在重复序列，或者基因家族成员之间具有高度的序列相似性导致的。由于这些 reads 无法明确其来源，在后续分析中使用时会带来一定的困难。例如在基因表达定量分析中，很难确定这些 reads 究竟是来自哪个基因。 在代码示例中，以 BWA（Burrows - Wheeler Aligner）比对工具的输出为例，SAM（Sequence Alignment/Map）格式文件中，对于 unique reads 和 multi - mapping reads 有不同的标记。以下是简单的 Python 代码示例，用于统计 unique reads 和 multi - mapping reads 的数量： ```python unique_count = 0 multi_count = 0 with open('alignment.sam', 'r') as f: for line in f: if line.startswith('@'): continue fields = line.strip().split('\t') # 这里简单假设 MAPQ 值为 0 表示 multi - mapping reads if int(fields[4]) == 0: multi_count += 1 else: unique_count += 1 print(f"Unique reads count: {unique_count}") print(f"Multi - mapping reads count: {multi_count}") ```
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

生信与基因组学

每一份鼓励是我坚持下去动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值