[samtools]flagstat命令简介

samtools flagstat命令简介：

统计输入文件的相关数据并将这些数据输出至屏幕显示。每一项统计数据都由两部分组成，分别是QC pass和QC failed，表示通过QC的reads数据量和未通过QC的reads数量。以“PASS + FAILED”格式显示。还可以根据samtools的标志位显示相应的内容，但是这里不做讨论。

命令格式：

samtools flagstat <in.bam> |<in.sam> | <in.cram>

运行flagstat命令得到的结果如下图所示。

从第一行至第十一行分别表示：

1. QC pass的reads的数量为2499917，未通过QC的reads数量为0，意味着一共有2499971条reads；

2. 重复reads的数量，QC pass和failed

3. 比对到参考基因组上的reads数量；

4. paired reads数据数量；

5. read1的数量；

6. read2 的数量；

7. 正确地匹配到参考序列的reads数量；

8. 一对reads都比对到了参考序列上的数量，但是并不一定比对到同一条染色体上；

9. 一对reads中只有一条与参考序列相匹配的数量；

10. 一对reads比对到不同染色体的数量；

11. 一对reads比对到不同染色体的且比对质量值大于5的数量。

如有谬误，望请指出。

[samtools]flagstat命令简介相关推荐

[samtools]index命令简介
samtools index命令的功能描述: 为了能够快速访问bam文件,可以为已经基于坐标排序后bam或者cram的文件创建索引,生成以.bai或者.crai为后缀的索引文件.必须使用排序后的文件, ...
[samtools]mpileup命令简介
samtools的mpileup命令是一个samtools中一个很重要的命令.它的主要功能主要是生成BCF.VCF文件或者pileup一个或多个bam文件.比对记录以在@RG中的样本名作为区分标识符. ...
samtools常用命令详解
samtools的说明文档:http://samtools.sourceforge.net/samtools.shtml samtools是一个用于操作sam和bam文件的工具合集.包含有许多命令.以 ...
samtools 常用命令详解
samtools是一个用于操作sam和bam文件的工具合集.包含有许多命令.以下是常用命令的介绍 1. view view命令的主要功能是:将sam文件转换成bam文件:然后对bam文件进行各种操作, ...
linux wc 命令简介
此wc命令不是让大家没有食欲的地方.而是linux下一个简单的小命令. NAME wc - word, line, character, and byte count SYNOPSIS wc [-cl ...
WINCE cvrtbin命令简介
********************************LoongEmbedded******************************** 作者:LoongEmbedded(kandi ...
find与grep命令简介及正则表达式（转）
find与grep命令简介及正则表达式两个更为有用的命令和正则表达式在我们开始学习新的Shell编程知识之前,我们先来看一下两个更为有用的两个命令,这两个命令虽然并不是Shell的一部分,但是在进 ...
NoSQL数据库Redis使用命令简介
NoSQL 数据库Redis使用命令简介 NOSQL数据库可以按照它们的数据模型分成4类: 1.键-值对存储库(Key-Value) redis-- 2.BigTable实现(BigTable-imp ...
快速上手Linux核心命令（一）：核心命令简介
Linux核心命令系列文章目录快速上手Linux核心命令(一):核心命令简介快速上手Linux核心命令(二):关机.重启快速上手Linux核心命令(三):文件和目录操作命令快速上手Linux核 ...

[samtools]flagstat命令简介

[samtools]flagstat命令简介相关推荐

最新文章

热门文章