seqretsplit
Peter Rice
pmr at ebi.ac.uk
Fri Jan 10 14:36:53 UTC 2003
Ciao Graziono,
> 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"
This must be a GCG format database, indexed in SRS ... you have gaps in
the ".." format.
Ideally SRS would fix that ... but EMBOSS could cope with some small
changes. Anyway, an SRS fix is non-trivial because it is reporting
exactly the text that GCG stores.
> Anyone knows an application or a script that allows me to split
> sequences (embl format) in different files without losing FT lines?
Simplest in your case would be a script that changes ". ." to ".."
before passing the data to EMBOSS.
Meanwhile, I will take a look at fixing this for a future release.
regards,
Peter
More information about the EMBOSS
mailing list