[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