snakemake学习笔记003:stringtie合并转录本
2022-05-05 本文已影响0人
小明的数据分析笔记本
SRR, = glob_wildcards("output.gtf/"+"{srr}.gtf")
#SRR = ["ERR188401","ERR204916"]
#rule all:
# input:
# expand("output.gtf/"+"{srr}.gtf",srr=SRR),
rule gtflist:
input:
gtffiles = expand("output.gtf/"+"{srr}.gtf",srr=SRR)
output:
output_txt = "MergedList.txt"
run:
with open(output.output_txt,'w') as f:
for gtf in input.gtffiles:
print(gtf,file=f)
第一行SRR后面有一个逗号
最开始没写这个逗号,一直遇到报错
image.pngBuilding DAG of jobs... MissingInputException in line 3 of /mnt/shared/scratch/myan/private/practice_data/RNAseq/snakemake.rnaseq/gtf_list.py: Missing input files for rule all: output.gtf/['ERR188337', 'ERR188245', 'ERR188428', 'ERR188401', 'ERR204916', 'ERR188383', 'ERR188104', 'ERR188257', 'ERR188273', 'ERR188044', 'ERR188234', 'ERR188454'].gtf