用snakemake更改目录名和fastq文件名

问题描述

我想更改文件夹名称和读取文件的名称。从这里我发现了类似的东西 (Move and rename files from multiple folders using snakemake):

workdir: "/path/to/workdir/"

import pandas as pd
from io import StringIO

sample_file = StringIO("""fastq sampleID
BOB_1234/fastq/BOB_1234.R1.fastq.gz TAG_1/fastq/TAG_1.R1.fastq.gz
BOB_1234/fastq/BOB_1234.R2.fastq.gz TAG_1/fastq/TAG_1.R2.fastq.gz
BOB_3421/fastq/BOB_3421.R1.fastq.gz TAG_2/fastq/TAG_2.R1.fastq.gz
BOB_3421/fastq/BOB_3421.R2.fastq.gz TAG_2/fastq/TAG_2.R2.fastq.gz""")

df = pd.read_table(sample_file,sep="\s+",header=0)
sampleID = df.sampleID
fastq = df.fastq

rule all:
    input:
       expand("{sample}",sample=df.sampleID)

rule move_and_rename_fastqs:
    input: fastq = lambda w: df[df.sampleID == w.sample].fastq.tolist()
    output: "{sample}"
    shell:
        """echo mv {input.fastq} {output}"""

我可以出错:

MissingOutputException in line 19 of snakefile:
Job Missing files after 5 seconds:
TAG_1/fastq/TAG_1.R1.fastq.gz

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)