Skip to content
Snippets Groups Projects
Commit 3ca02902 authored by Sascha Meiers's avatar Sascha Meiers
Browse files

Use more recent SV plot script + do not black list counts based initially

parent 8916a913
No related branches found
No related tags found
No related merge requests found
...@@ -5,7 +5,7 @@ ...@@ -5,7 +5,7 @@
"mosaicatcher" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/build/mosaicatcher", "mosaicatcher" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/build/mosaicatcher",
"plot_script" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/R/qc.R", "plot_script" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/R/qc.R",
"plot_segments" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher/R/plot_segments.R", "plot_segments" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher/R/plot_segments.R",
"sv_plot_script": "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher/R/svs.R", "sv_plot_script": "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/R/chrom.R",
"norm_script" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/R/norm.R", "norm_script" : "/g/korbel/meiers/tools/mosaicatcher/mosaicatcher_experimental/R/norm.R",
"samtools" : "samtools", "samtools" : "samtools",
"bcftools" : "bcftools", "bcftools" : "bcftools",
......
...@@ -179,9 +179,11 @@ rule plot_SV_calls: ...@@ -179,9 +179,11 @@ rule plot_SV_calls:
segments = "segmentation2/{sample}/{windows}.{bpdens}.txt" segments = "segmentation2/{sample}/{windows}.{bpdens}.txt"
output: output:
"sv_calls/{sample}/{windows}.{bpdens}/{method}.{chrom}.pdf" "sv_calls/{sample}/{windows}.{bpdens}/{method}.{chrom}.pdf"
params:
sv_plot_script = config["sv_plot_script"]
shell: shell:
""" """
Rscript utils/chrom.R segments={input.segments} strand={input.strand} calls={input.calls} {input.counts} {wildcards.chrom} {output} Rscript {params.sv_plot_script} segments={input.segments} strand={input.strand} calls={input.calls} {input.counts} {wildcards.chrom} {output}
""" """
rule plot_SV_calls_simulated: rule plot_SV_calls_simulated:
...@@ -193,9 +195,11 @@ rule plot_SV_calls_simulated: ...@@ -193,9 +195,11 @@ rule plot_SV_calls_simulated:
truth = "simulation/variants/genome{seed}-{window}.txt" truth = "simulation/variants/genome{seed}-{window}.txt"
output: output:
"sv_calls/simulation{seed}-{window}/{window}_fixed.{bpdens}/{method}.{chrom}.pdf" "sv_calls/simulation{seed}-{window}/{window}_fixed.{bpdens}/{method}.{chrom}.pdf"
params:
sv_plot_script = config["sv_plot_script"]
shell: shell:
""" """
Rscript utils/chrom.R segments={input.segments} strand={input.strand} truth={input.truth} calls={input.calls} {input.counts} {wildcards.chrom} {output} Rscript {params.sv_plot_script} segments={input.segments} strand={input.strand} truth={input.truth} calls={input.calls} {input.counts} {wildcards.chrom} {output}
""" """
...@@ -248,6 +252,7 @@ rule mosaic_count_fixed: ...@@ -248,6 +252,7 @@ rule mosaic_count_fixed:
""" """
{params.mc_command} count \ {params.mc_command} count \
--verbose \ --verbose \
--do-not-blacklist-hmm \
-o {output.counts} \ -o {output.counts} \
-i {output.info} \ -i {output.info} \
-x {input.excl} \ -x {input.excl} \
......
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