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]) |
