0

我尝试从标头中具有增量 ID 的 .csv 文件创建多 fasta。但是我的脚本效果不佳。有什么建议吗?

M=open('/home/anna/Scrivania/db_mat/matk_2db_c.csv','r')
M_out=open('/home/anna/Scrivania/db_mat/db_matk_bronx.fas', 'w')

for i in range[1,92]:
    num = "m" + str(i)
for line in M:
    line1=line.split(',')
    ID = line1[0]
    SEQ = line1[1]
    seq =line1[2]
M_out.write('>'+ num +'_' +ID +'_'+ SEQ +'\n' + seq )

M.close()
M_out.close()
4

1 回答 1

1
for i,line in enumerate(M,1):
    num = "m"+str(i)
    line1=line.split(',')
    ID = line1[0]
    SEQ = line1[1]
    seq =line1[2]
    M_out.write(...)

这只是枚举行并为您提供 1..N 的索引,其中 N 是行数

在你迭代所有的数字之前,所以当你到达for line in Mnum 时,所有行都等于“m91”

于 2013-06-18T16:07:25.683 回答