Hi Roi You could try cutseq to remove the middle bit of the sequence: http://emboss.sourceforge.net/apps/cvs/emboss/apps/cutseq.html
extractseq might also be worth a look: http://emboss.sourceforge.net/apps/cvs/emboss/apps/extractseq.html Have a play and let us know how you get on. Cheers Jon > > Hello, > > > > Suppose I have some circular DNA RefSeq in a Genbank format. I would > like to create a new (smaller) Genbank file (both features and sequence) > containing only a range of nucleotides from the original file, but a > range that overlaps the "end". For example, if my original Genbank spans > 1,000,000 base (1-1,000,000) I would like to get join(800,000-1,000,000, > 1-100,000), so the new Genbank file will actually contain 300,000 bp > from the circular DNA. > > > > How can I do that? Again, my input is a Genbank file and output should > also be a Gebnank file. > > > > Thanks! > > - Roi > _________________________________________________________________ > Hotmail is redefining busy with tools for the New Busy. Get more from your > inbox. > http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_2 > _______________________________________________ > EMBOSS mailing list > EMBOSS@lists.open-bio.org > http://lists.open-bio.org/mailman/listinfo/emboss > _______________________________________________ EMBOSS mailing list EMBOSS@lists.open-bio.org http://lists.open-bio.org/mailman/listinfo/emboss