ChIP-Seq分析

Introduction to ChIP-Seq using high-performance computing https://hbctraining.github.io/Intro-to-ChIPseq/schedule/2-day.html
ATAC-seq分析实操生信技能树健明教程 https://www.jianshu.com/p/09e05bcd6981
ChIP-seq流程结果文件解读 https://www.cnblogs.com/leezx/p/14897113.html
Sample GSM4212624 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM4212624
Chip-seq的22个疑问——【Jimmy】 https://www.jianshu.com/p/e10b951959bb
ATAC-seq 分析(上) https://www.jianshu.com/p/b7bd87315b6e

注意:MultiQC合并超过100个report时,网页版不再是鼠标放置显示名称及信息
Bowtie2使用方法与参数详细介绍 https://blog.csdn.net/herokoking/article/details/77847384

如何得到unique mapping reads

Bowtie2 manul也说了,MAPQ筛选即可
HISAT2 Bowtie2 提取唯一比对 unique mapping reads https://www.jianshu.com/p/bc3751b72f0c
怎么获取unique mapping read? http://www.biotrainee.com/thread-1115-1-3.html
如何从bam文件提取unique mapping reads呢? http://www.360doc.com/content/21/0714/12/76149697_986499744.shtml

去除PCR重复

samtools官网manual http://www.htslib.org/doc/samtools-markdup.html

samtools collate -o namecollate.bam example.bam
samtools fixmate -m namecollate.bam fixmate.bam
samtools sort -o positionsort.bam fixmate.bam
samtools markdup positionsort.bam markdup.bam                  # 加上参数-r直接去掉duplicate,或者-F 0x400筛选也行

Sambamba也是一个简便的工具,但现在github打不开,根据位置sort后去重和positionsort.bam去重,最后得到的reads数不一样,前者和picard MarkDuplicates一样,后者和samtools markdup一样,可能与sort前对bam文件的处理有关。

关于duplicate是什么的一些理解

MarkDuplicates 的意义与作用 https://blog.csdn.net/weixin_43569478/article/details/108079103
NGS测序中PCR重复序列的判定方法 https://cloud.tencent.com/developer/article/1624512
去重工具到底哪家强 https://www.baidu.com/link?url=MsFZBf4aY9tTitKMN3DX75Fo_J_mI620xmZGYdeSJ9iRgsaucUXQGdSyr0xbCDk-IMGlW45qao5QbOe4YpBxndqpjMiH0MEMNVrIA-UWj9_&wd=&eqid=da5d7cce0006e415000000056124a11b
NGS测序结果中duplicate序列问题的理解 https://blog.csdn.net/ncyujijian/article/details/107371264
NGS数据的Duplication问题 http://blog.sina.com.cn/s/blog_69e75efd0102wu57.html
测序的PCR duplicates - I(转载) https://www.jianshu.com/p/73483070379b

MACS2 callpeak

github上的一些问题讨论 https://github.com/macs3-project/MACS/issues
-f BAMPE 输出无_model.r

image.png

Q: Peakcalling using -f BAMPE vs -f BAM on paired-end reads https://github.com/macs3-project/MACS/issues/456
BAMPE #92 https://github.com/macs3-project/MACS/issues/92

IDR (Irreproducible Discovery Rate)

使用IDR软件处理生物学重复样本的peak calling https://blog.csdn.net/weixin_43569478/article/details/108079783
CHIP-seq流程学习笔记(9)-使用IDR 软件对生物学重复样本间的差异peak进行提取 https://blog.csdn.net/xiaomotong123/article/details/116535067

deeptools

如何使用deeptools处理BAM数据 https://www.jianshu.com/p/b494426ecd32
使用deeptools生成ChIP-seq信号热图与谱图 https://www.jianshu.com/p/1f38674b2475
bigwig归一化方式详解 https://blog.csdn.net/weixin_43569478/article/details/108079478
使用deeptools将bam文件转换为bw文件 https://www.jianshu.com/p/491557586118

ChIPseeker

【YuLabSMU】ChIPseeker被引超过1千
"input object should be a named list" error in ChIPseeker https://support.bioconductor.org/p/106903/
非模式基因GO富集分析:以玉米为例+使用OrgDb https://guangchuangyu.github.io/cn/2017/07/clusterprofiler-maize/#disqus_thread
植物构建txdb对象 https://www.jianshu.com/p/4a92ce5174fb
学习一遍ChIPseeker的使用 https://www.jianshu.com/p/c76e83e6fa57

分析流程

ChIP实验操作以及Chip-seq数据分析 https://www.jianshu.com/p/4f29f2afe3c0?utm_campaign=haruki
MACS2学习笔记 https://www.jianshu.com/p/ceace75a3dcf?utm_campaign=maleskine&utm_content=note&utm_medium=seo_notes&utm_source=recommendation
[软件使用 3] 使用MACS2分析ChIP-seq数据,快速入门! https://zhuanlan.zhihu.com/p/90180058
第3篇:用MACS2软件call peaks https://www.jianshu.com/p/21e8c51fca23
一篇文章学会ChIP-seq分析(上) https://mp.weixin.qq.com/s?__biz=MzAxMDkxODM1Ng==&mid=2247484370&idx=1&sn=3e29df3b621076bd1d33a8ec157858e5&scene=21#wechat_redirect
一篇文章学会ChIP-seq分析(下) https://mp.weixin.qq.com/s?__biz=MzAxMDkxODM1Ng==&mid=2247484372&idx=1&sn=30ad9c4c54f3967000af059eb7a38ace&scene=21#wechat_redirect
一文学会ChIP-Seq数据分析(想想也不可能) https://mp.weixin.qq.com/s?__biz=MzI1MjU5MjMzNA%3D%3D&mid=2247484614&idx=1&sn=196eda169967ef66c62f1bcc61d57194&scene=45#wechat_redirect
ChIP-Seq数据中包含了spike-in怎么分析
Chip-Seq 分析过程二 https://www.jianshu.com/p/66a6eabd2519
每日一生信--FastUniq去除pairedreads的duplicates http://blog.sina.com.cn/s/blog_670445240101lqat.html

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

推荐阅读更多精彩内容