GCG V. 9.0-GAP

Francois JEANMOUGIN pingouin at crystal.u-strasbg.fr
Wed May 28 08:40:33 EST 1997


In article <5mgpf6$ivr at lyra.csx.cam.ac.uk>,
	tjrc1 at mole.bio.cam.ac.uk (Tim Cutts) writes:
[...]
> I don't think it's a bug, per se, more a mis-feature.

	Sorry, but if I implement such a "mis-feature" in Clustal,
you will mail me it's a bug ;-)...Also, if the order of the 
sequences matter, it is also a corruption of the gap algorithm.

>  By default, in
> GCG9, gap (and pileup too) does not penalise inserting gaps at the
> ends of the sequences.  Consequently putting huge gaps in at the ends
> like this appears to be a good match (duh!).  If you use the -endw
> option as well, which enables penalisation of gaps inserted at the
> ends, you might get a more sensible answer.

	I will probably (after some tests) alias gap to gap -endw
for my users. This will save some time for me. 
	Even with the non-weighting gaps, the alignment shown 
can't be the highest scoring alignment, or a representant of
the highest scoring family of alignments (look at the quality!).

			There is a bug, sure.
						Francois.
-- 
Francois Jeanmougin     | groupe de bioinformatique / bioinformatics groupe
tel:(+33) 3 88 65 32 71 | IGBMC BP 163 67404 Illkirch France
e-mail : jeanmougin at igbmc.u-strasbg.fr
"C'est pas parcequ'on monte au banc, qu'il faut descendre a jeun."(Thiefaine)




More information about the Bio-soft mailing list