areagp61 at yahoo.it
Fri Jan 10 13:30:37 UTC 2003
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?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the EMBOSS