Greetings - I am studying the evolution of several families of g protein coupled receptors and I would like to identify positively selected sites and determine if lineages (subtrees) have evidence for positive selection.  One of the families, in particular, is quite diverse and parts of the alignment are 'questionable.'  I was eager to try out datamonkey so I plugged in the alignment and ran the SLAC analysis.  Out of 1701 sites, this analysis identified 419 positively selected sites and 45 under negative selection.  However, the analysis results also report:
Quote:


There are ~100 sequences in the file.  Well, I know the alignment still needs work, but this brings up the question of how close the sequences need to be in order to get meaningful results from this and other analyses.  I can remove distant and unalignable sequences from the file, and divide it into clusters of sequences that align well within the cluster but not between clusters. Is there some optimum or maximum substitutions per site that I should try to achieve?
Any comments or suggestions are welcome.
Mike

Hi Sergei -

I realized that deleting columns from the alignment was going to be difficult to do quickly without disrupting the reading frame so instead I deleted about 10 sequences that were causing the most gaps.  These sequences appear to have large insertions that may be incorrectly annotated exons. After removing these sequences I have 93 sequences and 734 codons. the SLAC analysis using the Resolve ambiguities option results in 619 positively selected sites.

Thanks
Mike