diff --git a/Snakefile b/Snakefile index 36d36ea618db2b026d4542c27292ec46c8b2acd3..b79e7e10b66c4c614f560e041d210d54041faee7 100644 --- a/Snakefile +++ b/Snakefile @@ -326,6 +326,7 @@ rule segmentation: shell: """ {params.mc_command} segment \ + --remove-none \ -m 0.33 \ -M 50000000 \ -o {output} \ diff --git a/utils/helper.prepare_segments.R b/utils/helper.prepare_segments.R index 0ee3ae156b4dcc668ba522ea9240c94b8d64a446..fc92d8c69bc28f881bc12d7351d3b6fb1c7d4bc6 100755 --- a/utils/helper.prepare_segments.R +++ b/utils/helper.prepare_segments.R @@ -5,6 +5,6 @@ qu d = fread(snakemake@input[[1]]) d # type = 1 is important to get discrete values! -d = d[, .SD[k == quantile(1:max(k), qu, type = 1)], by=chrom][,.(k, chrom, bps = breakpoint)] +d = d[, .SD[k == quantile(1:max(k), qu, type = 1)], by=chrom][,.(k, chrom, bps)] d write.table(d, file = snakemake@output[[1]], row.names=F, quote=F, sep="\t")