Added reference related logic.

config/default.yml

@@ -46,7 +46,7 @@ reference:
   file:
     - pc_transcripts
     - lncRNA_transcripts
-  spikeIns: 'ref/humanDux.fasta'
+  #spikeIns: 'path/to/spikeIns.fa'
 
 #
 # The salmon section

sc-preprocess.snake

@@ -30,7 +30,7 @@ SAMPLE_NAMES = sorted(SAMPLES.keys())
 #-------------------------------------------------------------------------------#
 
 output_files = [
-	#join(config['dirs']['ref'], 'tx2gene', basename(REFERENCE_FASTA).rstrip(".fa"))
+	join(config['dirs']['ref'], 'tx2gene', basename(REFERENCE_FASTA).rstrip(".fa"))
 ]
 
 if config["action"]["demultiplex"]:
@@ -47,7 +47,7 @@ if config["action"]["clean"]:
 
 include: "src/auxiliary.snake"
 #include: "src/demultiplex.snake" # not yet implemented
-#include: "src/reference.snake"
+include: "src/reference.snake"
 include: "src/clean.snake"
 
 if config["debug"]:

src/reference.snake

@@ -0,0 +1,52 @@
+# vim: syntax=python tabstop=4 expandtab
+# coding: utf-8
+
+'''
+@author: jpreuss
+
+Provides rules for reference related modifications
+'''
+
+from snakemake.remote.FTP import RemoteProvider as FTPRemoteProvider
+FTP = FTPRemoteProvider()
+
+rule reference_download:
+  input:
+    FTP.remote(expand(config['reference']['URL'], **config['reference']), keep_local = True)
+  output:
+    GENCODE_FASTA
+  threads: 1
+  message:
+    'Downloading gencode transcriptome reference.'
+  shell:
+    """
+    zcat {input} > {output}
+    """
+
+if 'spikeIns' in config.get('reference', {}):
+  rule reference_addSpikes:
+    input:
+      ref = GENCODE_FASTA,
+      spike = config['reference']['spikeIns']
+    output:
+      join(config['dirs']['ref'], config['reference']['organism'], config['reference']['release'] + '-withSpikeIns.fa')
+    threads: 1
+    message:
+      'Adding spike-in sequences to gencode transcriptome reference'
+    shell:
+      """
+      cat {input.ref} {input.spike} > {output}
+      """
+
+rule tx2gene_from_fasta:
+  input:
+    REFERENCE_FASTA
+  output:
+    join(config['dirs']['ref'], 'tx2gene', basename(REFERENCE_FASTA).rstrip(".fa"))
+  threads: 1
+  message:
+    'Creating tx2gene table from reference {input}.'
+  shell:
+    """
+    grep ">" {input} | tr -d '>' | cut -d'|' --output-delimiter=$'\t' -f1,2 | sed -e 's/\([A-Z]\+[[:digit:]]\+\)\(\.[[:digit:]]\+\)$/\\1/g' > {output}
+    """