将人类mtdna序列转化为变异位点,反之亦然。

oldowan.mtconvert的Python项目详细描述


mtconvert是一个小型、纯python的生物信息实用程序,用于(1)将人类有丝分裂软骨dna序列数据转换为相对于修订的剑桥参考序列(rcrs)的变异位点,以及(2)将变异位点数据转换为dna序列。关于人类mtdna序列的rcrs和变异位点的更多信息,请访问mtconvert网站。

最简单的安装方法是使用
setuptools\u包中的“easy\u install”。这通常是这样的::

$easy_install oldowan.mtconvert


或在类unix系统上,假设您是以非特权用户的身份安装到主python
``site packages``目录,则此::

$sudo easy_install oldowan.mtconvert

您还可以使用标准python distutils安装方法。从本页底部的文件列表中下载
当前源文件存档,将其取消存档,然后安装。在Mac OS X和许多其他类似Unix的系统上,下载了存档文件并将其更改到shell中包含此存档文件的目录后,可能会出现以下情况:


$tar xvzf oldowan.mtconvert*
$cd oldowan.mtconvert*
$python setup.py install

quick start
==``以及来自oldowan.mtconvert的“sites2seq”:

>;>;来自oldowan.mtconvert导入seq2sites,sites2seq

>;>;将序列转换为sites::

>;seq=“”ttcttgaggaagttggtaccacca
gtattgactcaccaccaccaccaccacacacacaccctatgttcgtacattactgcc
agccaccatgaatactagtaccatacacacacccattagt
acataaaacacacatccaaaccccccctcataccataccaaagca
agtaccaaacccataccacccaaccataccaaacagt
acatagtataaaagccattacgtacatagcatacatacattacagtcaatc
ccttctcgtcc“
>>>seq2站点(seq)

序列必须连续!必须单独分析序列的单独运行,如hvr1和hvr2
而无中间序列间隔。


序列中允许的周围位点(n)的数量也有一个截止值。默认情况下,这是10-但这是一个可以设置的选项:

>;>seq2站点(seq,ambig_cutoff=20)

返回
的默认序列区域是高变量区域1(hvr1),即rcr的位置16024到16365
(在基于生物的编号中):

>;sites2seq('16129a 16223t')


预定义的序列区域是:

-hvr1:16024-16365
-hvr2:73-340
-hvr1到2:16024-340
-编码:577-15992
-全部:1-16559

so,将hvr2站点列表转换为序列::

>;>;sites2seq('73g',region='hvr2')


站点也可以在列表中提供:

>;>;sites2seq(['16129a','16223t','73g',region='hvr1 to 2')


rcrs序列将在给定空字符串、空列表或
字符串“rcrs”的情况下返回。所有这些都是等效的:

>;>sites2seq(''
>;>sites2seq([])
>;>sites2seq('rcrs')

>可以通过将站点列表传递给
``region``选项来选择任意位置:

>;
python range函数对此很方便,但您必须记住
范围不包括其结束位置::

>;>;sites2seq('',region=range(73341));包括340,但不包括341

2009年)
模块的初始版本。

1.0.1(2009年3月25日)
小版本修复

1.0.2(2009年5月27日)
部分RFLP实现

<1.0.3(2015年6月22日)
在查询结束时添加虚假删除修复


1.0.4(2015年6月22日)
在查询结束时改进虚假删除修复


1.0.5(2015年6月22日)
请求区域外的站点应无提示通过;插入修复

1.0.6(8月4日,2015)
修复版本号安装问题



…_设置工具:http://peak.telecommunity.com/devcenter/easyinstall

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

推荐PyPI第三方库


热门话题
安卓中通过TCP发送时java文件损坏   Java测试预期=异常失败,断言错误   java ssh4 LazyInitializationException   java是否可以在集合中添加重复项?   java是一个免费的开源数据库管理工具   java是否可以在导出的html中嵌入图像   编译器构造解释一些东西,并用Java运行生成的字节码?   java KeyPairGenerator未生成随机密钥   java使用正则表达式生成字符串而不是匹配字符串   java中的多线程线程间通信   具有复合密钥的java Hibernate合并问题   java不能在MainActivity类之外使用SharedReferences变量   正在获取切割/零件异常消息。使用java Spring ControllerAdvice   java tomcat漏洞讨论   javascript在Java中实现“system”命令   java如何在JUnit5*中加载*Spring测试上下文之前获取回调?   java bluej关于缺少返回语句   JavaSpringBoot@ModelAttribute包含使用字符串键的其他模型   从可选中抛出异常。ifPresent()Java 8