小男孩‘自慰网亚洲一区二区,亚洲一级在线播放毛片,亚洲中文字幕av每天更新,黄aⅴ永久免费无码,91成人午夜在线精品,色网站免费在线观看,亚洲欧洲wwwww在线观看

分享

IGV基因組瀏覽器可視化高通量測序數(shù)據(jù)

 yjt2004us 2018-06-22

IGV是本地瀏覽測序數(shù)據(jù)功能最為強(qiáng)大的基因組瀏覽器,支持多種不同類型的輸入格式和不同的顯示方式,如峰圖、線圖、柱狀圖、Sashimi-plot。同時(shí)還可以配合bedtools使用。

新版的Windows IGV分發(fā)包中包含了Java運(yùn)行環(huán)境,點(diǎn)擊鏈接http://data./igv/projects/downloads/IGV_Win_2.3.93.zip即可下載,解壓后,點(diǎn)擊igv.bat文件即可啟動(dòng)。

若啟動(dòng)失敗,使用記事本打開并編輯igv.bat文件,在文件的最后新起一行輸入pause,保存后,再嘗試打開,就可以在Windows下的命令行界面(黑色背景的小框中)看到錯(cuò)誤信息,根據(jù)信息提示去解決問題。



IGV具體使用如下:



如果基因組未存在于當(dāng)前列表中,則選擇更多從Broad IGV獲取對應(yīng)的基因組信息。若所研究的物種沒有被Broad IGV收錄,則需要自己構(gòu)建基因組信息。具體參見前面發(fā)過的文章 (在公眾號后臺輸入 IGV基因組索引 即可獲?。?。



IGV展示所用的bw文件可以直接從bam文件得來。下面列出幾個(gè)工具和相應(yīng)的代碼輔助轉(zhuǎn)換??紤]到不同樣品的測序深度不同,需要根據(jù)樣品的測序量做一個(gè)標(biāo)準(zhǔn)化處理,建議使用deepToolsbamCoverage命令。



IGV是一個(gè)界面可操作軟件,學(xué)習(xí)的方式就是每個(gè)菜單點(diǎn)開看看,點(diǎn)選下看下效果,比如修改下Track的顏色,Track的高度,縱軸的坐標(biāo)范圍Auto scale, Log scale 或 Data range等方便比較。



如果想比較多個(gè)基因的表達(dá)峰圖,可以先在IGV中添加相應(yīng)的基因列表,然后針對特定的列表選擇view,就可以在一個(gè)屏幕顯示所有加入的基因。



如果我們要查詢的區(qū)域很多,可以配合bedtools igv -i input.bed,生成一個(gè)腳本,批量捕獲多個(gè)區(qū)域的比對信息或coverage信息。使用Tools - Run Batch Script - 選擇bedtools igv的輸出結(jié)果,即可以看到動(dòng)畫般的截屏界面。完成后,就可以一張張的瀏覽圖片或者傳給沒有安裝IGV的老師或朋友查看。



構(gòu)建自己的基因組索引

  1. 在IGV的菜單欄點(diǎn)擊 Genomes>Create .genome File,會(huì)彈出一個(gè)窗口填寫需要的信息。

  2. 輸入一個(gè)ID和基因組的描述信息,這個(gè)ID等同于常見的mm10,會(huì)作為構(gòu)建的基因組索引的名字。

  3. 使用 samtools faidx genome.fasta索引基因組文件,然后輸入基因組FASTA格式文件的路徑。

  4. 提供基因注釋的GTF文件,這個(gè)文件可以是拼裝基因組預(yù)測的基因,也可以是cufflinks 等拼裝的轉(zhuǎn)錄組的GTF文件。

  5. 其它的選項(xiàng)可以忽略。

  6. 點(diǎn)擊 Save,選擇基因組索引文件*.genome的存儲(chǔ)路徑。


    本站是提供個(gè)人知識管理的網(wǎng)絡(luò)存儲(chǔ)空間,所有內(nèi)容均由用戶發(fā)布,不代表本站觀點(diǎn)。請注意甄別內(nèi)容中的聯(lián)系方式、誘導(dǎo)購買等信息,謹(jǐn)防詐騙。如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容,請點(diǎn)擊一鍵舉報(bào)。
    轉(zhuǎn)藏 分享 獻(xiàn)花(0

    0條評論

    發(fā)表

    請遵守用戶 評論公約

    類似文章 更多