一看就能学会的，不同基因组版本坐标转换方法
欢迎关注”生信修炼手册”!
对于同一个物种而言，会存在不同的基因组组装版本，以human为例，UCSC有以下多个版本
hg38/GRCH38
hg19/GRCH19
hg18/NCBI18
由于版本太多，这里就没有一一列举完，完整的列表可以查看以下链接
http://hgdownload.soe.ucsc.edu/downloads.html#human
在使用基因组浏览器时，必须保证导入的基因组版本和数据所用的基因组版本一致，比如对于一个基于hg19的vcf文件，参考基因组也必须是hg19。在实际使用过程中，会遇到基因组版本不一致的问题，此时就需要进行基因组版本之间的转换，最常用的工具就是UCSC出品的liftover, 该工具既有在线服务，也有命令行版本。
命令行版本集成在UCSC Utilities工具集中，linux版本的链接如下
http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64.v385/liftOver
使用liftover有一个前提，需要从UCSC下载参考基因组对应的over.chain文件，在每个基因组组装的版本中都提供了liftover对应的over.chain文件，以human hg19为例，对应的文件如下
http://hgdownload.soe.ucsc.edu/goldenPath/hg19/liftOver/
除了在不同基因组版本之间转换，还提供了不同基因组之间的转换功能。liftover的基本用法如下
liftOver hg19.bed hg19ToHg38.over.chain.gz map.bed unmap.bed第一个参数为输入的bed文件，格式为bed3, 内容示意如下
第二个参数为对应的over.chain文件，第三个参数为输出的转换之后的bed文件，第四列为转换失败的bed文件。
对应的在线服务网址如下
https://genome.ucsc.edu/cgi-bin/hgLiftOver
在线版本内置了ove.chain文件，只需选择需要转换的基因组版本，提供需要转换的bed坐标文件即可。
·end·
—如果喜欢，快分享给你的朋友们吧—
原创不易，欢迎收藏，点赞，转发！生信知识浩瀚如海，在生信学习的道路上，让我们一起并肩作战！
本公众号深耕耘生信领域多年，具有丰富的数据分析经验，致力于提供真正有价值的数据分析服务，擅长个性化分析，欢迎有需要的老师和同学前来咨询。
更多精彩
KEGG数据库，除了pathway你还知道哪些
全网最完整的circos中文教程
DNA甲基化数据分析专题
突变检测数据分析专题
mRNA数据分析专题
lncRNA数据分析专题
circRNA数据分析专题
miRNA数据分析专题
单细胞转录组数据分析专题
chip_seq数据分析专题
Hi-C数据分析专题
HLA数据分析专题
TCGA肿瘤数据分析专题
基因组组装数据分析专题
CNV数据分析专题
GWAS数据分析专题
2018年推文合集
2019年推文合集
写在最后
转发本文至朋友圈，后台私信截图即可加入生信交流群，和小伙伴一起学习交流。
扫描下方二维码，关注我们，解锁更多精彩内容！
一个只分享干货的
生信公众号
本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！