問題：假設(shè)存在一個(gè)蛋白序列seq 'MPKRKKMRLFEEDDEIRESLLGADNKDKDEDEDLQSDTENRTFLEDTDTV'闰渔，計(jì)算各個(gè)氨基酸出現(xiàn)的頻率虑灰。
方法1：
seq = 'MPKRKKMRLFEEDDEIRESLLGADNKDKDEDEDLQSDTENRTFLEDTDTV'
ani_result = {}
for ani in range(len(seq)):
if seq[ani] in ani_result:
ani_result[seq[ani]] += 1
else:
ani_result[seq[ani]] = 1
for each in sorted(ani_result.keys()):
print(each + " " + str(ani_result[each]))
#或者是列表格式
print(sorted(ani_result.items(), key=lambda x: x[1]))
運(yùn)行結(jié)果
方法2：利用count函數(shù)，統(tǒng)計(jì)20種氨基酸在蛋白序列中出現(xiàn)的次數(shù)
seq = 'MPKRKKMRLFEEDDEIRESLLGADNKDKDEDEDLQSDTENRTFLEDTDTV'
for ani in 'ARNDCQEGHIKMFPSTWYV':
num_result = seq.count(ani)
print(ani, num_result)
