Protein rate shift analysis server

About this service

This is a server that provides rate shift analyses for proteins. It takes a number of aligned proteins as input. The two families to be analysed should be inserted in the first and second text areas. The outgroup sequences should be inserted in the last text area. All the sequences should be part of one big alignment, and should be inserted in fasta format. If a phylogeny is input, it should be in the newick format. The algorithm is described in the reference below.

Comments are very much appreciated.

At the moment, the server has a maximum of 45 sequences in each family and 10 sequences in the outgroup. The maximum alignment length is 600.

Bjarne Knudsen (bk@daimi.au.dk)

Sequence form

When you press `Calculate' your data will be processed.

E-mail address:

Family 1:


Family 2:


Optional outgroup:


Protein rate matrix:
     Jones Taylor Thornton (recommended)
     Jones Taylor Thornton: transmembrane
     Jukes Cantor

Tree:
    Optional tree: (Neighbour joining is used if not specified)
    

     Maximum likelihood optimize branchlengths

If an outgroup is not specified, the root is used as the split point between the two families, if it is binary, does indeed split the tree into the appropriate sub-families and branch length optimization is not chosen. Otherwise, the root of the tree is put in the center of the branch connecting the two families.

Reference

  • Knudsen, B. and Miyamoto, M. M. (2001) A likelihood ratio test for evolutionary rate shifts and functional divergence among proteins. Proc. Natl. Acad. Sci. USA, 98(25), 14512-14517.