这个步骤推荐在R里面做，载入表达矩阵，然后设置好分组信息，统一用DEseq2进行差异分析，当然也可以走走edgeR或者limma的voom流程。 基本任务是得到差异分析结果，进阶任务是比较多个差异分析...
转录组入门（6）： reads计数
要求 实现这个功能的软件也很多，还是烦请大家先自己搜索几个教程，入门请统一用htseq-count，对每个样本都会输出一个表达量文件。 需要用脚本合并所有的样本为表达矩阵。参考：生信编程直播第四题：多...
转录组入门（5）： 序列比对
比对软件很多，首先大家去收集一下，因为我们是带大家入门，请统一用hisat2，并且搞懂它的用法。 直接去hisat2的主页下载index文件即可，然后把fastq格式的reads比对上去得到sam文件...
转录组入门（4）：了解参考基因组及基因注释
准备工作 参考基因组 测序得到的是几百bp的短read， 相当于把拼图打散了给你。如果没有参考基因组，从头(de novo)组装等于是重走人类基因组计划的老路，也就是打散了拼图，却不告诉你原来是什么样...
转录组入门（3）：质量控制
需要用安装好的sratoolkit把sra文件转换为fastq格式的测序文件，并且用fastqc软件测试测序文件的质量！ 数据解压 之前下载了所有的数据，但只有样本915才是mRNA-Seq测序结果，...
转录组入门（2）：读文章拿到测序数据
本系列课程学习的文章是：AKAP95 regulates splicing through scaffolding RNAs and RNA processing factors. Nat Commu...
RNA-seq数据综合分析教程
2 RNA-seq数据分析 mRNA-seq是目前最常用的高通量测序技术，一般的用法就是看看基因表达谱，寻找差异表达的基因。我和高通量测序数据分析结缘，也是因为RNA-seq。 一开始我对mRNA-s...
搞懂illumina nextera Tn5 和ATAC seq adaptor 序列
在准备NGS文库的时候，会有用到转座酶Tn5，Nextera DNA Library Preparation Kit ，比如ATAC-seq就有用到这个Tn5。转座酶携带有特定的序列称为转座子Tran...
本体论和功能分析
为什么需要本体论 作为一位大学统计棉花表皮毛的苦逼生物狗，深刻体会什么叫做经验，也就是人类模式识别能力的强大和不精确性。当时的导师教我如何根据表皮毛的长短和浓密进行基因型的判定，但是我一直纠结长和短，...
如何进行变异检测
什么是基因组变异 基因组变异是一个定义比较模糊的概念. 所谓的变异是相对于一个完美的“参考基因组”而言。但是其实完美的“参考基因组”并不存在，因为我们只是选择某一个物种里的其中似乎比较正常的个体进行测...