Add mutation calculation functions in the enzy_htp/mutation module

[shaoqx]

Combine experience from KE_accre scripts about getting unique mutation and get relative mutation:

• calculate overall combination of mutations with given lists of mutaflag lists.
• calculate relative mutations with given a reference mutaflag list.

[shaoqx]

The operation are all targeting a list of mutaflag, maybe a variant class is useful