Python--找基因

在这里插入代码片
```def test(str,list1=[]):
    position=str.find('ATG')

    if (len(str)<3) or (position==-1):
        return list1
    else:
        p1=str[position:].find('TAG')
        p2=str[position:].find('TAA')
        p3=str[position:].find('TGA')
        if p1!=-1 and p2!=-1 and p3!=-1:
            p1=p1 if p1<p2 else p2
            p3=p2 if p3<p2 else p2
            p1=p1 if p1<p3 else p3
        elif p1==-1 and p2!=-1 and p3!=-1:
            p1=p2 if p2<p3 else p3
        elif p2 == -1 and p1 != -1 and p3 != -1:
            p1=p1 if p1<p3 else p3
        elif p3 == -1 and p2 != -1 and p1 != -1:
            p1=p1 if p1<p2 else p2
        elif p1==-1 and p2==-1 and p3!=-1:
            p1=p3
        elif p1==-1 and p3==-1 and p2!=-1:
            p1=p2
        else:
            p1=p1
        p1+=position
        if len(str[position+3:p1])%3==0:
            list1.append(str[position+3:p1])
            return test(str[p1+3:],list1)
    return list1
strr='TTATGTTTTAAGGATGGGGCGTTAGTT'
str2='TGTGTGTATAT'
list=[]
llist=[]
mylist=test(strr,list)
mylist2=test(str2,llist)
print(mylist)

print(mylist2)

猜你喜欢

转载自blog.csdn.net/weixin_41508608/article/details/89363310