Skip to content

Commit

Permalink
Use augur curate to generate intermediate Nextclade metadata file [#21]
Browse files Browse the repository at this point in the history
  • Loading branch information
genehack committed Dec 20, 2024
1 parent 2239f4a commit 0d6acd2
Showing 1 changed file with 25 additions and 0 deletions.
25 changes: 25 additions & 0 deletions ingest/rules/nextclade.smk
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,31 @@ rule run_nextclade:
"""


rule nextclade_metadata:
input:
nextclade="results/nextclade.tsv",
output:
nextclade_metadata=temp("results/nextclade_metadata.tsv"),
params:
nextclade_id_field=config["nextclade"]["id_field"],
nextclade_field_map=[f"{old}={new}" for old, new in config["nextclade"]["field_map"].items()],
nextclade_fields=",".join(config["nextclade"]["field_map"].values()),
log:
"logs/nextclade_metadata.txt",
benchmark:
"benchmarks/nextclade_metadata.tsv",
shell:
r"""
augur curate rename \
--metadata {input.nextclade:q} \
--id-column {params.nextclade_id_field:q} \
--field-map {params.nextclade_field_map:q} \
--output-metadata - \
| csvtk cut --tabs --fields {params.nextclade_fields:q} \
> {output.nextclade_metadata:q} 2> {log:q}
"""


rule join_metadata_and_nextclade:
input:
nextclade="results/nextclade.tsv",
Expand Down

0 comments on commit 0d6acd2

Please sign in to comment.