【Python小试】计算蛋白序列中指定氨基酸所占的比例（可以根据蛋白质的氨基酸序列计算出）

【Python小试】计算蛋白序列中指定氨基酸所占的比例（可以根据蛋白质的氨基酸序列计算出）

编码

from __future__ import divisiondef get_aa_percentage(protein, aa_list=['A','I','L','M','F','W','Y','V']):    protein = protein.upper()    protein_length = len(protein)    total = 0    for aa in aa_list:        aa = aa.upper()        aa_count = protein.count(aa)        total += aa_count    percentage = total * 100 / protein_length    return percentageassert get_aa_percentage("MSRSLLLRFLLFLLLLPPLP", ["M"]) == 5assert get_aa_percentage("MSRSLLLRFLLFLLLLPPLP", ['M', 'L']) == 55assert get_aa_percentage("MSRSLLLRFLLFLLLLPPLP", ['F', 'S', 'L']) == 70assert get_aa_percentage("MSRSLLLRFLLFLLLLPPLP") == 65

解析

Python assert（断言）用于判断一个表达式，在表达式条件为 False 的时候触发异常。 断言可以在条件不满足程序运行的情况下直接返回错误，而不必等待程序运行后出现崩溃的情况。

assert expression#等价于if not expression:    raise AssertionError

Ref：​​https://runoob.com/python3/python3-assert.html​​

作者：Bioinfarmer

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。