diff --git a/phylogenetic/rules/merge_sequences_usvi.smk b/phylogenetic/rules/merge_sequences_usvi.smk index ffc7a50..1b6c838 100644 --- a/phylogenetic/rules/merge_sequences_usvi.smk +++ b/phylogenetic/rules/merge_sequences_usvi.smk @@ -57,10 +57,10 @@ rule append_usvi: metadata = "data/metadata_all.tsv" shell: """ - seqkit rmdup {input.usvi_sequences} {input.sequences} > {output.sequences} - augur merge \ --metadata ingest={input.metadata} usvi={input.usvi_metadata} \ + --sequences ingest={input.sequences} usvi={input.usvi_sequences} --metadata-id-columns accession \ - --output-metadata {output.metadata} + --output-metadata {output.metadata} \ + --output-sequences {output.sequences} """ \ No newline at end of file