This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Script/Application to create a genbank format sub-sequence

I previously asked this question here:

Creating a 'subsetted' genbank of an operon from a whole genome

And accepted an answer using CLC which did work for the couple that I was doing at the time.

I'm now wondering if anyone has a commandline script/app/pipeline that would work for this as I need to do a fair few more, and doing them by hand will be a less than pleasant experience.

I need to do this because I need all the locus tags/gene IDs to be preserved - but not really sure how to go about it?

genome
• 2,424 views • link
written by Joe Joe

1 answer

Wrote my own in the end! Fairly novice biopython-er so suggestions/feedback welcomed.

Here: https://github.com/jrjhealey/bioinfo-tools/blob/master/Genbank_slicer.py

• 0 views • link
written by Joe Joe

Log in to answer this question.

More posts like this