这篇文章主要介绍MCscan如何做基因组共线性分析，文中介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们一定要看完！
多基因组基因水平共线性分析
在 多基因间共线性及局部基因水平的共线性 和 基因组间共线性分析 里讲了MCscan如何做基因组共线性分析，这里再说说多基因组基因水平的共线性分析。
方法是两两基因组做共线性分析。还是以桃子，葡萄和可可为例，方法如下:
$ python -m jcvi.compara.synteny mcscan grape.bed grape.peach.lifted.anchors --iter=1 -o grape.peach.i1.blocks $ python -m jcvi.compara.synteny mcscan grape.bed grape.cacao.lifted.anchors --iter=1 -o grape.cacao.i1.blocks $ python -m jcvi.formats.base join grape.peach.i1.blocks grape.cacao.i1.blocks --noheader | cut -f1,2,4,6 > grape.blocks
筛选关注的区域进行展示，这里展示前50列：
$ head -50 grape.blocks > blocks2
blocks2文件格式如下：
GSVIVT01012261001 . . GSVIVT01012259001 ppa005716m . GSVIVT01012258001 . . GSVIVT01012257001 ppa002846m . GSVIVT01012255001 ppa000919m Thecc1EG011472t1 GSVIVT01012253001 ppa010733m Thecc1EG011473t1 GSVIVT01012252001 ppa015194m Thecc1EG011474t1
配置文件如下：
# x, y, rotation, ha, va, color, ratio, label 0.5, 0.6, 0, center, top, , 1, grape Chr1 0.3, 0.4, 0, center, bottom, , .5, peach scaffold_1 0.7, 0.4, 0, center, bottom, , .5, cacao scaffold_2 # edges e, 0, 1 e, 0, 2
运行画图：
$ cat grape.bed peach.bed cacao.bed > grape_peach_cacao.bed $ python -m jcvi.graphics.synteny blocks2 grape_peach_cacao.bed blocks2.layout
以上是“MCscan如何做基因组共线性分析”这篇文章的所有内容，感谢各位的阅读！希望分享的内容对大家有帮助，更多相关知识，欢迎关注亿速云行业资讯频道！
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
