## script to reverse orientation of a sequence, including all feature annotations from sys import argv from backbonomist.libs.loaders import load_genbank from backbonomist.libs.writers import write_genbank infile = "data/" + argv[1] outfile = "data/" + argv[2] record = load_genbank(infile) # make a genbank file of the contig new_record = record.reverse_complement() write_genbank(outfile, new_record)
## script to reverse orientation of a sequence, including all feature annotations import re from sys import argv from backbonomist.libs.loaders import load_genbank from backbonomist.libs.writers import write_genbank infile = "data/"+argv[1] outfile = "data/"+argv[2] record = load_fasta(origin_dir+"/"+filename) # make a genbank file of the contig gbk_file = "".join([destin_dir, rec_name, ".gbk"]) record.name = rec_name record.id = rec_name record.seq.alphabet = generic_dna write_genbank(gbk_file, record)
## script to reverse orientation of a sequence, including all feature annotations from sys import argv from backbonomist.libs.loaders import load_genbank from backbonomist.libs.writers import write_genbank infile = "data/"+argv[1] outfile = "data/"+argv[2] record = load_genbank(infile) # make a genbank file of the contig new_record = record.reverse_complement() write_genbank(outfile, new_record)