[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