[BioPython] gapped consensus?
Sebastian Bassi
sbassi at asalup.org
Thu Jun 12 15:06:46 EDT 2003
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Sebastian Bassi wrote:
| I would like to have the chance to have "gapped consensus".
| Now the dumb_consensus (on Align.AlignInfo) work like this:
I reply yo myself in some extent. Commenting out this line do the work:
#if record.seq[n] != '-' and record.seq[n] != '.':
I didn't see the results at first because there was a AlignInfo.pyc file
and the Python interpreter didin't look for my modified AlignInfo.py file.
This was a quick hack, now I'm modifing the function in order to accept
as a parameter to display gaps or not. When I finish, I will submit it
for consideration of inclusion on the CVS.
- --
Best regards,
//=\ Sebastian Bassi - Diplomado en Ciencia y Tecnologia, UNQ //=\
\=// IT Manager Advanta Seeds - Balcarce Research Center - \=//
//=\ Pro secretario ASALUP - www.asalup.org - PGP key available //=\
\=// E-mail: sbassi at genesdigitales.com - ICQ UIN: 3356556 - \=//
~ http://Bioinformatica.info
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.0.6 (MingW32)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org
iEYEARECAAYFAj7osyYACgkQ6lc7ixf6gKp+5QCfSi6xUBnj38NkIJgeheVDgaBf
NnMAoIsaP56/+NQTI08FtZl2W9Xg5cJW
=ao1K
-----END PGP SIGNATURE-----
More information about the BioPython
mailing list