1.首先在Ensemblefungi獲取該物種整個基因組的gff注釋文件(或gtf注釋文件,有誰就下載誰)
2.指定基因id的獲取
可以用表格形式打開gff文件雕拼,一般指定基因的id命名如下:其實就是你目標基因前面還有一些信息需要被補上,TBtools才能識別曙强。最后整理一下的基因id迟蜜,如后圖。
3.直接利用TBtools的以下功能泳秀,只需要在藍色標記的地方讀入gff文件,和基因id文件榄攀,點擊start即可嗜傅。結(jié)果如后圖。
4.到這里小伙伴們是不是覺得空空的柱形圖不怎么好看檩赢,所以我看了大佬得到基因密度提嚷类帧(http://www.reibang.com/p/801807865864),再把相應(yīng)的文件讀入(藍色標記處)贞瞒,就得到了一張好看的圖偶房。
后記:本人今天看了各位大佬的介紹第一次做成功,希望可以幫到小白军浆。(特別感謝各位大佬的分享的解析棕洋,入行半年,從R語言開始瘾敢,導師僅僅提供大的方向拍冠,很多技術(shù)手段我均從各位大大那里獲取,感謝你們一步步成就我簇抵。)