Added flexibility for gz/plain reference fasta files

jenzopr · Apr 26, 2018 · dea6443 · dea6443
1 parent 820f9a8
commit dea6443
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/reference.snake b/src/reference.snake
@@ -20,7 +20,7 @@ rule reference_download:
     'Downloading gencode transcriptome reference.'
   shell:
     """
-    zcat {input} > {output}
+    zcat -f {input} > {output}
     """
 
 if 'spikeIns' in config.get('reference', {}):
@@ -35,7 +35,7 @@ if 'spikeIns' in config.get('reference', {}):
       'Adding spike-in sequences to gencode transcriptome reference'
     shell:
       """
-      cat {input.ref} {input.spike} > {output}
+      zcat -f {input.ref} {input.spike} > {output}
       """
 
 rule tx2gene_from_fasta: