In article <443gdk$2gk at rc1.vub.ac.be> Tylzanowski Przemko <przemko> writes:
>I have to search genpept for a pattern. Well, everything worked, I got my
>findpatterns.find file and now I am stuck. I have about 140 hits in about as
>many genes. I would like to pull all the sequences out and have a quick look
>but I don't knwo how to do that. Of coursde, I could do it by hand, but it
>would be a bit tedious, especially if with some other patterns I find more
>hits.
>So, my question is this, how could I, easily and painlessly, retrieve the
>sequences and comments in question. Subsequently I would like to put them into
>pileup or smething like that.
The solution depends on what you want to do: extract all the complete
sequences, or extract all the hit regions (plus perhaps a few residues/bases
either side).
--
------------------------------------------------------------------------
Peter Rice | Informatics Division
E-mail: pmr at sanger.ac.uk | The Sanger Centre
Tel: (44) 1223 494967 | Hinxton Hall, Hinxton,
Fax: (44) 1223 494919 | Cambs, CB10 1RQ
URL: http://www.sanger.ac.uk/~pmr/ | England