ucsc ``liftover``基因组坐标转换的纯python实现。

pyliftover的Python项目详细描述


https://travis-ci.org/konstantint/pyliftover.png?branch=master:target:https://travis-ci.org/konstantint/pyliftover

Pyliftover是一个库,用于在不同程序集之间快速轻松地转换基因组(点)坐标。

它使用与ucscliftOver tool相同的逻辑和坐标转换映射。

从当前版本(0.2)开始,pyliftover只进行点坐标的转换,即, 与liftOver不同,它不转换范围,也不提供任何处理床文件的特殊工具。 对于单点坐标,它会产生与liftOver完全相同的输出(目前至少使用hg17ToHg18.over.chain.gz文件进行了验证)。

安装

安装包的最简单方法是通过easy_installpip

$ easy_install pyliftover

用法

库支持的主要用法示例如下:

from pyliftover import LiftOver
lo = LiftOver('hg17', 'hg18')
lo.convert_coordinate('chr1', 1000000)

第一行将自动从ucsc下载hg17-to-hg18坐标转换chain file, 除非它已经缓存或在当前目录中可用。或者,您可以提供自己的链文件:

lo = LiftOver('hg17ToHg18.over.chain.gz')
lo.convert_coordinate('chr1', 1000000, '-')

lo.convert_coordinate调用的结果是None(如果无法识别源染色体名称)或 新组件。列表可以是空的(在新程序集中删除了轨迹)、具有单个元素(轨迹唯一匹配),或者原则上, 具有多个元素(尽管这对于大多数默认的种内基因组转换来说可能是罕见的)。 注意工具中的坐标是基于0的。也就是说,你可以通过^{tt8}在基因组浏览器中引用一个位置$ 对应于pyliftover术语中的坐标9

尽管您可以尝试使用任意链文件(如原始的liftOver工具)应用该工具,但对于 同一物种的不同集合之间的坐标。

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java我的应用程序崩溃,试图添加菜单   java如何创建TBase类实例   如何在java小程序中使用JMF功能?   数据结构Java中是否有用于双循环链表的内置接口?   java在使用JAXB读取XML文件时如何使用继承   在Java中获取和存储子进程的输出   从Java执行Windows命令   HQL中带日期的java限制   安卓如何使用GSON将单个JSON字段映射到多个JAVA字段?   std:random shuffle()和Collections之间的java差异。洗牌   Java正则表达式更高效的动态正则表达式   来自JSON的java POJO表示   active directory无法在Java中更改AD服务器中的密码   DataOutputStream出现java错误   scala如何指定JVM最大堆大小“Xmx”,以便在SBT中运行具有“run”操作的应用程序?   添加了`<f:actionListener>`之后,java`ViewScoped` bean会被重新创建吗?   python找不到满足java要求的版本(来自版本:),找不到与java匹配的发行版   java为什么刷新jdbc swing表中的按钮后数据会多次出现?   java在编程语言(尤其是dart)中,死代码是什么意思?   java为什么我必须在JDBC中使用回滚?