[Bioperl-l] protal2dna and Bio::SimpleAlign
Jason Stajich
jason.stajich at duke.edu
Sun Jan 23 09:19:13 EST 2005
I'm not familiar with the script.
Bio::Align::Utilities does protein to DNA mapping for an alignment with
the aa_to_dna_aln function.
-jason
On Jan 22, 2005, at 4:07 PM, Maureen L Coleman wrote:
> Hi.
> I'm trying to use the protal2dna script (downloaded from Pasteur site)
> to convert protein alignments back to DNA alignments. It works in some
> cases but not in others. In the cases where it doesn't work, it pulls
> out the same sequence twice instead of pulling out seq1 and seq2 from
> my protein alignment. Then when it tries to match it up with the
> corresponding DNA sequence, it doesn't work - it matches prot1 with
> dna1 (correctly) and prot1 with dna2 (incorrectly).
>
> I suspect this might be related to the name,start,end (nse) method in
> Bio::SimpleAlign. Any suggestions?
>
> Thanks,
> Maureen
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
>
--
Jason Stajich
jason.stajich at duke.edu
http://www.duke.edu/~jes12/
More information about the Bioperl-l
mailing list