廢話不多說烦却,直接簡單粗暴的開始吧
導(dǎo)入基因組
導(dǎo)入基因組文件:Genomes→load geome from file→選擇"reference.fa"文件
導(dǎo)入注釋文件:file →load geome from file →選擇om"reference.gff"辆苔,由于gff文件已經(jīng)是按坐標(biāo)順序排列的,因此不需要重新排序乓序。
導(dǎo)入甲基化數(shù)據(jù)文件
這是拿到的結(jié)果披坏,標(biāo)注了甲基化的類型以及所在位點。
IGV導(dǎo)入的數(shù)據(jù)需要改為IGV格式廉邑,建議先學(xué)習(xí)下IGV文件的格式
可以看到,官網(wǎng)中給了我們IGV文件的格式例子倒谷,patient就可以改為不同的甲基化類型蛛蒙,由于IGV格式的文件不支持正負(fù)鏈的標(biāo)注,所以在轉(zhuǎn)換成IGV格式文件時恨锚,我把負(fù)鏈的甲基化值變?yōu)樨?fù)值
最后修改的IGV文件如下(tab分割):
然后通過file→load from file導(dǎo)入宇驾,結(jié)果如下圖:
接下來倍靡,我們可以修改圖片的顏色和字體大小猴伶,右鍵點擊相應(yīng)的行即可
點擊file→save image即可保存最終圖片
轉(zhuǎn)載請注明:周小釗的博客>>>使用IGV簡單繪制甲基化分布圖