Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from Bio import SeqIO
- o = open("locus-tags.csv",'w')
- o.write("old-tag,new-tag\n")
- for seq in SeqIO.parse(open("file.gbk",'r'),'genbank'):
- for f in seq.features:
- if f.type == "gene":
- o.write("{},{}\n".format(f.qualifiers["old_locus_tag"][0],f.qualifiers["locus_tag"][0])
- o.close()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement