Skip to content

Commit

Permalink
change the output name of cnv modules
Browse files Browse the repository at this point in the history
  • Loading branch information
zzgeng committed Jul 17, 2023
1 parent 6676cf8 commit 3bde20f
Show file tree
Hide file tree
Showing 16 changed files with 515,811 additions and 247,331 deletions.
Empty file modified analyses/copy_number_consensus_call/README.md
100644 → 100755
Empty file.
16 changes: 8 additions & 8 deletions analyses/copy_number_consensus_call/Snakefile
100644 → 100755
Original file line number Diff line number Diff line change
Expand Up @@ -8,9 +8,9 @@ wildcard_constraints:
## Define the first rule of the Snakefile. This rule determines what the final file is and which steps to be taken.
rule all:
input:
"results/cnv_consensus.tsv",
"results/cnv_neutral.bed",
"results/cnv-consensus.seg.gz"
"results/cnv_consensus-freec-ckit-gatk.tsv",
"results/cnv_neutral-freec-ckit-gatk.bed",
"results/cnv-consensus-freec-ckit-gatk.seg.gz"


scratch_loc = config["scratch"] + '/'
Expand Down Expand Up @@ -290,7 +290,7 @@ rule merge_neutral:
expand(scratch_loc + "endpoints/{sample}.neutral.bed",
sample = config["samples"])
output:
"results/cnv_neutral.bed"
"results/cnv_neutral-freec-ckit-gatk.bed"
shell:
"cat {input} > {output}"

Expand All @@ -315,18 +315,18 @@ rule clean_output:
script=os.path.join(config["scripts"], "remove_dup_NULL_overlap_entries.py"),
cnv_file = scratch_loc + "endpoints/all_CNVs_combined.tsv"
output:
"results/cnv_consensus.tsv"
"results/cnv_consensus-freec-ckit-gatk.tsv"
shell:
"python3 {input.script} --file {input.cnv_file} > {output}"

rule make_segfile:
input:
script= os.path.join(config["scripts"], "bed_to_segfile.R"),
consensus = "results/cnv_consensus.tsv",
neutral = "results/cnv_neutral.bed",
consensus = "results/cnv_consensus-freec-ckit-gatk.tsv",
neutral = "results/cnv_neutral-freec-ckit-gatk.bed",
uncalled = "results/uncalled_samples.tsv"
output:
"results/cnv-consensus.seg.gz"
"results/cnv-consensus-freec-ckit-gatk.seg.gz"
shell:
"Rscript {input.script}"
" -i {input.consensus}"
Expand Down
Empty file modified analyses/copy_number_consensus_call/merge_cnv_file.R
100644 → 100755
Empty file.
Binary file not shown.
Empty file modified analyses/copy_number_consensus_call_manta/README.md
100644 → 100755
Empty file.
16 changes: 8 additions & 8 deletions analyses/copy_number_consensus_call_manta/Snakefile
100644 → 100755
Original file line number Diff line number Diff line change
Expand Up @@ -8,9 +8,9 @@ wildcard_constraints:
## Define the first rule of the Snakefile. This rule determines what the final file is and which steps to be taken.
rule all:
input:
"results/cnv_consensus.tsv",
"results/cnv_neutral.bed",
"results/cnv-consensus.seg.gz"
"results/cnv_consensus-freec-ckit-manta.tsv",
"results/cnv_neutral-freec-ckit-manta.bed",
"results/cnv-consensus-freec-ckit-manta.seg.gz"


scratch_loc = config["scratch"] + '/'
Expand Down Expand Up @@ -287,7 +287,7 @@ rule merge_neutral:
expand(scratch_loc + "endpoints/{sample}.neutral.bed",
sample = config["samples"])
output:
"results/cnv_neutral.bed"
"results/cnv_neutral-freec-ckit-manta.bed"
shell:
"cat {input} > {output}"

Expand All @@ -312,18 +312,18 @@ rule clean_output:
script=os.path.join(config["scripts"], "remove_dup_NULL_overlap_entries.py"),
cnv_file = scratch_loc + "endpoints/all_CNVs_combined.tsv"
output:
"results/cnv_consensus.tsv"
"results/cnv_consensus-freec-ckit-manta.tsv"
shell:
"python3 {input.script} --file {input.cnv_file} > {output}"

rule make_segfile:
input:
script= os.path.join(config["scripts"], "bed_to_segfile.R"),
consensus = "results/cnv_consensus.tsv",
neutral = "results/cnv_neutral.bed",
consensus = "results/cnv_consensus-freec-ckit-manta.tsv",
neutral = "results/cnv_neutral-freec-ckit-manta.bed",
uncalled = "results/uncalled_samples.tsv"
output:
"results/cnv-consensus.seg.gz"
"results/cnv-consensus-freec-ckit-manta.seg.gz"
shell:
"Rscript {input.script}"
" -i {input.consensus}"
Expand Down
Binary file not shown.
Binary file not shown.
Loading

0 comments on commit 3bde20f

Please sign in to comment.