Multi-records genbank to CDS

Question

[hazirliver]

Hi! I have a file containing several genebank records written one after the other. I need to extract CDS (protein sequnce(/translation), /locus_tag, /inference, /product and contig id) from all contigs. How can i do it?
The input format looks like this
And the result looks like this

How can i do this?

[husensofteng]

Since you are analyzing data, it would be helpful if you make some effort to write a small script to read a file line by line and process it.

Answer 1

[Joe]

To clarify, you want all proteins/products, from all the entries in the file?

If so, take a look here: https://warwick.ac.uk/fac/sci/moac/people/students/peter_cock/python/genbank2fasta/

[hazirliver]

Yes, thanks! The code in this article was giving me an error, but this article got me on the right way to find the answer. I found the right solution using SeqIO.InsdcIO.GenBankCdsFeatureIterator.