methylKit 是一个用于分析甲基化测序数据的R包,不仅支持WGBS
,RRBS
和目的区域甲基化测序,还支持oxBS-sq
, TAB-seq
等分析5hmc
的数据。 其核心功能是差异甲基化分析和差异甲基化位点和区域的注释。
安装过程如下:
source(“http://bioconductor.org/biocLite.R“)
biocLite(“methylKit”)
推荐使用最新版本的R进行安装,这样可以使用最新版本的methylKist
。
利用methylKit
做差异分析包括3步
1. 读取原始数据
每个样本一个原始数据,methylKit
支持两种格式的methylation calling
文件
纯文本格式
内容如下
每一行是一个甲基化位点,coverage
代表覆盖这个位点的reads数,freqC
代表甲基化C的比例,freqT
代表非甲基化C的比例。这种纯文本格式内容非常直观,文件大小相比bam 文件小很多,读取的速度更快。
纯文本格式的读取过程如下