Run workflow snakemake -c10 -pr pick_ancestral_sequence The output sequence will be in output/ref_to_coord.fasta