This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Multi-records genbank to CDS

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 enter image description here
And the result looks like this enter image description here


How can i do this?

cds genbank biopython
• 2,093 views • link
written by hazirliver • hazirliver

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.

• 0 views • link
written by husensofteng husensofteng

1 answer

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/

• 0 views • link
written by Joe Joe

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.

• 0 views • link
written by hazirliver • hazirliver

Log in to answer this question.

More posts like this