我正在尝试设置一个脚本,它将根据我在NGS数据中找到的变体将密码子序列转换为另一个密码子序列。目前,我的脚本创建了一个由制表符分隔的输出文件,包含6列.每一栏代表以下内容:2:核苷酸碱基第四,密码子中基因组的位置顺序。0.0346505 A
是否可以在输出文件中设置第7列,根据第4列的数字将密码子的核苷酸序列更改为例如,密码子<
我写了一个小循环,一次浏览一个密码子站点。从这里我生成了一个字典,它告诉我每个密码子的位置都有哪些密码子,代码看起来像这样: for i in range(0, len(sequencesCombined[0]), 3):
codons = {任何随机站点的输出示例可能如下所示: 'ATG':5 这告诉我在5个序列中有5个ATG密码子。或者,它可能看起来像这样: 'CCT':2, 'CAA':3 因此,这是一种非同义密码子替换,因为两个序列表达编