基因在染色體上位置可視化

感謝CJ大神的反復(fù)修改,這個功能越變越好了象迎!
本文旨在復(fù)現(xiàn)CJ大神推文內(nèi)容丹擎,具體請參考“省心炸雞”,手動狗頭(http://www.reibang.com/p/b9034ce82ec2

首先需要獲取的是染色體位置刽宪,我的辦法是在GFF3文件里取得。



復(fù)制出來的文件都在同一個單元格里,這個時候只需要使用分列功能就能把信息分開瓜富。




保存為txt文件
Chr Length
1A 594102056
1B 689851870
1D 495453186
2A 780798557
2B 801256715
2D 651852609
3A 750843639
3B 830829764
3D 615552423
4A 744588157
4B 673617499
4D 509857067
5A 709773743
5B 713149757
5D 566080677
6A 618079260
6B 720988478
6D 473592718
7A 736706236
7B 750620385
7D 638686055

之后需要準(zhǔn)備基因位置文件,首先用tbtools提取所有基因的ID和位置降盹。這里使用GFF3 gene position parse与柑。正經(jīng)的IOS(input,output,start)流程沒什么可說的价捧。



之后再input之前獲得的文件丑念。在set condition text填入需要提取的gene ID。模式選擇包含(contain)干旧,過濾模式選擇提惹邸(extract)。



需要注意的是椎眯,提取出來的文件行列不能直接使用挠将,需要設(shè)置table行列。

第一列g(shù)ene或者transcript ID编整,第二列chromosome舔稀,第三列起始位置,第四列結(jié)束為止掌测,第五列圖形分組(從0開始排起内贮,目前只試過四種最多幾種不清楚)。需要注意的是汞斧,該功能里merge會把不同顏色的基因堆疊在一起夜郁,而如果ID來自不同的圖形分組(第五列信息)則不會。
填入對應(yīng)信息即可


默認(rèn)全部分成一組的話粘勒,你就能得到



顯然沒什么意義竞端,分成三組,再加入顏色信息庙睡。



說個題外話事富,關(guān)于顏色的批量填充,在excel里按照組分類可以寫個if套娃乘陪,例如下面這樣
=if(c1=0,"253,180,98",if(c1=1,"228,26,28","31,120,180"))

也可以寫個index统台,然后用match去返回需要的顏色標(biāo)簽(R,G,B)

=INDEX({"253,180,98","228,26,28","31,120,180"},MATCH(C1,{0,1,2},0))

同理,你可以用類似方法給自己的基因批量分類啡邑。
最后的最后贱勃,一定要換形狀,否則merge的時候信息會出錯谣拣。


錯誤的merge

正確的merge.png

哦對了 基因密度信息請參考這里:http://www.reibang.com/p/801807865864
最后募寨,祝大家磕鹽順利

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市森缠,隨后出現(xiàn)的幾起案子拔鹰,更是在濱河造成了極大的恐慌,老刑警劉巖贵涵,帶你破解...
    沈念sama閱讀 218,682評論 6 507
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件列肢,死亡現(xiàn)場離奇詭異恰画,居然都是意外死亡,警方通過查閱死者的電腦和手機瓷马,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,277評論 3 395
  • 文/潘曉璐 我一進店門拴还,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人欧聘,你說我怎么就攤上這事片林。” “怎么了怀骤?”我有些...
    開封第一講書人閱讀 165,083評論 0 355
  • 文/不壞的土叔 我叫張陵费封,是天一觀的道長。 經(jīng)常有香客問我蒋伦,道長弓摘,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,763評論 1 295
  • 正文 為了忘掉前任痕届,我火速辦了婚禮韧献,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘研叫。我一直安慰自己锤窑,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 67,785評論 6 392
  • 文/花漫 我一把揭開白布嚷炉。 她就那樣靜靜地躺著果复,像睡著了一般。 火紅的嫁衣襯著肌膚如雪渤昌。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,624評論 1 305
  • 那天走搁,我揣著相機與錄音独柑,去河邊找鬼。 笑死私植,一個胖子當(dāng)著我的面吹牛忌栅,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播曲稼,決...
    沈念sama閱讀 40,358評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼索绪,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了贫悄?” 一聲冷哼從身側(cè)響起瑞驱,我...
    開封第一講書人閱讀 39,261評論 0 276
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎窄坦,沒想到半個月后唤反,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體凳寺,經(jīng)...
    沈念sama閱讀 45,722評論 1 315
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,900評論 3 336
  • 正文 我和宋清朗相戀三年彤侍,在試婚紗的時候發(fā)現(xiàn)自己被綠了肠缨。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 40,030評論 1 350
  • 序言:一個原本活蹦亂跳的男人離奇死亡盏阶,死狀恐怖晒奕,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情名斟,我是刑警寧澤脑慧,帶...
    沈念sama閱讀 35,737評論 5 346
  • 正文 年R本政府宣布,位于F島的核電站蒸眠,受9級特大地震影響漾橙,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜楞卡,卻給世界環(huán)境...
    茶點故事閱讀 41,360評論 3 330
  • 文/蒙蒙 一霜运、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧蒋腮,春花似錦淘捡、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,941評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至作彤,卻和暖如春膘魄,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背竭讳。 一陣腳步聲響...
    開封第一講書人閱讀 33,057評論 1 270
  • 我被黑心中介騙來泰國打工创葡, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人绢慢。 一個月前我還...
    沈念sama閱讀 48,237評論 3 371
  • 正文 我出身青樓灿渴,卻偏偏與公主長得像,于是被迫代替她去往敵國和親胰舆。 傳聞我的和親對象是個殘疾皇子骚露,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,976評論 2 355

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