在这里插入代码片
```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)
Python--找基因
最新推荐文章于 2024-08-17 21:30:16 发布