FastQC的基本介绍

FastQC的基本介绍:

FastQC是一款基于Java的软件,它可以快速地对测序数据进行质量评估,其官网为:http://www.bioinformatics.babraham.ac.uk/projects/fastqc/

主要参数:

-o --outdir FastQC生成的报告文件的储存路径,生成的报告的文件名是根据输入来定的

-t --threads 选择程序运行的线程数,

-q --quiet 安静运行模式,一般不选这个选项的时候,程序会实时报告运行的状况

主要格式:

fastqc -o outdir -t threads fastq1 fastq2 ...

image.png

运行结束后生成两个文件一个.html网页文件,一个是.zip压缩文件。

我们使用浏览器打开html文件,然后显示这样的:

image.png

结果介绍:

Summary信息:

 上图中Summary的部分就是整个报告的目录,整个报告分成若干个部分。合格会有个绿色的对勾,警告是个黄色的叹号,不合格是个红色的叉子。
  • 基本信息
image.png

Filename:指的是进行质控的文件名

Encoding:指测序平台的版本和相应的编码版本号

# Total Sequences:指reads的数量

Sequence length:指测序的长度

# %GC 指整体序列中的GC含量

  • 序列测序质量统计
image.png

此图中的横轴是测序序列第1个碱基到第151个碱基

纵轴是质量得分,Q = -10*log10(error P)即20表示0.01的错误率,30表示0.001

图中红线表示中值

图中蓝色的细线是各个位置的平均值的连线

# Warning 报警 如果任何碱基质量低于10,或者是任何中位数低于25

# Failure 报错 如果任何碱基质量低于5,或者是任何中位数低于20

  • 每条序列的测序质量统计
image.png

序列长度为151bp,那么这151个位置每个位置Q值的平均值就是这条reads的质量值

# 该图横轴是0-40,表示Q值

纵轴是每个值对应的reads数目

这个样本数据,测序结果主要集中在30-40中,证明测序质量很好!

  • 碱基分布
image.png

横轴是1 - 151 bp;纵轴是百分比

# 图中四条线代表A T C G在每个位置平均含量

理论上来说,A和T应该相等,G和C应该相等,但是一般测序的时候,刚开始测序仪状态不稳定,很可能出现上图开头的情况。

  • 序列平均GC含量分布图

  • image.png

# 横轴是0 - 100%; 纵轴是每条序列GC含量对应的数量

蓝色的线是程序根据经验分布给出的理论值,红色是真实值,两个应该比较接近才比较好

  • N统计含量
image.png

当测序仪器不能辨别某条reads的某个位置都是ATCG哪个碱基时,就会产生"N",对所有reads的每个位置统计N的比率。

Warning 报警 如果任意位置的N比例超过5%

Failure 报错 如果任意位置的N比例超过20%

  • 序列测序长度统计
image.png

# 每次测序仪测出来的长度在理论上应该是完全相等的,但是总会有一些偏差

比如此图中,151bp是主要的,但是还是有少量的150和152bp的长度

# 当测序的长度不同时,如果很严重,则表明测序仪在此次测序不成功

  • 重复序列
image.png

统计序列完全一样的reads的频率。横坐标是duplication的次数,纵坐标是duplicated reads的数目。

它的说明文档看的不是特别明白,大家可以自己看看

#http://www.bioinformatics.babraham.ac.uk/projects/fastqc/Help/3%20Analysis%20Modules/8%20Duplicate%20Sequences.html

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,590评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 86,808评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,151评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,779评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,773评论 5 367
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,656评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,022评论 3 398
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,678评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 41,038评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,659评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,756评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,411评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,005评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,973评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,203评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,053评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,495评论 2 343

推荐阅读更多精彩内容