刚写题目时发现,现在的分析应该是RNA-seq与ChIP-seq杂合在一起的分析了,为了省事还接着在这个系列写吧。最近有站友发私信建议说把代码跑出来的结果放上,方便学习。这里说明一下:关于学习笔记我尽量把图片放上来,对于分析的记录,因为涉及到实验室的课题,为了避免不必要的麻烦就不放了,本意这个记录就是自己的试验记录,记录下代码过程,为以后学生发表文章参考。
1. 根据gene_list从hg19_genes.gtf文件中提取gene_TSS信息
也可以使用biomaRt包进行注释,中间出错尚未解决,有机会再尝试。
2. 根据gene_TSS信息利用bedtools getfast从hg19_genes.fasta文件中提取promoter的序列
输出结果如下:

分析软件meme使用conda install命令进行安装【这个软件安装费了很大劲】
建立python=2.7的新环境
安装python=2.7的版本
启动相关环境并安装特定版本的meme软件和支持软件icu
查看软件帮助以检查安装效果
声明:本站部分文章及图片源自用户投稿,如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢!