comparison Snakefile @ 12:849ab444b4f1 draft

Uploaded
author chrisw
date Wed, 13 Feb 2019 15:51:45 -0500
parents 35ed7314038d
children af1f49e5242d
comparison
equal deleted inserted replaced
11:35ed7314038d 12:849ab444b4f1
11 'bamcount_auc.tsv', 11 'bamcount_auc.tsv',
12 'bamcount_frag.tsv', 12 'bamcount_frag.tsv',
13 'Chimeric.out.junction.zst', 13 'Chimeric.out.junction.zst',
14 'all.exon_bw_count.zst', 'unique.exon_bw_count.zst', 14 'all.exon_bw_count.zst', 'unique.exon_bw_count.zst',
15 'manifest'] 15 'manifest']
16
17 import os
18 config['output']=os.path.abspath(config['output'])
16 19
17 import subprocess 20 import subprocess
18 def run_command(cmd_args): 21 def run_command(cmd_args):
19 cmd_args = ' '.join(cmd_args) 22 cmd_args = ' '.join(cmd_args)
20 try: 23 try:
24 raise cpe 27 raise cpe
25 28
26 29
27 import re 30 import re
28 FASTQ_PATT=re.compile(r'([^_\.]+)(_(\d+))?\.((fastq)|fq)(\.gz)?$') 31 FASTQ_PATT=re.compile(r'([^_\.]+)(_(\d+))?\.((fastq)|fq)(\.gz)?$')
29 import os
30 def prep_for_galaxy_run(): 32 def prep_for_galaxy_run():
31 try: 33 try:
32 os.mkdir(config['temp']) 34 os.mkdir(config['temp'])
33 except OSError as ose: 35 except OSError as ose:
34 pass 36 pass
35 fastqs = config['inputs'].split(',') 37 fastqs = config['inputs'].split(',')
36 m = FASTQ_PATT.search(fastqs[0]) 38 m = FASTQ_PATT.search(fastqs[0])
37 run_acc = 'sample' 39 run_acc = 'sample'
38 if m is not None: 40 if m is not None:
39 run_acc = m.group(1) 41 run_acc = m.group(1)
40 study_acc = run_acc 42 study_acc = 'study'
41 if 'study' in config: 43 if 'study' in config:
42 study_acc = config['study'] 44 study_acc = config['study']
43 genome = 'hg38' 45 genome = 'hg38'
44 if 'genome' in config: 46 if 'genome' in config:
45 genome = config['genome'] 47 genome = config['genome']