问答文章1 问答文章501 问答文章1001 问答文章1501 问答文章2001 问答文章2501 问答文章3001 问答文章3501 问答文章4001 问答文章4501 问答文章5001 问答文章5501 问答文章6001 问答文章6501 问答文章7001 问答文章7501 问答文章8001 问答文章8501 问答文章9001 问答文章9501

ChIP-seq数据质控与过滤

发布网友 发布时间:2022-11-06 17:34

我来回答

1个回答

热心网友 时间:2023-11-17 22:44

早期生成的fastq数据,它的质量值是基于Phred 64(Illumina 1.3 and 1.4),现在的版本使用的Phred 33 (Illumina 1.8+)

使用fastqc进行质控,它的结果有几个关注点:

主要针对低质量reads和接头,注意:如果要比较不同的样本,要保持它们各自过滤前后一样的reads长度,避免给比对率引入人为因素(长短不一的两条reads,本身就不能放在一起比较)

如果DN*段比测序读长还短,那么得到的reads就会包含下游的接头序列。而接头序列可能会影响比对结果(不过一般接头都是重复序列,有的也很难比对到基因组,但是去掉总比带着好)

不同的过滤软件都会有和接头 stringency 相关的参数设置,比如reads和接头最小的重叠碱基数、最多的错配数。当设置一个比较小的stringency值,就保证最为严格,能检测绝大多数的接头。

比如 trim_galore 的这个参数(默认是非常严格:数值1):

大部分的ChIP-seq数据都是短读长,去低质量不是必须的。但是fastqc图中大量碱基质量如果存在明显的下降,那么就需要去掉

常规的去除方法是:对每条read,从检测到Q值低于某个阈值(比如Q<20)开始去除,直到恢复正常的Q值。但如果只是某个碱基的质量低比较低,这个还能忍受,可以设置一个滑窗而不用一个碱基一个碱基地去看,通过检测滑窗的平均质量值是否低于阈值来判断是不是要去掉这个滑窗中的碱基。

另外,如果存在不同长度的几组数据,可以直接指定固定长度(比如:为了比较50nt和100nt的单端ChIP-seq数据,可以将所有的reads长度都剪到50nt)

如果是分析X或Y染色体的allele-specific binding,它需要许多可靠的SNVs,因此需要更高的Q值来过滤,并且过滤标准更严格,不建议使用滑窗式的过滤

运行时会把所有的参数列出来,方便参考学习,看到其中就有 --fastqc 的命令,另外还有一个参数: -j 它是多线程运行,默认一个线程(需要基于python3)

另外会对每个样本自动构建一个同名目录,存放它相应的数据

看看过滤后数据量的变化,因为原来数据质量就不错,所以也没:

看看过滤后10个样本的GC变化:

声明声明:本网页内容为用户发布,旨在传播知识,不代表本网认同其观点,若有侵权等问题请及时与本网联系,我们将在第一时间删除处理。E-MAIL:11247931@qq.com
太平盛世专辑曲目 梦见自家灶台缺个角 如何看待宝马3系? - 知乎 为什么缺爱的女人婚姻难幸福快乐 求一些好看的架空言情小说,最好是以第一人称为视角的。小白文请绕。 带双目相机的无人机有哪些 禄莱1982年后 工业3d相机国内有哪些靠谱的厂商,特别是在阳光和弱光环境 海康全新双目单线相机 | 消盲区、抑杂光,引领3D检测新突破 10米范围内实现零盲区大FoV深度测量,奥比中光发布Gemini 2双目... 猫眼我的日程在哪 基金持仓成本价是买入时的时间点价格还是当天收盘时的价格? 万德福超市圆形贴纸什么意思 维生素ce什么时间吃最好 维生素ce什么时候吃好些 台式电脑,dj,声音太大怎么调? DJ包房怎样调音量 instruct sb how to do 他命令我完成这项工作.(instruct sb to do sth) 装修的时候需要注意什么事项? 留学保证金是多少?要存银行多久?英国留学 出国留学,在中国银行存,最低要存够多少天,才给开存款证明?去英国留学。 万合集团医疗保险交了吗? 邯郸万合公司的保险有保障吗? 邯郸市万合集团倒闭了吗? 2021宁夏退休人员创城奖金发放时间 职工医疗保险只能报销住院的医疗费用吗? 做梦梦到老情人 Microsoft Office 2010能不能免费发传真 XP怎么样免费发传真啊 ChIP-Seq数据挖掘系列-5.1: ngs.plot 可视化ChIP-Seq 数据 中北大学科研助理工资 河北大学科研助理岗能不能转正 事业单位社会招聘科研助理有编制吗 梓字取名的寓意是什么好不好?梓字取名的寓意女孩名字最佳组合 带木字旁的字男孩名字寓意好的 带木字旁的字适合男孩取名的有哪些 足少阴肾经的经脉循行 55斤汽油等于多少升 想问50升汽油相当于多少斤 ping域名拼不出来,提示找不到主机,是什么原因 伯换一个偏旁组成新字,再组词. 做梦梦见小猫是什么意思 梦见猫的解说 请问百词斩属于腾讯软件吗?联通大王卡能免流量吗? 在微信上听百词斩需要费用吗? 湖州埭溪老虎滩钓鱼免费吗?那里可以自助烧烤吗?烧烤工具自己带。 怎么把花的香味提炼出来 60秒避难所电脑版怎么下载 12123驾驶证锁定怎么解锁? 再婚一方去世后财产怎么分 再婚夫妻共同财产一方去世怎么分割