是用pip下载的pysam，为了做index文件用brew install下载了samtools。
import pysam
读取bam文件如下。
n1 = pysam.AlignmentFile('./atac_bam/scN1.pe.q10.sort.rmdup.bam', "rb")
为了能够读取特定位置的read，用samtools做index文件。
samtools index scN1.pe.q10.sort.rmdup.bam
然后用fetch的方法读取行。fetch返回的是一个迭代器(iterator),可以迭代读取内容。
fetch的API如下：
fetch(self, reference=None, start=None, end=None, region=None, tid=None, until_eof=False, multiple_iterators=False)
for read in n1.fetch('chr6'):
print read
n1.close()
是用pip下载的pysam，为了做index文件用brew install下载了samtools。 读取bam文件...
Pysam操作BAM文件 Pysam包是一个处理基因组数据的python模块，它打包了htslib-1.3、sam...
Pysam 是一个python的模块，可用于处理bam文件 安装 使用 Pysam的函数有很多，主要的读取函数有：...
pysam包使用AlignmentFile对象存储BAM/SAM文件，而通过AlignmentSegment对象存...
pysam是python下的处理bam、sam、bcf等序列比对文件格式的功能强大的包！linux下的安装过程： ...
Pysam作为一个python模块可以读取与处理SAM/BAM文件，并提供了方便的对象与丰富的操作接口，让我们了解...
在python中读取、处理文件可以用pysam这个包。以下简单介绍一下这个包的使用。 读取文件 仅读取某条染色体某...
使用pysam包读取BAM/SAM文件返回的AlignmentFile对象包含丰富的属性。 header属性 he...
如何使用deeptools处理BAM数据 总体介绍 deeptools是基于Python开发的一套工具，用于处理诸...
现在我们已经具备采用pysam模块统计比对信息的知识了，这里我们编写一个程序来统计BAM/SAM文件的信息并使用s...
本文标题：用pysam处理bam数据
本文链接：https://www.haomeiwen.com/subject/qumdvqtx.html
网友评论