Extract Certain Info From Multiple Files In Genbank Format

Question

[Redse]

hi,

i'm new in bioinformatics and python. my first project is to extract /gene, /product, /locustag and /ECnumber qualifiers info from Genbank file. i need the info for all these qualifiers only and save them in 1 specific file. is there any specific script to do this? thanks..

[Redse]

Hi,

thanks guys...will look into those links right now..

Answer 1

[Lars Juhl Jensen]

This should be easily doable using BioPython. I have not used it much myself, but take a look at the examples here.

Answer 2

[Pierre Lindenbaum]

I found tutorial about biopython and Genbank at http://www2.warwick.ac.uk/fac/sci/moac/students/peter_cock/python/genbank/

[Lars Juhl Jensen]

Ha, I think we posted answers with the exact same links within seconds of each other :)

[Pierre Lindenbaum]

I was the first one ! ;-) I was the first one ! I was the first one ! Vote for me ! Vote for me ! Vote for me ! :-))

Answer 3

[Paulo Nuin]

You can definitely use the BioPython approach, but if you want to take a look a "raw" python way check this out.

You can modify the code to fit your needs.