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

oldowan.mitomotifs的Python项目详细描述


mitotimotis是一个小型、纯python的生物信息实用程序,用于(1) 将人类有丝分裂软骨DNA序列数据转换为与 修正剑桥参考序列(rcrs)和(2)变换变异位点 数据转换成DNA序列。关于RCR和变型站点的进一步信息 人类mtdna序列的命名可在MitoMotifs网站上找到。

安装说明

这个包是纯python的,没有标准之外的依赖项 图书馆。最简单的安装方法是使用 setuptools包。通常是这样的:

$ easy_install oldowan.mitomotifs

或者在类unix系统上,假设您正在安装到主python site-packages目录作为非特权用户,此:

$ sudo easy_install oldowan.mitomotifs

您还可以使用标准的python distutils设置方法。下载 从文件列表到本页底部的当前源存档, 把它拆开,然后安装。在MacOSX和许多其他类似Unix的系统上, 下载了存档并更改为包含此存档的目录 你的外壳,可能会像这样:

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

快速启动

从oldowan.mitotimotis导入seq2sitessites2seq

>>> from oldowan.mitomotifs import seq2sites, sites2seq

将序列转换为站点:

>>> seq = """TTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAA
GTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTACTGCC
AGCCACCATGAATATTGTACAGTACCATAAATACTTGACCACCTGTAGT
ACATAAAAACCCAATCCACATCAAAACCCCCTCCCCATGCTTACAAGCA
AGTACAGCAATCAACCTTCAACTATCACACATCAACTGCAACTCCAAAG
CCACCCCTCACCCACTAGGATACCAACAAACCTACCCACCCTTAACAGT
ACATAGTACATAAAGCCATTTACCGTACATAGCACATTACAGTCAAATC
CCTTCTCGTCCC"""
>>> seq2sites(seq)

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

也有一个边界的数量的矛盾地点(n)允许在 顺序。默认情况下,这是10-但这是一个可以设置的选项:

>>> seq2sites(seq, ambig_cutoff=20)

将变量站点列表转换为序列。默认的序列区域 返回的是超变量区域1(hvr1),即位置16024到16365 (以生物为基础的编号):

>>> sites2seq('16129A 16223T')

预定义的序列区域是:

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

因此,要将hvr2站点列表转换为序列:

>>> sites2seq('73G', region='HVR2')

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

>>> sites2seq(['16129A', '16223T', '73G'], region='HVR1to2')

在给定空字符串、空列表或 字符串“rcrs”。以下各项都是等效的:

>>> sites2seq('')
>>> sites2seq([])
>>> sites2seq('rCRS')

通过将站点列表传递给 region选项:

>>> sites2seq('', region=[1,2,3])

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

>>> sites2seq('', region=range(73,341))  # include 340, but not 341

发布历史记录

1.0.0(2008年8月16日)
模块的初始版本。
1.0.1(2008年8月22日)
为简化的hvr1站点(即16129a为129a)在站点2seq中添加新的“add16k”选项

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

推荐PyPI第三方库


热门话题
java Listview在setcontentview()中变得模棱两可   将Java集合转换为Clojure数据结构   java应用程序发布(process com.example.ship)意外停止。请再试一次   java折叠工具栏不可滚动   java struts2<s:选择   java子集数   单例“getInstance”函数的未解决Java编译问题   java如何在SWT中显示React UI   java如何将文件写入PC上可见的Android SD卡   从导航抽屉中打开新片段时,java将content_main替换为   java如何作为web服务安全地执行shell命令   java Web视图不显示HTML   java如何通过soap发送单向请求?   未创建java共享首选项