IGV-基因組瀏覽器-改造記錄(一)

我一直都相信啃炸,
只要有足夠長(zhǎng)的時(shí)間崭放,沒(méi)有什么事情是我們做不到的。
我也明白担平,
人的生命是有限的李剖,所以追求的東西應(yīng)該盡可能少芒率。
我只是沒(méi)想到,
很多時(shí)候篙顺,不是你沒(méi)時(shí)間偶芍,而是你身邊的人,等不了德玫。而事實(shí)上匪蟀,你也等不了。
感傷是沒(méi)有意義的宰僧。
無(wú)論何時(shí)何地萄窜,保持往前行進(jìn),總是對(duì)的。

寫(xiě)在前面

在前述的推文查刻,我已經(jīng)分享了《任何人键兜,如何改造IGV基因組瀏覽器》。但是似乎穗泵,并沒(méi)有太多人感興趣普气。當(dāng)然,這個(gè)也并不值得Care佃延,畢竟不是每個(gè)人都有時(shí)間做一個(gè)自己需要的工具现诀。
前述我已經(jīng)提過(guò),設(shè)計(jì)履肃,實(shí)現(xiàn)了兩個(gè)我自己其實(shí)并不滿意的瀏覽器仔沿,所以我最終還是決定直接改造IGV。
浪費(fèi)了兩天時(shí)間尺棋,終于基本可以認(rèn)為是封锉,完成了一個(gè)兩年前就計(jì)劃做的事情
這個(gè)事情膘螟,也是現(xiàn)在的博導(dǎo)與我剛認(rèn)識(shí)不久成福,就討論到的一個(gè)事情。
可能時(shí)間確實(shí)能解決以前荆残。兩年前奴艾,做不到的事情;現(xiàn)在并沒(méi)有覺(jué)得太過(guò)復(fù)雜就折騰好了内斯。

改造IGV的初衷

關(guān)注過(guò)sRNA測(cè)序數(shù)據(jù)分析蕴潦,尤其是相對(duì)高水平的植物小RNA相關(guān)論文,那么應(yīng)該見(jiàn)過(guò)類似下圖俘闯。



這三個(gè)小圖潭苞,展示的都是sRNA中特殊的一類,即phasiRNAs(往往與育性有關(guān)备徐,如水稻的熱敏雄性不育)萄传。三個(gè)圖甚颂,均來(lái)自植物sRNA研究大牛Blake課題組的網(wǎng)頁(yè)瀏覽器(傳言蜜猾,瀏覽數(shù)據(jù)和輸出可發(fā)表的圖片,需要的是兩個(gè)不同版本的瀏覽器)振诬。呃....各種原因蹭睡,我從來(lái)沒(méi)有使用過(guò)...

兩年前,博導(dǎo)與我商量赶么,猜想著肩豁,是否有可能,在IGV上直接展示這類點(diǎn)圖,那么一方面高效清钥,另一方面直接得到矢量圖琼锋。那會(huì)他的設(shè)想是改造IGV,我的設(shè)想是自己從頭寫(xiě)一個(gè)祟昭。中間的插曲是缕坎,我花了大半天,確實(shí)是改造了一個(gè)看起來(lái)能用的IGV篡悟,其實(shí)有問(wèn)題谜叹。后續(xù),我們都沒(méi)有再提起這個(gè)事情搬葬。但這個(gè)想法荷腊,我從來(lái)就沒(méi)有放下過(guò)。為此急凰,我在寫(xiě)JIGplot圖形引擎的時(shí)候女仰,都是不時(shí)地會(huì)考慮一些可以用在基因組瀏覽器的部件。
總的來(lái)說(shuō)香府,兩戰(zhàn)兩敗董栽,無(wú)論我怎么寫(xiě),再有看起來(lái)不錯(cuò)的特性企孩,始終感覺(jué)無(wú)法超過(guò)IGV锭碳。只能認(rèn)慫,重新回來(lái)改造IGV勿璃。于是有了以下功能擒抛。

sRNA測(cè)序數(shù)據(jù)的readLen點(diǎn)圖

在IGV中,每一個(gè)BAM或者SAM文件導(dǎo)入之后补疑,可以實(shí)時(shí)生成一個(gè)Coverage Track歧沪,如果是RNAseq數(shù)據(jù)的話,那么也還可以實(shí)時(shí)生成一個(gè)Juncitons Track莲组。這兩類Track在大多數(shù)數(shù)據(jù)分析中都是有用的诊胞,如


可以很好的看出每個(gè)位點(diǎn)的覆蓋率,可變剪接等锹杈。
但是撵孤,對(duì)于小RNA數(shù)據(jù),則并沒(méi)有太多有用的插件竭望。尤其是當(dāng)我們需要前述Blake課題組的瀏覽器類似的功能時(shí)邪码。于是,改之可得咬清。

改造闭专,并使IGV對(duì)于BAM文件的展示奴潘,即AlignmentTrack可以跟支持CoverageTrack一樣,支持一種新的Track影钉。
其中画髓,每一個(gè)位置會(huì)有0~5個(gè)點(diǎn),不同點(diǎn)的顏色對(duì)應(yīng)不同的read長(zhǎng)度平委。如橙色對(duì)應(yīng)的是24nt的sRNA雀扶,而淺藍(lán)色對(duì)應(yīng)的是21nt的sRNA,綠色對(duì)應(yīng)的是22nt的sRNA肆汹,紫色對(duì)應(yīng)23nt的sRNAs愚墓,其他長(zhǎng)度對(duì)應(yīng)的是剩余長(zhǎng)度的sRNAs。點(diǎn)的高度對(duì)應(yīng)的是該sRNAreads的豐度昂勉。
所以整體上浪册,可以辨識(shí)出,右側(cè)即為一個(gè)21nt的PHAS locus岗照。

由于村象,我們?cè)黾拥囊粋€(gè)新的Track類型,所以整體上攒至,他的使用與CoverageTrack是類似的厚者。于是可以在不查看Alignment的時(shí)候,同時(shí)可視化多個(gè)BAM文件迫吐,進(jìn)行不同樣本之間的比較库菲。如



可以看出,第一個(gè)樣品右側(cè)的locus上24nt的reads相對(duì)于其他樣品多志膀,或許是甲基化高頻區(qū)域熙宇。

功能的調(diào)出

開(kāi)啟這個(gè)Track的方式,與Coverage的開(kāi)啟和關(guān)閉是類似的溉浙,


另外的烫止,當(dāng)然也只是AutoScale(上圖均用的SetDataRange...)

寫(xiě)在后面

事實(shí)上,也有不少其他增加的功能戳稽。留待后續(xù)再推馆蠕。
距離人生大事,有十天時(shí)間惊奇』ス總的來(lái)說(shuō),整體上這個(gè)事情赊时,麻煩了雙方家人較多事情吨铸,后續(xù)也要麻煩很多的親朋好友行拢。
這幾天邀請(qǐng)了一些朋友祖秒,也忘了邀請(qǐng)一些朋友,也不敢邀請(qǐng)一些朋友。
這似乎是一個(gè)比較復(fù)雜的事情竭缝。不邀請(qǐng)的房维,或許會(huì)擔(dān)心朋友覺(jué)得你沒(méi)把他放在心上;邀請(qǐng)的抬纸,則會(huì)擔(dān)心朋友不好意思拒絕而不得不(請(qǐng)假并趕去)參與咙俩。我最后是更多的選擇不邀請(qǐng)。
當(dāng)然湿故,這個(gè)過(guò)程中阿趁,我也重新審視了自己過(guò)去的時(shí)間里面,與認(rèn)識(shí)的人之間的交往方式坛猪。事實(shí)上脖阵,我發(fā)現(xiàn),我的朋友確實(shí)很少墅茉;尤其是碩士以前命黔,這主要是我個(gè)人始終沒(méi)搞明白,到底人與人之間的交往就斤,如果才是正確悍募。或者直白一些說(shuō)洋机,我似乎很少真誠(chéng)地對(duì)待確實(shí)將我視作朋友的人坠宴。
碩士以后,我認(rèn)識(shí)了現(xiàn)在的女朋友绷旗,我才真正明白啄踊,其實(shí)每個(gè)人都可以過(guò)得很真實(shí),并沒(méi)有必要去擔(dān)心太多其他人的想法和看法刁标〉咄ǎ或者簡(jiǎn)單地說(shuō),還是要活出自己膀懈。更或者通俗地說(shuō)顿锰,“人活著,最重要的是開(kāi)心”启搂。
總的來(lái)說(shuō)硼控,或許我認(rèn)知內(nèi),朋友之間的交往胳赌,應(yīng)該是一個(gè)真心換真心的過(guò)程牢撼。而是否真誠(chéng),其實(shí)一兩次交談便可以知曉疑苫。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末熏版,一起剝皮案震驚了整個(gè)濱河市纷责,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌撼短,老刑警劉巖再膳,帶你破解...
    沈念sama閱讀 218,122評(píng)論 6 505
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異曲横,居然都是意外死亡喂柒,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,070評(píng)論 3 395
  • 文/潘曉璐 我一進(jìn)店門(mén)禾嫉,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)灾杰,“玉大人,你說(shuō)我怎么就攤上這事熙参】月叮” “怎么了?”我有些...
    開(kāi)封第一講書(shū)人閱讀 164,491評(píng)論 0 354
  • 文/不壞的土叔 我叫張陵尊惰,是天一觀的道長(zhǎng)讲竿。 經(jīng)常有香客問(wèn)我,道長(zhǎng)弄屡,這世上最難降的妖魔是什么题禀? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,636評(píng)論 1 293
  • 正文 為了忘掉前任,我火速辦了婚禮膀捷,結(jié)果婚禮上迈嘹,老公的妹妹穿的比我還像新娘。我一直安慰自己全庸,他們只是感情好秀仲,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,676評(píng)論 6 392
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著壶笼,像睡著了一般神僵。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上覆劈,一...
    開(kāi)封第一講書(shū)人閱讀 51,541評(píng)論 1 305
  • 那天保礼,我揣著相機(jī)與錄音,去河邊找鬼责语。 笑死炮障,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的坤候。 我是一名探鬼主播胁赢,決...
    沈念sama閱讀 40,292評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼白筹!你這毒婦竟也來(lái)了智末?” 一聲冷哼從身側(cè)響起谅摄,我...
    開(kāi)封第一講書(shū)人閱讀 39,211評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎吹害,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體虚青,經(jīng)...
    沈念sama閱讀 45,655評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡它呀,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,846評(píng)論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了棒厘。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片纵穿。...
    茶點(diǎn)故事閱讀 39,965評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖奢人,靈堂內(nèi)的尸體忽然破棺而出谓媒,到底是詐尸還是另有隱情,我是刑警寧澤何乎,帶...
    沈念sama閱讀 35,684評(píng)論 5 347
  • 正文 年R本政府宣布句惯,位于F島的核電站,受9級(jí)特大地震影響支救,放射性物質(zhì)發(fā)生泄漏抢野。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,295評(píng)論 3 329
  • 文/蒙蒙 一各墨、第九天 我趴在偏房一處隱蔽的房頂上張望指孤。 院中可真熱鬧,春花似錦贬堵、人聲如沸恃轩。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 31,894評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)叉跛。三九已至,卻和暖如春蒸殿,著一層夾襖步出監(jiān)牢的瞬間昧互,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 33,012評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工伟桅, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留敞掘,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,126評(píng)論 3 370
  • 正文 我出身青樓楣铁,卻偏偏與公主長(zhǎng)得像玖雁,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子盖腕,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,914評(píng)論 2 355

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