如何轻松打造个人专属的IGV Browser

Integrative Genomics Viewer(IGV)是一个强大的基因组可视化工具,可以交互式的查看大多数的基因组相关数据,并且支持多种NGS测序数据类型。在大多数NGS数据分析中都会涉及到分析数据的可视化展示,IGV无疑为我们提供了一个便利的可视化途径。IGV不仅提供了基于Java的本地版集成工具,还可以定制化的构建web版的IGV Browser,方便地存储并展示个人的专属数据。

安装Apache服务器

首先需要下载并安装Apache服务器搭建web环境,这里以CentOS系统为基础安装apache服务器,安装完成后会生成/var/www/html文件夹,默认开启80端口,可以通过服务器IP(http://xxx.xxx.xxx.xx)在Chrome浏览器上进行访问。

# 利用CentOS系统的包管理器yum直接安装httpd服务器
yum install -y httpd
cd /var/www/html
cat >index.html <<EOF
<html>
<head>
</head>
hello world
</html>
EOF

下载并安装igv-webapp

image.png

IGV Web App提供了一个官方的示例网站,可以通过https://igv.org/app/网址进行访问。
image.png

igv-webapp的网站源代码可以在GitHub上进行下载安装。

cd /var/www/html
git clone https://github.com/igvteam/igv-webapp
cd igv-webapp/
# 新建data和log文件夹,用于存放基因组数据和web访问日志文件
mkdir data log
# 创建web配置文件igv.conf,开放1220端口,ServerName设置为自己服务器的IP
cat >igv.conf <<EOF
Listen 1220
<VirtualHost *:1220>
    ServerName xxx.xxx.xxx.xx
    ServerAdmin 1369852697@qq.com
    ProxyRequests Off
    <Proxy *>
        Order deny,allow
        Allow from all
    </Proxy>
    DocumentRoot /var/www/html/igv-webapp
    ErrorLog /var/www/html/igv-webapp/log/igv.error.log
    LogLevel info
    CustomLog /var/www/html/igv-webapp/log/igv.access.log combined
</VirtualHost>
EOF
# 将igv.conf文件放置到/etc/httpd/conf.d/文件夹中
mv igv.conf /etc/httpd/conf.d/

安装完成后在Chrome浏览器中输入服务器IP地址和端口http://xxx.xxx.xxx.xx:1220,即可成功访问自己的IGV Web browser,默认IGV browser提供了多个物种的基因组数据信息。

image.png

IGV browser默认提供的hg19的参考基因组信息,可以通过igvwebConfig.js配置文件进行更改

image.png

配置个人专属基因组数据

IGV browser的基因组数据是通过resources文件下的genomes.json文件进行读入访问的,可视化的track数据信息存放在resources文件下的tracks文件夹中,通过trackRegistry.json文件进行读入访问,默认会从中读取数据展示,访问数据的加载速度可能会比较慢。
查看genomes.json文件

genomes.json

查看trackRegistry.json文件
trackRegistry.json

下载并存放个人的基因组数据到data文件中,更改resources文件中的genomes.json文件

cd /var/www/html/igv-webapp/data
# 创建个人参考基因组文件夹
mkdir hg19 && cd hg19
# 下载参考基因组相关文件
wget https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg19/hg19.fasta
wget https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg19/hg19.fasta.fai
wget https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg19/cytoBand.txt
wget https://s3.amazonaws.com/igv.org.genomes/hg19/refGene.sorted.txt.gz
wget https://s3.amazonaws.com/igv.org.genomes/hg19/refGene.sorted.txt.gz.tbi
# 更改genomes.json文件
cd ../../resources
vim genomes.json

将URL后的地址改为自己的服务器IP端口和数据存放的地址


image.png

同样的,track的相关信息(如bam、bed,vcf和gtf等文件)也可以存放在data文件夹中,添加或更改tracks文件夹下对应的json文件。除了在服务器中存放相关的数据,IGV browser还可以直接从本地,远端服务器或链接地址导入相关数据进行可视化展示。


image.png

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

推荐阅读更多精彩内容

  • ORA-00001: 违反唯一约束条件 (.) 错误说明:当在唯一索引所对应的列上键入重复值时,会触发此异常。 O...
    我想起个好名字阅读 5,167评论 0 9
  • 国家电网公司企业标准(Q/GDW)- 面向对象的用电信息数据交换协议 - 报批稿:20170802 前言: 排版 ...
    庭说阅读 10,863评论 6 13
  • 现在的自己到底算是在经历什么 应该说现在的自己一直是活在自己的恐惧当中 仔细想来的话应该是自己对自己没有自信 你害...
    百花梨落阅读 307评论 2 3
  • 广州易常青企业宣传片
    卖报小话家阅读 74评论 0 0
  • 今天是冬至节,要不是各大群的饺子汤圆红包雨,我还真忘记了这么一个节日。 群也许是现在跟生活联系最紧密的一部分了,...
    悠然花开73阅读 141评论 3 1