Skip to content
Snippets Groups Projects
Commit 56e57773 authored by Tobias Marschall's avatar Tobias Marschall
Browse files

Update summary_statistics to use --merged-file option

parent 889bfebc
No related branches found
No related tags found
No related merge requests found
...@@ -1046,12 +1046,13 @@ rule summary_statistics: ...@@ -1046,12 +1046,13 @@ rule summary_statistics:
input: input:
segmentation = 'segmentation2/{sample}/{windows}.{bpdens}.txt', segmentation = 'segmentation2/{sample}/{windows}.{bpdens}.txt',
strandstates = 'strand_states/{sample}/{windows}.{bpdens}/intitial_strand_state', strandstates = 'strand_states/{sample}/{windows}.{bpdens}/intitial_strand_state',
sv_calls = 'sv_calls/{sample}/{windows}.{bpdens}/{method}.txt', sv_calls = 'sv_calls/{sample}/{windows}.{bpdens}/{method}_filter{filter}.txt',
complex = "sv_calls/{sample}/{windows}.{bpdens}/{method}.complex.tsv", complex = "sv_calls/{sample}/{windows}.{bpdens}/{method}_filter{filter}.complex.tsv",
merged = "postprocessing/merge/{sample}/{windows}.{bpdens}/{method}.txt",
output: output:
tsv = 'stats/{sample}/{windows}.{bpdens,selected_j[0-9\\.]+_s[0-9\\.]+_scedist[0-9\\.]+}/{method}.tsv', tsv = 'stats/{sample}/{windows}.{bpdens,selected_j[0-9\\.]+_s[0-9\\.]+_scedist[0-9\\.]+}/{method}_filter{filter,(TRUE|FALSE)}.tsv',
log: log:
'log/summary_statistics/{sample}/{windows}.{bpdens}/{method}.log' 'log/summary_statistics/{sample}/{windows}.{bpdens}/{method}_filter{filter}.log'
run: run:
p = [] p = []
try: try:
...@@ -1068,6 +1069,9 @@ rule summary_statistics: ...@@ -1068,6 +1069,9 @@ rule summary_statistics:
p.append(f) p.append(f)
except KeyError: except KeyError:
pass pass
if wildcards.filter == 'TRUE':
p.append('--merged-file')
p.append(input.merged)
additional_params = ' '.join(p) additional_params = ' '.join(p)
shell('utils/callset_summary_stats.py --segmentation {input.segmentation} --strandstates {input.strandstates} --complex-regions {input.complex} {additional_params} {input.sv_calls} > {output.tsv} 2> {log}') shell('utils/callset_summary_stats.py --segmentation {input.segmentation} --strandstates {input.strandstates} --complex-regions {input.complex} {additional_params} {input.sv_calls} > {output.tsv} 2> {log}')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment