Python - Rosalind 开放阅读框架问题

如何解决Python - Rosalind 开放阅读框架问题

一个关于 Rosalind 的开放阅读框架练习,我得到的结果与示例任务中得到的结果不同。可以在 here 中找到练习说明。

我有这个代码

gencode = {"GCT": "A","GCC": "A","GCA": "A","GCG": "A","TGT": "C","TGC": "C","GAT": "D","GAC": "D","GAA": "E","GAG": "E","TTT": "F","TTC": "F","GGT": "G","GGC": "G","GGA": "G","GGG": "G","CAT": "H","CAC": "H","ATA": "I","ATT": "I","ATC": "I","AAA": "K","AAG": "K","TTA": "L","TTG": "L","CTT": "L","CTC": "L","CTA": "L","CTG": "L","ATG": "M","AAT": "N","AAC": "N","CCT": "P","CCC": "P","CCA": "P","CCG": "P","CAA": "Q","CAG": "Q","CGT": "R","CGC": "R","CGA": "R","CGG": "R","AGA": "R","AGG": "R","TCT": "S","TCC": "S","TCA": "S","TCG": "S","AGT": "S","AGC": "S","ACT": "T","ACC": "T","ACA": "T","ACG": "T","GTT": "V","GTC": "V","GTA": "V","GTG": "V","TGG": "W","TAT": "Y","TAC": "Y","TAA": "_","TAG": "_","TGA": "_"}

seq = 'AGCCATGTAGCTAACTCAGGTTACATGGGGATGACCCCGCGACTTGGATTAGAGTCTCTTTTGGAATAAGCCTGAATGATCCGAGTAGCATCTCAG'
rev_seq = seq[::-1]


def get_orf_proteins(seq):
    proteins=[]
    for i in range(len(seq)-2):
        if gencode[seq[i:i+3]] == 'M':
            print(i)
            prot = ''
            k = i
            while gencode[seq[k:k+3]] != '_' and k < len(seq)-3:
              prot += gencode[seq[k:k+3]]
              k += 3
            proteins.append(prot)
    return(list(set(proteins)))


print(get_orf_proteins(seq))
print(get_orf_proteins(rev_seq))

返回以下蛋白质序列:

['MGMTPRLGLEsllE','MTPRLGLEsllE','M','MIRVAS']
['MY','MSLVSPNKVFSEIRFSAPVGVHWTQSMY']

我是否遗漏了什么,或者示例解决方案不正确?

解决方法

DNA 串的反向补码不仅仅是反向的 DNA 串。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?
Java在半透明框架/面板/组件上重新绘画。
Java“ Class.forName()”和“ Class.forName()。newInstance()”之间有什么区别?
在此环境中不提供编译器。也许是在JRE而不是JDK上运行?
Java用相同的方法在一个类中实现两个接口。哪种接口方法被覆盖?
Java 什么是Runtime.getRuntime()。totalMemory()和freeMemory()?
java.library.path中的java.lang.UnsatisfiedLinkError否*****。dll
JavaFX“位置是必需的。” 即使在同一包装中
Java 导入两个具有相同名称的类。怎么处理?
Java 是否应该在HttpServletResponse.getOutputStream()/。getWriter()上调用.close()?
Java RegEx元字符(。)和普通点?