Parsing clustal file for sequence variation analysis
import msaparser
parser = msaparser.Parser(blockstar=0.9, blocklen=10, checknum=5)
parser.parse('/path/to/protein.clustal', 'a')
parser.block # the list of parsed blocks
parser.mutnum # the number of point mutation
parser.mutprofile # the list of mutation profiles
parser.mutposlist # the list of mutation position
parser.get_clustal_html() # html output of parsed result