过程介绍
1, 安装snakemake
2, 新建文件
3, 新建一个简单的Snakemake参数文件
4, 扩展, 去关联输出文件
5, 使用全局变量, 关联文件
6, 批量运行
1, 安装snakemake...这里需要时python3, 不支持python2
pip3 install --user snakemake pyaml
2, 新建几个FASTQ文件
这里, 我们新建两个配对的RNA-seq数据,...格式是FASTQ的文件, 然后经过下面两步处理:
第一步: 数据质量控制
第二部: 将基因表达合并为一个文件
创建文件
创建genome.fa文件, 使用touch创建空文件即可
创建fastq文件夹..., SAMPLES,里面有两个元素: Sample1和Sample2
定义一个rule, 名称为all, input使用expand函数, 能够将数组的内容解析给{sample}
rule all:..., 使用命令:
snakemake -np
参数介绍
-n 或者—dryrun, 表示只生成命令, 但是不执行命令, 可以预览一下生成的命令.