斯奈克马克生物信息学图书馆。

SMBL的Python项目详细描述


https://travis-ci.org/karel-brinda/smbl.svg?branch=masterhttps://travis-ci.org/karel-brinda/smbl.svg?branch=devel

如果有任何问题,请随时与我联系。

简短说明

smbl是一些有用的规则和python函数库,可用于snakemake(https://bitbucket.org/johanneskoester/snakemake/)管道。它使得自动 安装各种生物信息学程序,如读取映射程序、读取模拟器、转换工具等。 它还支持下载和转换fasta格式的一些重要参考文献(如人类基因组)。

安装/升级

要安装smbl,您至少需要有类unix的操作系统(例如linux、macos)和python 3.3。 可以使用以下命令执行安装/升级。

pip3 install --upgrade smbl

如果还没有安装snakemake,它将 与SMBL一起自动安装。

git中smbl的当前版本可以通过以下方式安装

pip3 install --upgrade git+git://github.com/karel-brinda/smbl

要求

为了能够自动下载和安装软件,smbl要求unix系统中存在以下程序:

  • wget或curl
  • 合同一般条件4.7+
  • 吉特
  • 制造

用法

要使用smbl,您必须导入smblpython包,并包含一个包含所有规则的文件,使用:

importsmblinclude:smbl.include()

然后您可以使用所有支持的程序或数据。当它们作为规则的输入出现时,它们将被下载或编译。

所有程序都安装到~/.smbl/bin/,所有fasta文件都安装到~/.smbl/fa/

程序

ProgramVariable with its filenameLink
art_454^{tt3}$http://www.niehs.nih.gov/research/resources/software/biostatistics/art/
art_illumina^{tt4}$http://www.niehs.nih.gov/research/resources/software/biostatistics/art/
art_solid^{tt5}$http://www.niehs.nih.gov/research/resources/software/biostatistics/art/
bcftools^{tt6}$http://github.com/samtools/bcftools
bfast^{tt7}$http://github.com/nh13/bfast
bgzip^{tt8}$http://github.com/samtools/htslib
bowtie2^{tt9}$http://github.com/BenLangmead/bowtie2
bowtie2-build^{tt10}$http://github.com/BenLangmead/bowtie2
bowtie2-inspect^{tt11}$http://github.com/BenLangmead/bowtie2
bwa^{tt12}$http://github.com/lh3/bwa
curesim.jar^{tt13}$http://www.pegase-biosciences.com/tools/curesim/
curesim_eval.jar^{tt14}$http://www.pegase-biosciences.com/tools/curesim/
deez^{tt15}$http://github.com/sfu-compbio/deez
drfast^{tt16}$http://github.com/BilkentCompGen/drfast
dwgsim^{tt17}$http://github.com/nh13/dwgsim
dwgsim_eval.pl^{tt18}$http://github.com/nh13/dwgsim
freec^{tt19}$http://bioinfo-out.curie.fr/projects/freec/
gem-indexer^{tt20}$http://algorithms.cnag.cat/wiki/The_GEM_library
gem-mapper^{tt21}$http://algorithms.cnag.cat/wiki/The_GEM_library
gem-2-sam^{tt22}$http://algorithms.cnag.cat/wiki/The_GEM_library
gnuplot4^{tt23}$http://www.gnuplot.info/
gnuplot5^{tt24}$http://www.gnuplot.info/
kallisto^{tt25}$https://github.com/pachterlab/kallisto
lastal^{tt26}$http://last.cbrc.jp/
lastdb^{tt27}$http://last.cbrc.jp/
mason_frag_sequencing^{tt28}$http://packages.seqan.de/mason2/
mason_genome^{tt29}$http://packages.seqan.de/mason2/
mason_materializer^{tt30}$http://packages.seqan.de/mason2/
mason_methylation^{tt31}$http://packages.seqan.de/mason2/
mason_simulator^{tt32}$http://packages.seqan.de/mason2/
mason_splicing^{tt33}$http://packages.seqan.de/mason2/
mason_variator^{tt34}$http://packages.seqan.de/mason2/
mrfast^{tt35}$http://github.com/BilkentCompGen/mrfast
mrsfast^{tt36}$http://mrsfast.sourceforge.net/
perm^{tt37}$http://code.google.com/p/perm/
pbsim^{tt38}$https://code.google.com/p/pbsim
picard^{tt39}$http://broadinstitute.github.io/picard/
sambamba^{tt40}$http://lomereiter.github.io/sambamba/
samtools^{tt41}$http://github.com/samtools/samtools
sirfast^{tt42}$http://github.com/BilkentCompGen/sirfast
storm-color^{tt43}$http://bioinfo.lifl.fr/yass/iedera_solid/storm/
storm-nucleotide^{tt44}$http://bioinfo.lifl.fr/yass/iedera_solid/storm/
tabix^{tt45}$http://github.com/samtools/htslib
twoBitToFa^{tt46}$http://hgdownload.cse.ucsc.edu/admin/exe/
vcfutils.pl^{tt47}$http://github.com/samtools/bcftools
wgsim^{tt48}$http://github.com/lh3/wgsim
wgsim_eval.pl^{tt49}$http://github.com/lh3/wgsim
xs^{tt50}$http://bioinformatics.ua.pt/software/xs/

fasta文件

FASTA fileVariable with its filename
An example small FASTA file^{tt51}$
An example small FASTA file^{tt52}$
An example small FASTA file^{tt53}$
Human genome HG38 (GRCh38)^{tt54}$, ^{tt55}$
Mouse genome MM10^{tt56}$
Chimpanzee genome PANTR04^{tt57}$

示例

下面的示例演示如何使用smbl自动安装软件。

创建名为Snakefile的空文件,其内容如下:

importsmblinclude:smbl.include()ruleall:input:smbl.prog.DWGSIM,smbl.prog.BWA,smbl.fasta.EXAMPLEparams:PREF="simulated_reads",INDEX="bwa_index"output:"alignment.sam"run:# read simulationshell("{input[0]} -C 1 {input[2]} {params.PREF}")# creating BWA index of the reference sequenceshell("{input[1]} index {input[2]}")# mapping by BWAshell("{input[1]} mem {input[2]} {params.PREF}.bfast.fastq > alignment.sam")

运行脚本

snakemake

发生了什么:

  1. 下载了一个示例fasta文件
  2. dwgsim和bwa被下载、编译和安装
  3. dwgsim模拟对示例fasta文件的读取
  4. 这些读取被bwa映射回引用(alignment.sam被创建)

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

推荐PyPI第三方库


热门话题
java从JSP或HTML向servlet发送多个参数   java方法来查看字符是否在字符数组中   使用带有java的MAC地址连接到设备   java如何将csv文件中的数据打印到secondactivity?   java如何从netbean 7.0.1连接到数据库   java考虑所有可能的类值,用于输出测试分割的预测值。   java我的actionListener调用有什么问题   swing在Java中实现粒子过滤器最有效的方法是什么?   java运行。getFontFamily()为返回null。使用apachepoi的docx文件   一个事务中的java领域循环与每个步骤循环中的一个事务   java日期格式与Spring Boot不兼容   java类冲突。处理   java GridBagLayout不工作   java将图像发送到另一个应用程序