snakemake流程開發(fā)

官方文檔https://snakemake.readthedocs.io/en/stable/
使用場(chǎng)景:搭建生信分析流程,可重復(fù)性使用,支持分析中斷繼續(xù)分析的功能。
安裝:推薦使用conda安裝。
先安裝conda,然后使用conda安裝mamba軟件

$ conda install -n base -c conda-forge mamba

激活base環(huán)境,使用mamba安裝snakemake

$ conda activate base
$ mamba create -c conda-forge -c bioconda -n snakemake snakemake

激活snakemake環(huán)境,測(cè)試是否安裝成功

$ conda activate snakemake
$ snakemake —help

怎么用?

需要先了解snakemake的語法結(jié)構(gòu)

rule map_reads:
    input:
        "data/genome.fa",
        "data/samples/A.fastq"
    output:
        "results/mapped/A.bam"
    shell:
        "bwa mem {input} | samtools view -b - > {output}"

主要需要三部分:input,output和run(執(zhí)行shell命令時(shí)用shell),在run部分可以編寫python代碼

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容