简易序列比对图

shabam的Python项目详细描述


沙巴姆

绘制BAM或CRAM序列读取的python工具。在

安装

安装cairo(如果尚未安装) 已安装:

# macOS via conda, or via homebrew (choose one, then set the library path)
conda install cairo pkgconfig
brew install cairo pkg-config

安装沙巴姆:

^{pr2}$

梦想

fromshabamimportseqplotseqplot('example.bam',chrom='1',start=30243,end=30321,fastafile='reference.fasta',out='plot.svg')

Reality

绘图选项

  • 使用by_strand=True对字符串进行着色读取
  • 将多个序列文件与路径列表一起绘制,例如。 ['child.bam', 'mom.bam', 'dad.bam']
  • 导出具有匹配文件扩展名的PDF、PNG、SVG或PS格式的绘图

命令行版本

shabam \
  --seqfiles tests/data/example.bam \
  --chrom 1\
  --start 30243\
  --end 30321\
  --fastafile tests/data/reference.fasta \
  --out plot.svg

梦想的进一步改善

  • []当我们没有提供参考序列时,使用一致性序列
  • []允许自定义颜色
  • []计算任何站点的变体比例
  • []绘图读取深度
  • []在变量比例大于阈值的位置,反映比例 带底色
  • []用于缩放打印基大小的选项,当前为每基10像素
  • []可选阴影图深变化样式:
    • 红通道:核苷酸
    • 蓝色通道:读串
    • 绿色通道:基础质量
    • alpha:base支持ref或alt
  • []在高深度序列数据中展平垂直绘图
  • []针对极高深度序列数据(大于1000倍)的向下采样读取

信贷

初始雪茄解析代码已从的权限中提取 pybamview。在

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

推荐PyPI第三方库


热门话题
java能否使用OpenJPA将一个实体持久化到多个持久化单元?   java如何访问当前按下的键   java singleton是否意味着哈希代码总是返回相同的值?   inputstream无法在Java中同时使用STDIN和STDERR读取程序输出   java运行时有多贵。getRuntime()与获取内存度量相结合   java AndroidStudio ListView适配器更新   java小程序请求未扩展实时会话   java是设置布尔值所必需的关键部分?   java比较两个字符串,并突出显示发现的不匹配项   java带抽屉布局,无论哪个屏幕处于活动状态,如何在按下后退键时关闭应用程序?   爪哇:颜色有什么区别。黑色和彩色。黑色   Velocity模板中子类的java访问方法   java如何快速学习Drools或其他规则引擎   从Java应用程序访问时缓存数据库查询结果   java cassandra nodetool JPLISAgent。c错误   java我正在解析一个没有pubDate的RSS提要,有没有其他方法可以确定一个项目是何时发布的?