This issue replaces #32, and is coded up in #179. However, my current implementation is overly complex because SampleData files have fixed ancestral states (i.e. states at the root), which are not currently accounted for in map_mutations.
Note that it should be possible to change the (future) Sankoff algorithm to allow for known ancestral states, but this may risk bloating the algorithm code for this specific corner case.