Image

Image How to use PASSML

PARAMETER SETTING MENU

*********************************************************

PaSSml --> Phylogeny and Secondary Structure using ML

*********************************************************

this program infers branch lengths of an evolutionary tree from a multialignment of amino acid sequences and predicts or improves secondary structure

- job file name: passml.inp (ASCII format)

- verbose file name: prolix (all the trees evaluated)

PARAMETERS:

a) secondary structure y/n:

b) accessibility y/n:

c) number sequences :

d) protein length :

e) user tree y/n:

f) use branch lengths y/n:

g) sequence filename (a short one):

h) sequential or aligned s/i:

i) output file (a short one):

j) eigen file: eigen

k) model file: rhoij38

l) aa stationary probabilities file: rev38

m) HMM (a), HMM &1st category is Coil exposed (b) [NULL model (c)]:

n) structure prediction: No (a), H(helix)_E(sheet)_O(others) (b), H_E_T(turn)_O (c), b(buried)_e(exposed) (d), H_E_T_O * b_e (e):

o) [ancestral node prediction]: n

p) [bootstrap analysis] y/n: n

q) quit y/n: n

choose a letter (. to accept) :>

INPUT FILE

  • The file should contain the set of sequences in a sequential or interleaved format.

  • each sequence has a 10 character name (spaces are allowed).

  • It is basically the same format accepted by the Phylip programs without the number of

    species and sites at the beginning.

  • Gaps and unknown amino acids are represented by the "-" character.

  • A tree with or without branch lengths can follow an empty line after the set of sequences.

    OUTPUT FILES

  • The user-named output file reports the chosen parameters, the tree topology (the branch lengths and links) and the maximum likelihood value; then the 'best' tree in a standard format, e.g.

    (species_A:8.084913,(species_B:0.904446,species_C:0.880132):23.084233,species_D:0.263436);

    Then, for all the sites, the secondary structure predictions in terms of posterior probabilities that each site is in each of the chosen categories, e.g.

    Secondary Structure prediction

    Number of predicted categories: 4

    siteAlpha helixBeta SheetTurnCoil+others
    10.7430700.1745810.0128320.069516
    20.7761100.1845170.0133440.026029
    30.7904180.1829280.0143390.012316
    ....
    ....
    n0.008849 0.1222150.5370830.331853

  • The program generates the "prolix" file, reporting all the likelihood values and branch lengths for all the trees evaluated.

    Image