Iterative GAP with tabulated output???

Arle Kruckeberg arle at biovax.leeds.ac.uk
Sun Jan 22 15:05:57 EST 1995

Dear GCG whizzes:

I am studying a family of homologous proteins from diverse organisms. I
like to produce a table that includes % identity and similarity,
Needleman and Wunsch (GAP) alignment quality and the quality +-
standard deviation of x randomizations, and the number of gaps
introduced to produce each alignment. In other words, to run GAP lots
of times, on a list of file names, and put the numerical output into a
separate file (and leave the sequence alignment to PILEUP). Is this
easy, hard, already code in someone's directory?

Also, does anyone have a preferred method of converting standard
deviations (of x randomized sequence comparisons, as done in GAP/RAND)
to probabilities that gapped alignment of two sequences yields a
quality by chance?


a.k.a. Dr. Arthur L. Kruckeberg
Department of Biochemistry and Molecular Biology
University of Leeds
Leeds LS2 9JT
Great Britain
+44 532 333172
+44 532 333167 FAX
arle at biovax.leeds.ac.uk

More information about the Info-gcg mailing list

Send comments to us at biosci-help [At] net.bio.net