Mercurial > repos > yating-l > rename_scaffolds
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]) |