resolved typo in input rule for single-end reads

workflow/rules/common.smk

@@ -377,7 +377,14 @@ def get_star_input(wildcards):
   elif config['data']['rnaseq_filetype'] == '.fq' or config['data']['rnaseq_filetype'] == '.fastq':
     if config['preproc']['activate']:
       if config['data']['rnaseq_readtype'] == 'SE':
-        return expand("results/{sample}/rnaseq/preproc/reads.fq.gz", **wildcards)
+        return dict(
+            zip(
+              ["fq1"],
+              expand("results/{sample}/rnaseq/reads/{group}_preproc.fq.gz",
+                      sample=wildcards.sample,
+                      group=wildcards.group)
+            )
+        )
       elif config['data']['rnaseq_readtype'] == 'PE':  # PE
         return dict(
             zip(
@@ -390,7 +397,7 @@ def get_star_input(wildcards):
         )
       else: # no pre-processing
         return config['data']['rnaseq'][wildcards.group]
-
+    
 # collect the individual alignments from splitted bamfiles
 def aggregate_aligned_rg(wildcards):
     # make sure that all samples are processed in checkpoint - split fastq file

workflow/rules/hlatyping.smk

@@ -30,6 +30,7 @@ rule index_hla_panel:
     "../envs/yara.yml"
   shell:
     """
+      mkdir -p resources/hla/yara_index/
       yara_indexer \
           -o resources/hla/yara_index/{wildcards.nartype} \
           {input.fa} > {log}

workflow/rules/quantification.smk

@@ -48,8 +48,7 @@ rule countfeatures_rnaseq:
     sample = "results/{sample}/rnaseq/align/{group}_final_STAR.bam",
     annotation_file = "resources/refs/genome.gtf"
   output:
-    table = "results/{sample}/rnaseq/quantification/{group}_counts.txt",
-    summary = "results/{sample}/rnaseq/quantification/{group}_counts.txt.summary",
+    "results/{sample}/rnaseq/quantification/{group}_counts.txt"
   log:
     "logs/{sample}/featurecounts/rnaseq_{group}.log",
   threads: 2