[Bioperl-l] concatenate two embl sequence files

jan aerts (RI) jan.aerts at bbsrc.ac.uk
Tue Jan 17 04:54:29 EST 2006


Hi all,

Does anyone know of an easy way to concatenate two sequences, including
recalculation of features positions of the second one? E.g.
  seq 1 = 100 bp
    feature A: 5..15
  seq 2 = 200 bp
    feature B: 20..30
  => concatenated sequence 3 = 300 bp
       feature A: 5..15
       feature B: 120..130  <<<<<<<<<<<

Annotations (features without range) should be transferred as well.

Of course, it must be possible to create a blank sequence and work my
way through all features, adding them to a new collection of features
and stuff. But I was wondering if a simpler technique is possible.

Many thanks,
Jan Aerts
Bioinformatics Department
Roslin Institute
Roslin, Scotland, UK

---------The obligatory disclaimer--------
The information contained in this e-mail (including any attachments) is
confidential and is intended for the use of the addressee only.   The
opinions expressed within this e-mail (including any attachments) are
the opinions of the sender and do not necessarily constitute those of
Roslin Institute (Edinburgh) ("the Institute") unless specifically
stated by a sender who is duly authorised to do so on behalf of the
Institute. 



More information about the Bioperl-l mailing list