Skip to content

Commit

Permalink
Update Snakefile_step2_single
Browse files Browse the repository at this point in the history
Rm singleton output for single read AdapterRemoval
  • Loading branch information
kdillmcfarland authored Jun 26, 2024
1 parent bbf32d4 commit f6b702c
Showing 1 changed file with 2 additions and 3 deletions.
5 changes: 2 additions & 3 deletions Snakefile_step2_single
Original file line number Diff line number Diff line change
Expand Up @@ -32,7 +32,6 @@ rule adapterremoval:
R1 = lambda wildcards: config["SampleList"][wildcards.sample]["R1"]
output:
R1 = "result/1_trim/{sample}_R1_trim.fastq.gz",
single = "result/1_trim/{sample}_singleton.truncated.gz",
discard = "result/1_trim/{sample}.settings",
set = "result/1_trim/{sample}_discarded.gz"
params:
Expand All @@ -44,9 +43,9 @@ rule adapterremoval:

run:
if params.trimAdapt == True:
shell("AdapterRemoval --file1 {input.R1} --output1 {output.R1} --singleton {output.single} --discarded {output.discard} --settings {output.set} --threads {threads} --gzip --maxns 1 --minlength 15 --trimqualities --minquality 30 --trim5p {params.trim5p} --adapter1 {params.adapter1}")
shell("AdapterRemoval --file1 {input.R1} --output1 {output.R1} --discarded {output.discard} --settings {output.set} --threads {threads} --gzip --maxns 1 --minlength 15 --trimqualities --minquality 30 --trim5p {params.trim5p} --adapter1 {params.adapter1}")
if params.trimAdapt == False:
shell("AdapterRemoval --file1 {input.R1} --output1 {output.R1} --singleton {output.single} --discarded {output.discard} --settings {output.set} --threads {threads} --gzip --maxns 1 --minlength 15 --trimqualities --minquality 30 --trim5p {params.trim5p}")
shell("AdapterRemoval --file1 {input.R1} --output1 {output.R1} --discarded {output.discard} --settings {output.set} --threads {threads} --gzip --maxns 1 --minlength 15 --trimqualities --minquality 30 --trim5p {params.trim5p}")
if params.trimAdapt != True and params.trimAdapt != False:
print("Error: Config trimAdapt parameter must be set to True or False.")

Expand Down

0 comments on commit f6b702c

Please sign in to comment.