从snakemake下载fastq

Question

在执行以下规则之前，我正在拼命尝试让一个规则下载我的fastq文件。我尝试了很多方法，包括这里建议的：http://ivory.idyll.org/blog/tag/snakemake.html

下面是我的snakemake的简化版本：

######### functions

def read_samplesTable(inputTable):
    data = pandas.read_csv(inputTable)
    # Verify column names
    if not {'run', 'organism', 'name', 'experiment_title', 'cell_line', 'rep', 'study_name', 'library_strategy', 'library_layout', 'study_title'}.issubset(data.columns.values):
            raise KeyError("The samples file must contain the following named columns: 'run', 'organism', 'name', 'experiment_title', 'cell_line', 'rep', 'study_name', 'library_strategy', 'library_layout', 'study_title'")
    return data

def retrieveName(description):
    result = []
    for items in description.iteritems():
        result.append(items[1].split(":")[0])
    return result 

######### Variables

input_table = config["samples"]["summaryFile"]
samplesData = read_samplesTable(input_table)

index_single = samplesData['library_layout'] == 'SINGLE - '
samplesData_single = samplesData[index_single]

gsm_single = retrieveName(samplesData_single["experiment_title"])
outputName_single = samplesData_single['name'] + "_" + samplesData_single['run'] + "_" + gsm_single + "_" + samplesData_single['study_name'] + "_" + samplesData_single['cell_line'] + "_" + samplesData_single['rep']
single_samples = outputName_single.tolist()

names_srrID_single = samplesData_single['run']

############ Rule

rule all:
  input:
    expand("data/single/{singleEndName}.fastq.gz", singleEndName = single_samples)

rule download_fastq_single:
  output:
    singleFastq = "data/single/{singleEndName}.fastq.gz"
  params:
    outputdirectory = config["rawdata"]["fastqrootfolder"]
    ssridsingle = lambda wildcards: samplesData_single.loc[wildcards.names_srrID_single, "run"]    
  shell:
    "fastq-dump --accession {params.srridsingle} --defline-seq '@$sn[_$rn]/$ri' --defline-qual \'+\' --gzip --outdir {params.outputdirectory}"

规则在一个单独的文件'download-fastq-snakefile‘中，我得到了错误SyntaxError in line 6 of download-fastq-snakefile。

正如我所料，问题来自ssridsingle = lambda wildcards: samplesData_single.loc[wildcards.names_srrID_single, "run"]行

如果你能帮我，那就太好了！

谢谢

Answer 1

我猜问题只是缺少了一个逗号：

  params:
    outputdirectory = config["rawdata"]["fastqrootfolder"], <-- add this comma
    ssridsingle = lambda wildcards: samplesData_single.loc[wildcards.names_srrID_single, "run"]    

Answer 2

是的，这是第一个问题，谢谢！

对于第二个问题，我必须修改对srr id的访问，如下所示：

rule download_fastq_single:
  output:
    singleFastq = "data/single/{singleEndName}.fastq.gz"
  params:
    outputdirectory = config["rawdata"]["fastqrootfolder"],
    ssridsingle = lambda wildcards: names_srrID_single    
  shell:
    "fastq-dump --accession {params.ssridsingle} --defline-seq '@$sn[_$rn]/$ri' --defline-qual \'+\' --gzip --outdir {params.outputdirectory}"