add a function to trim to reference

[roblanf]

The idea here is that we want to trim aligned readsets according to some DNA reference sequence, before calling the consensus.

It's not totally obvious whether it's better to trim reads individually at the start, or wait and do it later. I suspect it doesn't make much difference, in which case we are better doing it wherever makes the most sense for users.