Biostar Beta. Not for public use.
Trimming local GB file with Biopython
0
Entering edit mode
3.2 years ago
13dsc • 0

I was wondering if I could use Biopython in order to trim a local gb file.

Ex. I have a 100kb gb file but I am interested in a continuous 25kb portion. I know the location of the start and stop of the 25kb region. How do I create a gb file which contains only that 25kb portion from the larger file.

ADD COMMENTlink
1
Entering edit mode

You can. Look into into SeqIO, and string slicing.

ADD REPLYlink
0
Entering edit mode

will that keep the gene annotations associated with the gb file for that slice?

ADD REPLYlink
0
Entering edit mode

Yes, just edit the seq record, and output the Seq object to a new file.

ADD REPLYlink
2
Entering edit mode
9 months ago
Joe 12k
United Kingdom

I wrote a script that does exactly this a while back.

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

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1