|
整合基因組瀏覽器(IGV)是一種高性能的可視化工具,用來交互式地探索大型綜合基因組數(shù)據(jù)。它支持各種數(shù)據(jù)類型,包括array-based的和下一代測序的數(shù)據(jù)和基因注釋。IGV這個(gè)工具很牛,發(fā)了NB: 使用IGV的基本思路安裝和配置要用linux上的IGV(首先需要去官網(wǎng)下載一個(gè)Linux下的IGV),然后需要一個(gè)遠(yuǎn)程控制工具: 這里使用 VNC,網(wǎng)上下載破解版,安裝(輸入注冊碼即可破解,文章底部有注冊碼)。 運(yùn)行IGV開啟 linux 的 VNCserver,開啟了之后,linux上會顯示節(jié)點(diǎn)和端口,用于Windows端的連接。 vncserver -geometry 1366x768 #調(diào)節(jié)分辨率
運(yùn)行 VNC 客戶端(VNC viewer)里的 IGV.sh 開啟可視化程序(打開Windows的VNC會自動跳轉(zhuǎn)到Linux下的IGV安裝目錄,雙擊IGV.sh即可打開)。 1.載入?yún)⒖蓟蚪M數(shù)據(jù) 可以點(diǎn)擊選擇一個(gè),也可以通過菜單欄Genomes載入文件。 2.載入排序后的bam文件 通過菜單欄的File,load from file。。。 IGV可視化數(shù)據(jù)準(zhǔn)備不是什么數(shù)據(jù)都可以拿IGV看的,參考基因組必須為FASTA格式; IGV只是負(fù)責(zé)將比對結(jié)果可視化,并沒有比對過程,所以不能直接載入reads; 需要將待比對的reads與前面指定的參考基因組用bwa進(jìn)行比對; 比對后的sam文件也不能直接載入(麻煩),要轉(zhuǎn)bam;bam排序;bam建索引(可以一步完成); 命令流程如下:BWA出來的sam文件無法直接被 IGV 利用的,必須經(jīng)過 samtools 處理后才能被 IGV 顯示出來。 samtools view -@ 10 -bS -F 4 ./contigs_sequence_align_to_public_genome.sam > ./contigs_sequence_align_to_public_genome.bam samtools sort -@ 10 ./contigs_sequence_align_to_public_genome.bam contigs_sequence_align_to_public_genome.sorted samtools index contigs_sequence_align_to_public_genome.sorted.bam contigs_sequence_align_to_public_genome.sorted.bai
samtools depth contigs_sequence_align_to_public_genome.sorted.bam >depth_reads.txt wc -l depth_reads.txt > Coverage-aln_reads.txt更多的使用細(xì)節(jié)參照 samtools 使用方法。 IGV 需要設(shè)置內(nèi)存,在看大型基因組時(shí)很吃內(nèi)存. 示例: 下面顯示的是組裝異常區(qū),這是手動連接的地方,Pacbio全部不能順利跨過去,說明組裝有問題。 以下是單細(xì)胞轉(zhuǎn)錄組比對可視化:
IGV基本操作如何去掉右邊自動彈出的框框? 藍(lán)色框起來的reads是啥? …
IGV用戶指南目錄:
1.用戶界面1.1 Main Window(主窗口)
1.tool bar(工具欄),menu bar(菜單欄),pop-up menus(彈出式菜單)。 2.染色體上的紅色盒子表示顯示這部分染色體,顯示完整染色體是紅框會消失。 3.尺度顯示了染色體的可見部分,刻度線顯示了染色體的位置,跨度列表顯示了當(dāng)前顯示的堿基的數(shù)量。 4.IGV在水平行顯示的數(shù)據(jù)稱為tracks。通常,每個(gè)tracks代表一個(gè)樣本或?qū)嶒?yàn)。這個(gè)例子展示了甲基化、基因表達(dá)、拷貝數(shù),LOH和突變數(shù)據(jù)。 5.IGV也顯示某些特性,比如在tracks中的基因。默認(rèn)情況下,IGV在一個(gè)面板顯示數(shù)據(jù),在另一個(gè)面板顯示數(shù)據(jù)特性。拖放一個(gè)track名稱,將一個(gè)track從一個(gè)面板移動到另一個(gè)地方。 6.Track名稱列在最左邊面板。名字的易讀性取決于 tracks的高度,例如,track越小,它的名字的可讀性越小。 7.屬性名稱被列在頂部的屬性面板。彩色塊代表屬性值,每個(gè)獨(dú)特的值被都有一個(gè)獨(dú)特的顏色。鼠標(biāo)放在一個(gè)顏色塊的附近來查看其屬性值。
參考資料Viewing Alignments(查看比對,核心文檔) IGV軟件使用指南(博客)
軟件安裝: RealVNC是一款免費(fèi)的遠(yuǎn)程控制程序,它是VNC (Virtual Network Computing)眾多操作平臺版本中的一員,這里小編提供的RealVNC適用于windows操作平臺,軟件采用2048 位 RSA 服務(wù)器驗(yàn)證和 128 位 AES 會話加密術(shù),安裝此軟件的電腦可在全球各地遠(yuǎn)程操控自己的電腦,所操作平臺沒有任何限制。 |
|
|