GCG V. 9.0-GAP
arve at nada.kth.se
Mon May 26 06:38:28 EST 1997
Iddo Friedberg writes:
IF> I've encountered an interesting phenomenon in GCG v. 9.0, UNIX.
I don't know much about GCG, but want to suggest an explanation anyways.
IF> Try the following:
IF> gap -in1=sw:lyg_ansan -in2=sw:kad3_bovin -gap=5 -len=3 -out=out1.pair
IF> gap -in1=sw:kad3_bovin -in2=sw:lyg_ansan -gap=5 -len=3 -out=out2.pair
IF> Now look at the results in files out1.pair and out2.pair
IF> Seems like the results are dependent upon sequence order. AFAIK, that
IF> shouldn't be the case in pairwise sequence alignment.
My guess that you are seeing to two different alignments, but with the
same score? If there is more than one alignment with the optimal score, one
alignment is chosen arbitrarily. There can be many equivalent (with
respect to alignment score) alignments.
Lars Arvestad (MIME correspondence welcome)
Dept. of Numerical Analysis and Computing Science
Royal Institute of Technology (KTH), Stockholm, Sweden
More information about the Bio-soft