How to convert GTF format into BED12 or BIGBED format?

#see below for UPDATES that include more shorter ways of conversions

#How to convert GTF format into BED12 format (Human-hg19)?

#How to convert GTF or BED format into BIGBED format?

#Why BIGBED (If GTF or BED file is very large to upload in UCSC, you can use trackHubs. However trackHubs do not accept either of the formats. Therefore you would need bigBed format)

#First, download UCSC scripts

wget http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/gtfToGenePred

wget http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/genePredToBed

wget http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bedToBigBed

#Second, download chromosome sizes and filter out unnecessary chromosomes

wget http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/hg19.chrom.sizes

grep -v chrM hg19.chrom.sizes|grep -v _hap|grep -v Un_gl|grep -v random>hg19.chrom.filtered.sizes

rm hg19.chrom.sizes

#Third, make them executable

chmod +x gtfToGenePred genePredToBed bedToBigBed

#Convert Gtf to genePred

./gtfToGenePred 1st_53_tissues.combined.gtf 1st_53_tissues.combined.genePred

#Convert genPred to bed12

./genePredToBed 1st_53_tissues.combined.genePred 1st_53_tissues.combined.bed12

#sort bed12

sort -k1,1 -k2,2n 1st_53_tissues.combined.bed12>1st_53_tissues.combined.sorted.bed

#Convert sorted bed12 to bigBed (useful for trackhubs)

./bedToBigBed 1st_53_tissues.combined.sorted.bed hg19.chrom.filtered.sizes 1st_53_tissues.combined.bb

?著作權歸作者所有,轉載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末测摔,一起剝皮案震驚了整個濱河市置济,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌锋八,老刑警劉巖浙于,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異挟纱,居然都是意外死亡羞酗,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進店門紊服,熙熙樓的掌柜王于貴愁眉苦臉地迎上來檀轨,“玉大人,你說我怎么就攤上這事欺嗤〔翁眩” “怎么了?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵煎饼,是天一觀的道長讹挎。 經(jīng)常有香客問我,道長吆玖,這世上最難降的妖魔是什么筒溃? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任,我火速辦了婚禮沾乘,結果婚禮上怜奖,老公的妹妹穿的比我還像新娘。我一直安慰自己翅阵,他們只是感情好歪玲,可當我...
    茶點故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著掷匠,像睡著了一般读慎。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上槐雾,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天,我揣著相機與錄音幅狮,去河邊找鬼募强。 笑死株灸,一個胖子當著我的面吹牛,可吹牛的內(nèi)容都是我干的擎值。 我是一名探鬼主播慌烧,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼鸠儿!你這毒婦竟也來了屹蚊?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤进每,失蹤者是張志新(化名)和其女友劉穎汹粤,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體田晚,經(jīng)...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡嘱兼,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了贤徒。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片芹壕。...
    茶點故事閱讀 38,117評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖接奈,靈堂內(nèi)的尸體忽然破棺而出踢涌,到底是詐尸還是另有隱情,我是刑警寧澤序宦,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布睁壁,位于F島的核電站,受9級特大地震影響挨厚,放射性物質(zhì)發(fā)生泄漏堡僻。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一疫剃、第九天 我趴在偏房一處隱蔽的房頂上張望钉疫。 院中可真熱鬧,春花似錦巢价、人聲如沸牲阁。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽城菊。三九已至,卻和暖如春碉克,著一層夾襖步出監(jiān)牢的瞬間凌唬,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工漏麦, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留客税,地道東北人况褪。 一個月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓,卻偏偏與公主長得像更耻,于是被迫代替她去往敵國和親测垛。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,877評論 2 345

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

  • 本文資料來源 https://web.archive.org/web/20161125133249/http://...
    x2yline閱讀 2,772評論 0 4
  • ANNOVAR的安裝 ANNOVAR網(wǎng)址 log in之后才能download秧均,使用教育機構后綴的郵箱即可注冊食侮。 ...
    面面的徐爺閱讀 22,890評論 1 26
  • rljs by sennchi Timeline of History Part One The Cognitiv...
    sennchi閱讀 7,294評論 0 10
  • 穿過密林達頂端,景物濃縮眼底寬目胡。 汽車大小一厘米锯七,樓群各異積木盤。
    老槐樹閱讀 148評論 0 2
  • 禮貌是指塑造人際關系的能力讶隐。禮貌的表現(xiàn)形式是言行舉止起胰。一個人通過自己的言行舉止表現(xiàn)他對社會行為規(guī)范的承認,表現(xiàn)他對...
    寇廷聚閱讀 333評論 0 2