全基因组分析相关基因组学Omics Academy

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

2020-10-26  本文已影响0人  OmicsAcademy

GenomeView是什么?

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

如何安装?

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

上一篇 下一篇

猜你喜欢

热点阅读