seqretsplit

[Graziano P.]

Hi all,
I have a single file containing 164 sequences in embl format retrieved by SRS. I want to split all the entries in distinct files. I have tried to use "seqretsplit", but it do not report FT lines. For example:
My entry retrieved by means of SRS is:

ID   AY045754_4; parent: AY045754AC   AY045754;
FT   rRNA            610. .772
FT                   /product="5.8S ribosomal RNA"
SQ   Sequence   163 BP;
     TATATTCTAG CCTTATCGGT GGATCACTCG GCTCGTGGAT CGATGAAGAA CGCAGCCAAC        60
     TGCGATAATT AGTGTGAACT GCAGAAACCT TGAACACTGA ACTTTCGAAT GCACATTGCG       120
     CCATTGGAGT TATATCCATT GCACGCCTGG TTCAGGGTCG TAA                         163
//

but seqretsplit convert it in:

ID   AY045754_4 standard; DNA; UNC; 163 BP.
SQ   Sequence 163 BP; 43 A; 37 C; 39 G; 44 T; 0 other;
     TATATTCTAG CCTTATCGGT GGATCACTCG GCTCGTGGAT CGATGAAGAA CGCAGCCAAC        6
     TGCGATAATT AGTGTGAACT GCAGAAACCT TGAACACTGA ACTTTCGAAT GCACATTGCG       12
     CCATTGGAGT TATATCCATT GCACGCCTGG TTCAGGGTCG TAA                         16
//

Anyone knows an application or a script that allows me to split sequences (embl format) in different files without losing FT lines?

Best regards

Graziano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.open-bio.org/pipermail/emboss/attachments/20030110/952875b3/attachment-0001.html>