[Bioperl-l] (no subject)
Chris Fields
cjfields at uiuc.edu
Mon Jul 21 16:42:37 UTC 2008
Ya beat me to it! I agree, this is the best way to go, and it would
save an extra BLAST run (though if you need it Diego's bl2seq
suggestion is the way to go).
chris
On Jul 21, 2008, at 11:31 AM, Dave Messina wrote:
> You know which protein sequence corresponds to which cDNA sequence,
> so I
> think you could also:
>
> - take the start and end coordinates of the domain in the protein
> (from the
> HMMer output)
> - convert those into the start and end coordinates of the domain in
> the cDNA
> - use the subseq() method to extract the subsequence of the cDNA
>
>
> Dave
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
More information about the Bioperl-l
mailing list