将人类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第三方库


热门话题
java RxJava无重复合并   java在JADE中的同一容器中创建多个代理   带有嵌入式commandLink的java Primefaces数据表   当我必须使用文件读取器Java时,如何从txt文件填充数组??   JavaSpring从数据库和列表中删除记录   java将工作簿对象转换为文件对象   java组织。阿帕奇。火花sql。AnalysisException:无法解析''S.SID`'   java如何将int值从一个活动传递到另一个活动并在TextView中显示?   java如何将pom项目打包为可运行jar   java如何定义标准。不同的根实体实际工作吗?   如何在java中读取文件并将其存储在数组中   通过位操作将Java代码从Java转换为JavaScript代码   粘贴后,java JTextField密钥侦听器落后   java在socket和服务器socket之间发送和接收信息   在Java中处理大型输入流以列出字符串   java如何创建一个示例FastLoad1。csv文件看起来像什么?   可视化和实现椭圆曲线密码的java工具   MS Project Server 2010 Java API   java Jackson错误地序列化了XMLGregorianCalendar