对.txt文件,经过处理提取之后,存为.txt文件,可以采用以下代码:
file_1 = '/cancer_fasta_PA.txt' #读取的.txt文件
file_3 = '/cancerPA.txt' #要保存的.txt文件
# positive seq protein A
pos_seq_protein_A = read_file(file_1) #这里的读取是前一个文章写的
outfile1=open(file_3,'w')
#方法一
for line in pos_seq_protein_A:
outfile1.write(line+'\n')
outfile1.close
#方法二
for i in range(len(pos_seq_protein_A)):
outfile1.write(pos_seq_protein_A[i]+'\n')
if i>=91:
print(pos_seq_protein_A[i])
outfile1.close
方法1存在一个bug,我怎么看怎么对,需要大神指点,我的.txt里面是108个蛋白质序列,但是经过方法一只剩91个(0-90),但是line确实能输出到最后一个。
于是我就采用方法二,每一个都输出,这样总能进去吧,但是依然不能,最后,我只能输出91个之外的,就得过且过粘贴了。
但我依然觉得代码没错,哪位大神路过指点一上下~
ps: 今天的科研不开心,搞了点啥啊!难过……