阅读 78

GenomeView: 基于Python的基因组数据展示工具

GenomeView是什么?

GenomeView通过Python直接展示基因组数据。主要的特点包括:

  • 易扩展

  • 可与Jupyter notebook整合

  • 高质量图片生成

  • 兼容不用类型的数据:BAM文件(SE和PE read)和图形数据文件例如基因组Coverage,wiggle文件等等

    GenomeView的目的不是像IGV等基因组浏览器一样,产生可交互的数据展示。当然通过Jupyter notebook,用户也可以很简单的快速生成新的展示图。
    

如何安装?

GenomeView: Python 3.3 or greater.

下面的命令行可以帮助你安装GenomeView

pip install genomeview  ## --user 如果没有root权限</pre>

或者直接从GitHub的源代码安装:

pip install -U git+https://github.com/nspies/genomeview.git  ## --user 如果没有root权限

如果想展示bigWig 文件,pyBigWig 包也需要安装

pip install pyBigWig

快速入门

image
    通过下面的代码可以产生上面的图:
import genomeview

dataset_paths = ["data/pacbio.chr1.bam",
                 "data/illumina.chr1.bam",
                 "/Users/nspies/Downloads/hg19.refseq.sorted.bed.gz"]
reference = "ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/technical/reference/phase2_reference_assembly_sequence/hs37d5.fa.gz"

chrom = "chr1"
start = 224368899
end =   224398899

doc = genomeview.visualize_data(dataset_paths, chrom, start, end, reference)

如果你使用Jupyter notebook, 直接在最后一行输入doc

doc

你也可以使用 genomeview.save() 来将图片存入一个文件中。

genomeview.save(doc, "/path/to/output.svg")   # or .png/.pdf

更多细节请参考:http://genomeview.readthedocs.io/en/latest/index.html

Reference:

http://genomeview.readthedocs.io/en/latest/index.html

https://github.com/nspies/genomeview

作者:OmicsAcademy

原文链接:https://www.jianshu.com/p/ec14411b254f

文章分类
后端
版权声明:本站是系统测试站点,无实际运营。本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 XXXXXXo@163.com 举报,一经查实,本站将立刻删除。
相关推荐