Skip to content

Commit

Permalink
fixup: smk spacing
Browse files Browse the repository at this point in the history
  • Loading branch information
j23414 committed Feb 12, 2024
1 parent e9a58ce commit 4f2591d
Showing 1 changed file with 25 additions and 25 deletions.
50 changes: 25 additions & 25 deletions ingest/workflow/snakemake_rules/split_serotypes.smk
Original file line number Diff line number Diff line change
Expand Up @@ -14,30 +14,30 @@ Parameters are expected to be defined in `config.transform`.
"""

rule split_by_ncbi_serotype:
"""
Split the data by serotype based on the NCBI metadata.
"""
input:
metadata = "data/metadata.tsv",
sequences = "data/sequences.fasta"
output:
metadata = "results/metadata_{serotype}.tsv",
sequences = "results/sequences_{serotype}.fasta"
params:
ncbi_serotype = lambda wildcards: wildcards.serotype,
id_field = config["transform"]["id_field"]
shell:
"""
if [[ "{params.ncbi_serotype}" == "all" ]]; then
cp {input.metadata} {output.metadata}
cp {input.sequences} {output.sequences}
else
augur filter \
--sequences {input.sequences} \
--metadata {input.metadata} \
--metadata-id-columns {params.id_field} \
--query "ncbi_serotype=='{params.ncbi_serotype}'" \
--output-sequences {output.sequences} \
--output-metadata {output.metadata}
fi
Split the data by serotype based on the NCBI metadata.
"""
input:
metadata = "data/metadata.tsv",
sequences = "data/sequences.fasta"
output:
metadata = "results/metadata_{serotype}.tsv",
sequences = "results/sequences_{serotype}.fasta"
params:
ncbi_serotype = lambda wildcards: wildcards.serotype,
id_field = config["transform"]["id_field"]
shell:
"""
if [[ "{params.ncbi_serotype}" == "all" ]]; then
cp {input.metadata} {output.metadata}
cp {input.sequences} {output.sequences}
else
augur filter \
--sequences {input.sequences} \
--metadata {input.metadata} \
--metadata-id-columns {params.id_field} \
--query "ncbi_serotype=='{params.ncbi_serotype}'" \
--output-sequences {output.sequences} \
--output-metadata {output.metadata}
fi
"""

0 comments on commit 4f2591d

Please sign in to comment.