pysam操作sam文件
2024-09-01 04:59:20
pysam模块
因为要分析sam文件中序列的情况,因此要对reads进行细分,所以之前想用数据库将sam文件信息存储,然后用sql语句进行分类。后来发现很麻烦,pysam就是一个高效读取存储在SAM / BAM / CRAM格式文件中的映射短读序列数据信息的python模块,可以轻松地对reads进行操作。
1.安装Pysam
$ pip install pysam
2.检查是否安装成功
import pysam # 注意,此步是进入python交互环境
3.读取bam文件
import pysam
samfile = pysam.AlignmentFile("ex1.bam", "rb")
提取指定区域reads
for read in samfile.fetch('chr1', 100, 120):
print read samfile.close()
结果如下:
EAS56_57:6:190:289:82 0 99 <<<7<<<;<<<<<<<<8;;<7;4<;<;;;;;94<; 69 CTCAAGGTTGTTGCAAGGGGGTCTATGTGAACAAA 0 192 1
EAS56_57:6:190:289:82 0 99 <<<<<<;<<<<<<<<<<;<<;<<<<;8<6;9;;2; 137 AGGGGTGCAGAGCCGAGTCACGGGGTTGCCAGCAC 73 64 1
EAS51_64:3:190:727:308 0 102 <<<<<<<<<<<<<<<<<<<<<<<<<<<::<<<844 99 GGTGCAGAGCCGAGTCACGGGGTTGCCAGCACAGG 99 18 1
...
参考资料
最新文章
- c#编程基础之ref、out参数
- 修改sys密码与nbu备份脚本密码后,nbu备份报密码无效
- Mysql监控、优化
- css3学习总结4--CSS3背景
- wordpress无法安装这个包。: PCLZIP_ERR_MISSING_FILE (-4) : Missing archive file &#39;C:\WINDOWS\TEMP/wordpress-4.tmp&#39;
- Android 常用权限
- DB2 insert into 三种写法
- SP_CreateInsertScript 将表内的数据全部拼接成INSERT字符串输出
- 05-UIKit(UITableViewController)
- [ACM] POJ 3687 Labeling Balls (拓扑排序,反向生成端)
- LeetCode &; Q88-Merge Sorted Array-Easy
- HihoCoder1338 A Game (区间DP)
- Java编程基础篇第四章
- edit this cookie chrome插件 (HttpAnalyzerStdV3 获取Cookie 后,再用edit this cookie添加cookie)
- 【Tomcat】性能优化
- android--------性能优化之Allocation Tracker
- android--Activity有返回值的跳转
- 胡乱摸的NOIP2017游记和总结
- 使用stream(流)实现多表数据传输
- [CodePlus 2017 11月赛]晨跑 题解(辗转相除法求GCD)