ChIP-Seq數(shù)據(jù)挖掘系列-4: liftOver - 基因組坐標(biāo)在不同基因組注釋版本間轉(zhuǎn)換

UCSC Genome Browser

liftover 是UCSC提供的一個(gè)用于不同基因組注釋版本間基因組坐標(biāo)轉(zhuǎn)換的工具优妙。UCSC 提供了一個(gè)web版本程剥;數(shù)據(jù)過多時(shí)用起來很不方便,也可以下載程序到本地進(jìn)行批量操作而芥。

應(yīng)用: 已經(jīng)發(fā)表的chipseq peak 結(jié)果可能基于不同基因組版本可都,這樣就不好直接比較某些peak之間的關(guān)系了,重新分析一次數(shù)據(jù)需要時(shí)間成本婚陪;因此族沃,使用liftover將peak位置信息轉(zhuǎn)換到同一基因組版本不失為為一個(gè)簡單的方法。


#1. 網(wǎng)頁版 liftover在線使用

##1.1 UCSC 官網(wǎng)的web版liftover :

web版liftover

web版liftover

##1.2 使用方法

  1. 選定原始數(shù)據(jù)物種及版本信息
  2. 選定需要轉(zhuǎn)化的目標(biāo)物種及版本信息
  3. 上傳數(shù)據(jù)或直接黏貼數(shù)據(jù)泌参,提交

##1.3 示例

chr1    3505060 3505170
chr1    4504150 4504258
chr1    4791886 4792031
chr1    4792358 4792452
chr1    4797825 4797956
chr1    7079047 7079164
chr1    9550473 9550602
chr1    9615477 9615630
chr1    9616792 9616994
chr1    9623913 9624108
chr1    9624276 9624419
chr1    9685842 9685935
chr1    9929393 9929495
chr1    9932939 9933071
chr1    9974265 9974358

##1.4 示例操作結(jié)果和參數(shù)使用

Result and Parameters Used

#2 命令行版本liftover

##2.1 根據(jù)自己系統(tǒng)選擇程序版本

版本選擇
[linux.x86_64 liftOver]http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/liftOver

##2.2 使用方法

$ chmod +x liftOver
$ liftOver input.bed map.chain output.bed unMapped.bed
#input.bed: 輸入數(shù)據(jù)
#output.bed:轉(zhuǎn)換成功的坐標(biāo)
#unMapped.bed:轉(zhuǎn)換失敗的坐標(biāo)

##2.23 map.chain

? liftOver工作需要的一個(gè)轉(zhuǎn)換數(shù)據(jù)脆淹,此數(shù)據(jù)可以從UCSC下載;

  • map.chain數(shù)據(jù)下載沽一,以mm10ToHg19.over.chain 為例:
    mm10ToHg19.over.chain.gz包含了將mm10 (Mouse GRCm38) 坐標(biāo)轉(zhuǎn)換成 hg19 (Human GRCh37)坐標(biāo)需要的信息盖溺。


ChIP-Seq 數(shù)據(jù)挖掘系列文章目錄:
ChIP-Seq數(shù)據(jù)挖掘系列-1:Motif 分析(1)-HOMER 安裝
ChIP-Seq數(shù)據(jù)挖掘系列-2: Motif 分析(2) - HOMER Motif 分析基本步驟
ChIP-Seq數(shù)據(jù)挖掘系列-3: Motif 分析(3) - 利用ChIP-Seq結(jié)果在基因組區(qū)域中尋找富集的Motifs
ChIP-Seq數(shù)據(jù)挖掘系列-4: liftOver - 基因組坐標(biāo)在不同基因組注釋版本間轉(zhuǎn)換
ChIP-Seq數(shù)據(jù)挖掘系列-5.1: ngs.plot 可視化ChIP-Seq 數(shù)據(jù)
ChIP-Seq數(shù)據(jù)挖掘系列-5.2: ngs.plot 畫圖工具ngs.plot.r 和 replot.r 參數(shù)詳解

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市铣缠,隨后出現(xiàn)的幾起案子烘嘱,更是在濱河造成了極大的恐慌,老刑警劉巖蝗蛙,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件蝇庭,死亡現(xiàn)場離奇詭異,居然都是意外死亡歼郭,警方通過查閱死者的電腦和手機(jī)遗契,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來病曾,“玉大人牍蜂,你說我怎么就攤上這事√┩浚” “怎么了鲫竞?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長逼蒙。 經(jīng)常有香客問我从绘,道長,這世上最難降的妖魔是什么是牢? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任僵井,我火速辦了婚禮,結(jié)果婚禮上驳棱,老公的妹妹穿的比我還像新娘批什。我一直安慰自己,他們只是感情好社搅,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布驻债。 她就那樣靜靜地躺著乳规,像睡著了一般。 火紅的嫁衣襯著肌膚如雪合呐。 梳的紋絲不亂的頭發(fā)上暮的,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天,我揣著相機(jī)與錄音淌实,去河邊找鬼冻辩。 笑死,一個(gè)胖子當(dāng)著我的面吹牛翩伪,可吹牛的內(nèi)容都是我干的微猖。 我是一名探鬼主播谈息,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼缘屹,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了侠仇?” 一聲冷哼從身側(cè)響起轻姿,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎逻炊,沒想到半個(gè)月后互亮,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡余素,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年豹休,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片桨吊。...
    茶點(diǎn)故事閱讀 38,117評論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡威根,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出视乐,到底是詐尸還是另有隱情洛搀,我是刑警寧澤,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布佑淀,位于F島的核電站留美,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏伸刃。R本人自食惡果不足惜谎砾,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望捧颅。 院中可真熱鬧景图,春花似錦、人聲如沸隘道。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至忘晤,卻和暖如春宛蚓,著一層夾襖步出監(jiān)牢的瞬間拼缝,已是汗流浹背谅辣。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留带污,地道東北人闰蛔。 一個(gè)月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓痕钢,卻偏偏與公主長得像,于是被迫代替她去往敵國和親序六。 傳聞我的和親對象是個(gè)殘疾皇子任连,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,877評論 2 345