Skip to content

Commit

Permalink
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion phylogenetic/rules/annotate_phylogeny.smk
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,8 @@ rule ancestral:
"""Reconstructing ancestral sequences and mutations"""
input:
tree = "results/{gene}/tree_{serotype}.nwk",
alignment = "results/{gene}/aligned_{serotype}.fasta"
alignment = "results/{gene}/aligned_{serotype}.fasta",
root_sequence = lambda wildcard: "config/reference_{serotype}_genome.gb" if wildcard.gene in ['genome'] else "results/config/reference_{serotype}_{gene}.gb",
output:
node_data = "results/{gene}/nt-muts_{serotype}.json"
params:
Expand All @@ -36,6 +37,7 @@ rule ancestral:
--tree {input.tree} \
--alignment {input.alignment} \
--output-node-data {output.node_data} \
--root-sequence {input.root_sequence} \
--inference {params.inference}
"""

Expand Down

0 comments on commit 9431655

Please # to comment.