[Solved] how to write in python without using Biopython package
I recommend using biopython from Bio import SeqIO file = “file.gb” #gb = next(SeqIO.parse(open(file), “genbank”)) in python 3 gb = SeqIO.parse(open(file), “gb”).next() phosphorylation_list = [f for f in gb.features if f.type==”Site” and “phosphorylation” in f.qualifiers[‘site_type’]] for f in phosphorylation_list: print((int(f.location.start), int(f.location.end))) you get, (228, 229) (677, 678) (692, 693) (694, 695) (990, 991) (994, 995) … Read more