教訓(xùn)：
1、由于大小寫不分，導(dǎo)致劈分字符串得到空集，浪費(fèi)了很多時(shí)間?。?！
2、循環(huán)中忘記重置列表，導(dǎo)致列表重復(fù)！??！，浪費(fèi)了很多時(shí)間?。。?br> 知識點(diǎn)：
re.finditer的用法
f=re.finditer（‘a(chǎn)’，b）#在字符串b中尋找所有‘a(chǎn)’
for i in f：
print（f.start()） #打印每個(gè)‘a(chǎn)’出現(xiàn)的位置
如果尋找的字符串有多個(gè)字節(jié)，則可以print（f.span()）

'''
提取陽性位點(diǎn)思路：
1、把a(bǔ)cc提取出來，放進(jìn)一個(gè)列表
2、把position提取出來，放進(jìn)一個(gè)列表
3、把sequence提取出來，放進(jìn)一個(gè)列表
4、遍歷position列表，遍歷position列表元素，在sequence中定位position位置，提取上下7個(gè)
5、寫入文件
'''
acc=[]
position_=[]
sequence_=[]
import re
with open('fasta-r.txt','r') as f:
    for line in f:
        if re.match(">", line):                                   #若該行開頭匹配到>，提取acc和position
            acc.append(re.search('(>)(.*)(@)',line).group(2))
            position_.append(re.search('(@)(.*)(;)',line).group(2))
        else:                                                     #該行沒有匹配到>，提取sequence
            line=line.rstrip()
            sequence_.append(line)

position=[]
for i in position_:           #用;分割position_放進(jìn)position列表
    position.append(i.split(';'))

sequence=[]                  #在sequence首尾都加上7個(gè)*
for i in sequence_:
    d='*******'+i+'*******'
    sequence.append(d)

positive=[]
for i in range(len(position)):
    c=sequence[i]
    positive_=[]    #每次循環(huán)要重制放每個(gè)序列positive的列表?。。。。。?！
    for k in position[i]:
        s=c[int(k)-1:int(k)+14]
        positive_.append(s)
    positive.append(positive_)

with open('positive.txt','w+') as g:
    g.write('Positive位點(diǎn)：'+'\n')
    for i in range(len(acc)):
        g.write('acc:'+acc[i]+'\n')
        x=position[i]
        y=positive[i]
        for k in range(len(position[i])):
            g.write('position:'+' '+x[k]+' '+y[k]+'\n')


'''
提取陰性位點(diǎn)思路：
1、提取最初文件中的code列
2、在每個(gè)acc中，查找每個(gè)code位置，刪除陽性位置，放進(jìn)列表negative_position
'''
all_data=[]                            #提取最初文件中的code列和acc列
data=[]
with open('example.elm','r') as f:
    all_data=f.readlines()
    for i in range(len(all_data)):
        data1=[]
        for element in all_data[i].split():
            data1.append(str(element))
        data.append(data1)
code_=[]
acc_=[]
for i in data:
    code_.append(i[3])
    acc_.append(i[0])
from collections import defaultdict    #key=acc,value=code,一個(gè)acc對應(yīng)多個(gè)value，將acc去重
d = defaultdict(list)
for key, value in zip(acc_, code_):
    d[key].append(value)
code=[]
for key in d:
    code.append(d.get(key))            #code列表每個(gè)元素對應(yīng)每個(gè)acc，一個(gè)元素是對應(yīng)acc所有的磷酸化位點(diǎn)對應(yīng)的氨基酸類型
#發(fā)現(xiàn)都是k，那只要找sequence中所有除psition中的k對應(yīng)的k的上下7個(gè)

#先找到每個(gè)序列所有k的位置，然后刪除陽性位點(diǎn)
import re
n_position=[]
for h in range(len(sequence)):
    f=re.finditer('K',sequence[h])
    all_position=[]
    n_position1=[]
    for i in f:
        a=i.start()
        b=a-6
        all_position.append(str(b))              #所有K在sequence中的位置（索引+1）
    for g in all_position:
        if g not in position[h]:
            n_position1.append(g)
    n_position.append(n_position1) #刪除陽性位點(diǎn),得到每個(gè)序列陰性位點(diǎn)的列表

negative=[]
for i in range(len(n_position)):
    c=sequence[i]
    negative_=[]
    for k in n_position[i]:
        s=c[int(k)-1:int(k)+14]
        negative_.append(s)
    negative.append(negative_)

with open('negative.txt','w+') as q:
    q.write('Negative位點(diǎn)：'+'\n')
    for i in range(len(acc)):
        q.write('acc:'+acc[i]+'\n')
        x=n_position[i]
        y=negative[i]
        for k in range(len(n_position[i])):
            q.write('position:'+' '+x[k]+' '+y[k]+'\n')