comparison rename.py @ 1:148e353417e7 draft default tip

planemo upload commit af6c4dc7f6d6795f2b2db7299eb14ae1d5622ac1
author yating-l
date Fri, 30 Dec 2016 17:50:54 -0500
parents 1a0e39acb62f
children
comparison
equal deleted inserted replaced
0:1a0e39acb62f 1:148e353417e7
12 i = 1 12 i = 1
13 for line in lines: 13 for line in lines:
14 if ">" in line: 14 if ">" in line:
15 oldname = line[1:].rstrip() 15 oldname = line[1:].rstrip()
16 newname = "scaffold_" + str(i) 16 newname = "scaffold_" + str(i)
17 line = ">" + newname + "\n" 17 line = ">" + newname
18 i = i+1 18 i = i+1
19 namemap[oldname] = newname 19 namemap[oldname] = newname
20 out.write(line) 20 #TODO: Add line breaks to chromosome sequences that are in a single line
21 #else:
22 #if (len(line) > 50):
23 #for
24 out.write(line.rstrip() + "\n")
21 with open(indexfile, 'w') as index: 25 with open(indexfile, 'w') as index:
22 for k in namemap: 26 for k in namemap:
23 index.write(k + "=>" + namemap[k] + "\n") 27 index.write(k + "\t" + namemap[k] + "\n")
24 28
25 def main(): 29 def main():
26 inputfile = str(sys.argv[1]) 30 inputfile = str(sys.argv[1])
27 outputfile = str(sys.argv[2]) 31 outputfile = str(sys.argv[2])
28 indexfile = str(sys.argv[3]) 32 indexfile = str(sys.argv[3])