fawnshao的个人博客分享 http://blog.sciencenet.cn/u/fawnshao

博文

UCSC genome browser 镜像搭建

已有 6920 次阅读 2013-1-3 15:28 |个人分类:基因组工具|系统分类:科研笔记| 搭建, Browser

实验室测了几个ChIP-seq,想要把自己的结果展示出来,用UCSC genome browser的自定义track毕竟还是不方便,于是就开始硬着头皮去看本地化。

http://genome.ucsc.edu/admin/mirror.html 这里开始。

具体参考http://genome-source.cse.ucsc.edu/gitweb/?p=kent.git;a=tree;f=src/product

主要问题:

1. MYSQLLIBS 按照文档中写的不一定能用。

2. trashlink很重要,位置错了就不能显示图。

3. mysql的权限要搞对,不然导入数据时总出错。

4. 已有的数据库用rsyncfetchFullGbdb.sh下载就好了,要注意位置,我把他们下载在服务器其他地方,做了个link,方便管理。

5. 搞不清楚数据库样子的时候,去看看源数据库是怎么搞的:mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A

6. 创建自定义的数据库后,先授权。

7. 如果是在已有的基因组中增加自定义的track,允许hgTrackDb更新后才会在网页中显示出来,config文件就参考已有的照猫画虎,多试几次才知道该写些什么。

8. 如果是创建新的基因组,在hgcentral表中新增自定义基因组的信息,最少需要grptrackDbhgFindSpecchromInfogoldgap这几个表,不知道这些表需要什么信息的时候,找个已知数据库看看就知道了。

9. 遇到个错误:mySQL error 1036: Table 'history' is read only,原来是mysql目录权限问题,chown -R mysql:mysql /var/lib/mysql; etc/init.d/mysqld restart就好了。

10. 我要把ChIP-seqpeak显示在服务器上,用wigEncodemacs生成的wig转成browserwig文件,用hgLoadWiggle导入对应的基因组数据库,在hgTrackDb配置文件中写好typecolor之类的参数,就好了。


自定义的基因组显示,就可以对外展示自己的ChIP-seq结果了



https://blog.sciencenet.cn/blog-824692-649161.html

上一篇:用R画多系列的bar
下一篇:基因组坐标
收藏 IP: 158.182.150.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-10-4 21:16

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部