如何輕松打造個(gè)人專屬的IGV Browser

Integrative Genomics Viewer(IGV)是一個(gè)強(qiáng)大的基因組可視化工具躁劣,可以交互式的查看大多數(shù)的基因組相關(guān)數(shù)據(jù)弹渔,并且支持多種NGS測序數(shù)據(jù)類型舆声。在大多數(shù)NGS數(shù)據(jù)分析中都會(huì)涉及到分析數(shù)據(jù)的可視化展示吞获,IGV無疑為我們提供了一個(gè)便利的可視化途徑供置。IGV不僅提供了基于Java的本地版集成工具谨湘,還可以定制化的構(gòu)建web版的IGV Browser,方便地存儲(chǔ)并展示個(gè)人的專屬數(shù)據(jù)芥丧。

安裝Apache服務(wù)器

首先需要下載并安裝Apache服務(wù)器搭建web環(huán)境紧阔,這里以CentOS系統(tǒng)為基礎(chǔ)安裝apache服務(wù)器,安裝完成后會(huì)生成/var/www/html文件夾续担,默認(rèn)開啟80端口擅耽,可以通過服務(wù)器IP(http://xxx.xxx.xxx.xx)在Chrome瀏覽器上進(jìn)行訪問。

# 利用CentOS系統(tǒng)的包管理器yum直接安裝httpd服務(wù)器
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提供了一個(gè)官方的示例網(wǎng)站物遇,可以通過https://igv.org/app/網(wǎng)址進(jìn)行訪問乖仇。
image.png

igv-webapp的網(wǎng)站源代碼可以在GitHub上進(jìn)行下載安裝。

cd /var/www/html
git clone https://github.com/igvteam/igv-webapp
cd igv-webapp/
# 新建data和log文件夾询兴,用于存放基因組數(shù)據(jù)和web訪問日志文件
mkdir data log
# 創(chuàng)建web配置文件igv.conf乃沙,開放1220端口,ServerName設(shè)置為自己服務(wù)器的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瀏覽器中輸入服務(wù)器IP地址和端口http://xxx.xxx.xxx.xx:1220蕉朵,即可成功訪問自己的IGV Web browser崔涂,默認(rèn)IGV browser提供了多個(gè)物種的基因組數(shù)據(jù)信息。

image.png

IGV browser默認(rèn)提供的hg19的參考基因組信息始衅,可以通過igvwebConfig.js配置文件進(jìn)行更改

image.png

配置個(gè)人專屬基因組數(shù)據(jù)

IGV browser的基因組數(shù)據(jù)是通過resources文件下的genomes.json文件進(jìn)行讀入訪問的冷蚂,可視化的track數(shù)據(jù)信息存放在resources文件下的tracks文件夾中,通過trackRegistry.json文件進(jìn)行讀入訪問汛闸,默認(rèn)會(huì)從中讀取數(shù)據(jù)展示蝙茶,訪問數(shù)據(jù)的加載速度可能會(huì)比較慢。
查看genomes.json文件

genomes.json

查看trackRegistry.json文件
trackRegistry.json

下載并存放個(gè)人的基因組數(shù)據(jù)到data文件中诸老,更改resources文件中的genomes.json文件

cd /var/www/html/igv-webapp/data
# 創(chuàng)建個(gè)人參考基因組文件夾
mkdir hg19 && cd hg19
# 下載參考基因組相關(guān)文件
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后的地址改為自己的服務(wù)器IP端口和數(shù)據(jù)存放的地址


image.png

同樣的隆夯,track的相關(guān)信息(如bam、bed别伏,vcf和gtf等文件)也可以存放在data文件夾中蹄衷,添加或更改tracks文件夾下對(duì)應(yīng)的json文件。除了在服務(wù)器中存放相關(guān)的數(shù)據(jù)厘肮,IGV browser還可以直接從本地愧口,遠(yuǎn)端服務(wù)器或鏈接地址導(dǎo)入相關(guān)數(shù)據(jù)進(jìn)行可視化展示。


image.png

image.png
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末类茂,一起剝皮案震驚了整個(gè)濱河市耍属,隨后出現(xiàn)的幾起案子托嚣,更是在濱河造成了極大的恐慌,老刑警劉巖厚骗,帶你破解...
    沈念sama閱讀 206,839評(píng)論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件示启,死亡現(xiàn)場離奇詭異,居然都是意外死亡领舰,警方通過查閱死者的電腦和手機(jī)夫嗓,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來提揍,“玉大人啤月,你說我怎么就攤上這事±驮荆” “怎么了谎仲?”我有些...
    開封第一講書人閱讀 153,116評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長刨仑。 經(jīng)常有香客問我郑诺,道長,這世上最難降的妖魔是什么杉武? 我笑而不...
    開封第一講書人閱讀 55,371評(píng)論 1 279
  • 正文 為了忘掉前任辙诞,我火速辦了婚禮,結(jié)果婚禮上轻抱,老公的妹妹穿的比我還像新娘飞涂。我一直安慰自己,他們只是感情好祈搜,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,384評(píng)論 5 374
  • 文/花漫 我一把揭開白布较店。 她就那樣靜靜地躺著,像睡著了一般容燕。 火紅的嫁衣襯著肌膚如雪梁呈。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,111評(píng)論 1 285
  • 那天蘸秘,我揣著相機(jī)與錄音官卡,去河邊找鬼。 笑死醋虏,一個(gè)胖子當(dāng)著我的面吹牛寻咒,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播颈嚼,決...
    沈念sama閱讀 38,416評(píng)論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼毛秘,長吁一口氣:“原來是場噩夢(mèng)啊……” “哼!你這毒婦竟也來了粘舟?” 一聲冷哼從身側(cè)響起熔脂,我...
    開封第一講書人閱讀 37,053評(píng)論 0 259
  • 序言:老撾萬榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎柑肴,沒想到半個(gè)月后霞揉,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,558評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡晰骑,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,007評(píng)論 2 325
  • 正文 我和宋清朗相戀三年适秩,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片硕舆。...
    茶點(diǎn)故事閱讀 38,117評(píng)論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡秽荞,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出抚官,到底是詐尸還是另有隱情扬跋,我是刑警寧澤,帶...
    沈念sama閱讀 33,756評(píng)論 4 324
  • 正文 年R本政府宣布凌节,位于F島的核電站钦听,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏倍奢。R本人自食惡果不足惜朴上,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,324評(píng)論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望卒煞。 院中可真熱鬧痪宰,春花似錦、人聲如沸畔裕。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽柴钻。三九已至淮韭,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間贴届,已是汗流浹背靠粪。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評(píng)論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留毫蚓,地道東北人占键。 一個(gè)月前我還...
    沈念sama閱讀 45,578評(píng)論 2 355
  • 正文 我出身青樓,卻偏偏與公主長得像元潘,于是被迫代替她去往敵國和親畔乙。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,877評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容

  • ORA-00001: 違反唯一約束條件 (.) 錯(cuò)誤說明:當(dāng)在唯一索引所對(duì)應(yīng)的列上鍵入重復(fù)值時(shí)翩概,會(huì)觸發(fā)此異常牲距。 O...
    我想起個(gè)好名字閱讀 5,182評(píng)論 0 9
  • 國家電網(wǎng)公司企業(yè)標(biāo)準(zhǔn)(Q/GDW)- 面向?qū)ο蟮挠秒娦畔?shù)據(jù)交換協(xié)議 - 報(bào)批稿:20170802 前言: 排版 ...
    庭說閱讀 10,869評(píng)論 6 13
  • 現(xiàn)在的自己到底算是在經(jīng)歷什么 應(yīng)該說現(xiàn)在的自己一直是活在自己的恐懼當(dāng)中 仔細(xì)想來的話應(yīng)該是自己對(duì)自己沒有自信 你害...
    百花梨落閱讀 307評(píng)論 2 3
  • 廣州易常青企業(yè)宣傳片
    賣報(bào)小話家閱讀 74評(píng)論 0 0
  • 今天是冬至節(jié)返咱,要不是各大群的餃子湯圓紅包雨,我還真忘記了這么一個(gè)節(jié)日牍鞠。 群也許是現(xiàn)在跟生活聯(lián)系最緊密的一部分了咖摹,...
    悠然花開73閱讀 141評(píng)論 3 1