使用BAM格式的成对读取

BamTyper的Python项目详细描述


88888888888
888"88b"888
888.888p 888
88888 k.888b.888b.d888b.8888888888B.D888B.888D888
888"Y88B"888"888"88B 888 888"888 888"88B D8P"Y8B 888P"
888{888.D88888888 888 888 888)888 888 888/888 888 888/888 888 888 888 888 88888)88888888888888 88888 Y8B.888
888 P""Y8888 888"888"Y888"888"Y888 888P""Y8888"888
888 888
Y8B D888
"Y888"888;概述




脚手架)使用管道:




/usr/bin/env python
来自bamtyper.utilities导入bamparser btbp
bp bp=btbp()
bam some types=bp types(bamfiles)
其中:
bamfiles>bam名称列表br brERT大小和BAM文件中读取的相对方向}
其中:
类型--方向类型:
0:输出:<;---->;
1:相同:-->;-->;
2:输入-(>;--stdev-插入片段大小的标准偏差

>35;获取链接对:
35;!/usr/bin/env python
from bamtyper.utilities import bamparbe btbp
bp bp=btbp()
(links,ref of lengthlengthlength,coverage)=bp.getlink(bamfiles,docoverage=true)

ct包含两个contig之间的链接信息
{c1:(c2,链接链接,类型gap),…和}
其中:
在链接上-悬停读取的次数确认链接
链接类型-两个控件(开始和结束)的相对方向
ss、se、es、ee或erroren te two contig s

coverage s——包含BAM文件中每个contig的片段覆盖率的dict
{c1[(cov1,cov2,-(covn),……}

ref lenges——包含所有contigs的长度的dict
{c1:len,…ref lenges}

注意:
typebamr将自动计算每个文件中读取的
的方向和插入大小,并在此基础上估计链接的类型和间隔

35;在命令上:
typer type>分析文件并确定文件键入

$bam typer type bam file.bam


确定bam"bamfile"的OT
方向:在insert:301中,stdev:29
usabr links:
br links>.bam

1.
contig2,[contig1,39,se,69]
contig1,[contig2,39,es,69]


],[contig1,4,ee,45]
contig1,[contig2,4,ee,45]




$bamtyper链接bamfile1.bam bamfile2.bam-c

contig2,[contig3,3,ss,58],[contig1,4,ee,45]
contig1,[contig2,4,ee,45]
contig3 0.6206 0.5234
contig2 0.6558 0.0123
contig1 0.6523 0.5634

f 62x和52x分别

管理

有关详细信息,请参见license.txt。

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

推荐PyPI第三方库


热门话题
在OSGI中使用cxf生成的客户端时出现Java类装入器问题和JaxB异常   java为什么要在javamail中迭代多部分电子邮件中的部分?   并发编程问题   JFileChooser&&System中未调用java windowClosing。退出功能不正常?   SQL查询的java语法分析   java如何使用AspectJ声明字段上的警告   什么是java向量。元素()C#等价物   java解析Android应用程序中tornado web服务中的CSV文件   java我试过c2dm,我需要服务器端   java调整JPanel大小以适应新的JLabel图标   Java与Python脚本的通信   java使用Saxon通过XSLT生成URL   java net::ERR_complete_CHUNKED_编码200(OK)来自struts应用程序中的tomcat   java如何为我的窗格设置不同的位置?   java使用Cypher Neo4j获取给定类型的所有节点(从SQL世界中的tablename中选择*)   nio使用Java解析文件值   java使用WSDL生成REST客户端会是错误的方向吗?   java如何在我的应用程序中构建类映射?   java按钮。setEnabled在第一个循环中不起作用   xPath适用于最后一页,但不适用于第一页