adding public function for creating identity substitution matrices

this is immediately getting a deprecation warning however, as we will soon be adding a SubstitutionMatrix class that will include this functionality but we need the functionality in the meantime for a few special cases.

fixes #535