突变氨基酸位点

D:\1CAAS\Lab\songqianlin\Cas新蛋白\Cas12\cas12_lmnopq_fasta\重复项去除\突变氨基酸.py

点击查看代码

# amino_acid_mutation.py

seq = input("请输入蛋白序列：").strip().upper()

pos = int(input("请输入突变位点（从1开始）："))

new_aa = input("请输入突变后的氨基酸：").strip().upper()

# 检查位置是否合法
if pos < 1 or pos > len(seq):
    print("错误：位点超出序列长度")
    exit()

# 检查氨基酸是否合法
aa_set = set("ACDEFGHIKLMNPQRSTVWY")
if new_aa not in aa_set:
    print("错误：请输入标准20种氨基酸单字母代码")
    exit()

wt_aa = seq[pos - 1]

mut_seq = seq[:pos - 1] + new_aa + seq[pos:]

print("\n===== 结果 =====")
print(f"原氨基酸: {wt_aa}")
print(f"突变位点: {pos}")
print(f"突变后氨基酸: {new_aa}")
print(f"突变名称: {wt_aa}{pos}{new_aa}")

print("\n突变后序列:")
print(mut_seq)