Skip to content
Snippets Groups Projects
Verified Commit c3fd3dfa authored by Mia Croiset's avatar Mia Croiset
Browse files

hicstuff fragment & filter plot option correction

parent a757671e
No related branches found
No related tags found
No related merge requests found
......@@ -538,13 +538,20 @@ if __name__ == "__main__":
prefix = args.prefix
pairs_out = prefix+"_"+pairs_out
plot = args.plot
plot.title()
if plot == "False":
if plot.lower() == "false":
plot = False
elif plot == "True":
elif plot.lower() == "true":
plot = True
dist_plot = args.dist_plot
if dist_plot == "None":
dist_plot = None
else:
dist_plot = prefix+"_"+dist_plot
pie_plot = args.pie_plot
if pie_plot == "None":
pie_plot = None
else:
pie_plot = prefix+"_"+pie_plot
uncut_thr, loop_thr = get_thresholds(
pairs_idx, plot_events=plot, fig_path=dist_plot, prefix=prefix
......
......@@ -19,9 +19,17 @@ if __name__ == "__main__":
index = args.index
min_size = args.min_size
circular = args.circular
if circular.lower() == "false":
circular = False
elif circular.lower() == "true":
circular = True
output_contigs = args.output_contigs
output_frags = args.output_frags
plot = args.plot
if plot.lower() == "false":
plot = False
elif plot.lower() == "true":
plot = True
fig_path = args.fig_path
......
......@@ -121,10 +121,10 @@ params {
// Hicstuff
hicstuff_bwt2_align_opts = '--very-sensitive-local'
hicstuff_min_size = 0
hicstuff_circular = 'False'
hicstuff_circular = 'frue'
hicstuff_output_contigs = 'info_contigs.txt'
hicstuff_output_frags = 'fragments_list.txt'
hicstuff_frags_plot = 'False'
hicstuff_frags_plot = 'false'
hicstuff_frags_plot_path = 'frags_hist.pdf'
hicstuff_valid_pairs = 'valid.pairs'
hicstuff_valid_idx = 'valid_idx.pairs'
......@@ -132,9 +132,9 @@ params {
hicstuff_matrix = 'abs_fragments_contacts_weighted.txt'
hicstuff_bin = 10000
hicstuff_valid_idx_filtered = 'valid_idx_filtered.pairs'
hicstuff_plot_events = 'False'
hicstuff_pie_plot = 'None'
hicstuff_dist_plot = 'None'
hicstuff_plot_events = 'false'
hicstuff_pie_plot = 'distrib'
hicstuff_dist_plot = 'dist'
//Hicstuff optional modules
filter_event = true
......
......@@ -10,6 +10,7 @@ process FILTER_EVENT {
output:
tuple val(meta1), path("*_valid_idx_filtered.pairs"), emit: idx_pairs_filtered
path("*.png"), optional: true
script:
......
......@@ -12,6 +12,7 @@ process FRAGMENT_ENZYME {
output:
tuple val(meta), path(params.hicstuff_output_contigs) , emit: info_contigs
tuple val(meta), path(params.hicstuff_output_frags), emit: fragments_list
path("*.pdf"), optional: true
script:
......
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