Commit d82790d4 authored by Rachel  LEGENDRE's avatar Rachel LEGENDRE
Browse files

include peaks from IDR in igv_session only if narrow mode is chosen

parent 5e495879
...@@ -970,8 +970,10 @@ if config["igv_session"]["do"]: ...@@ -970,8 +970,10 @@ if config["igv_session"]["do"]:
if len(IP_REP) > 1 and config["macs2"]["do"]: if len(IP_REP) > 1 and config["macs2"]["do"]:
if config["intersectionApproach"]["do"]: if config["intersectionApproach"]["do"]:
igv_session_input += expand(intersectionApproach_output, CALLER="macs2_" + model, IP_IDR=IP_REP) igv_session_input += expand(intersectionApproach_output, CALLER="macs2_" + model, IP_IDR=IP_REP)
else: elif model in ["narrow"] and not config["intersectionApproach"]["do"] and config["compute_idr"]["do"]:
igv_session_input += expand(select_peaks_output, CALLER="macs2", IP_IDR=IP_REP) igv_session_input += expand(select_peaks_output, CALLER="macs2", IP_IDR=IP_REP)
else :
igv_session_input += expand(macs2_output, SAMPLE=IP_ALL)
elif len(IP_REP) > 1 and config["seacr"]["do"]: elif len(IP_REP) > 1 and config["seacr"]["do"]:
igv_session_input += expand(intersectionApproach_output, CALLER="seacr_" + config["seacr"]["threshold"], IP_IDR=IP_REP) igv_session_input += expand(intersectionApproach_output, CALLER="seacr_" + config["seacr"]["threshold"], IP_IDR=IP_REP)
if len(IP_ALL) > 1 and config["macs2"]["do"]: if len(IP_ALL) > 1 and config["macs2"]["do"]:
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment