将人类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
最简单的安装方法是使用
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