提取子串模式

时间:2016-04-13 22:39:47

标签: python

我有1200个序列的长文件

>3fm8|A|A0JLQ2
CFLVNLNADPALNELLVYYLKEHTLIGSANSQDIQLCGMGILPEHCIIDITSEGQVMLTP
QKNTRTFVNGSSVSSPIQLHHGDRILWGNNHFFRLNLP


>2ht9|A|A0JLT0
LATAPVNQIQETISDNCVVIFSKTSCSYCTMAKKLFHDMNVNYKVVELDLLEYGNQFQDA
LYKMTGERTVPRIFVNGTFIGGATDTHRLHKEGKLLPLVHQCYL

我想阅读每个可能的模式在中间有半胱氨酸并且在开头有五个字符串并且跟随其他五个字符串如xxxxxCxxxxx

输出应该是这样的:

  • QDIQLCGMGIL
  • ILPEHCIIDIT
  • TISDNCVVIFS
  • FSKTSCSYCTM

这是pogram只给出C的位置。它不像我想要的那样工作

pos=[]

def find(ch,string1):

    for i in range(len(string1)):
        if ch == string1[i]:
            pos.append(i)
            return pos



z=find('C','AWERQRTCWERTYCTAAAACTTCTTT')

print z

0 个答案:

没有答案