大家好,我想检查氨基酸序列的分值,我为其编写了以下代码,效果非常好。但问题是每次我都要编辑文件。我能从命令界面给出氨基酸序列吗。
AA_seq='AVTLSPQRS' # this is the input variable
sum=0
value={"V": 3.1,"Y":3.5,"W":4.7,"T" :5.3,"S":5.1,"P":3.7,
"F":4.7,"M":1.5,"K":8.9,"L":6,"I":4.3,"H":3.3,"G":7.1,
"E":7,"Q":5.4,"C":0.6,"D":7.6,"N":6,"R":8.7,"A":3.4}
print("Total length of sequence is:", len(AA_seq))
for i in AA_seq:
sum+=value[i]
print("Total Score is :", sum)
