Author / Distributor: Alexandros Stamatakis, The Exelixis Lab1,
Ludwig-Maximilians-Universit ¨at M¨unchen
If you use RAxML please always cite the following paper: Alexandros Stamatakis : “RAxML-VI-HPC:
Maximum Likelihood-based Phylogenetic Analyses with Thousands of Taxa and Mixed Models”, Bioinformatics
22(21):2688–2690, 2006 [4].
RAxML (Randomized Axelerated Maximum Likelihood) is a program for sequential and parallel Maximum
Likelihood based inference of large phylogenetic trees.
Frist one, it will automatically determine the best-scoring AA substitution model on a fixed starting tree.
For unpartitioned datasets execute it as following. The outfile will then contain the best-scoring AA model to use with RAxML.
The second perl script will convert FASTA files to RAxML’s relaxed PHYLIP format, i.e., not truncate sequence names. It will generate a file alignment.fasta.phylip that can be read by RAxML.
Installation: source code downloaded from RAxML website.
Pthreads-parallelized version is compiled with the gcc compiler by typing make -f Makefile.MPI.
Version 7.0.4 is the last version, located at /usr/local/RAxML/RAxML-7.0.4
Shortcut /usr/local/RAxML-latest is linked with Version 7.0.4 now .
Please refer to pdf file for detail commands
/usr/local/RAxML-latest raxmlHPC-MPI -# numberOfRuns -s sequenceFileName -n outputFileName -m substitutionModel -b bootstrapRandomNumberSeed
The script submitted to the queue, also refer to submit jobs to queues at rcluster,IOB#!/bin/bash
Then submit the job as:cd working_directory
echo $LSB_HOSTS
cat /dev/null > mlist.$$
for variable in $LSB_HOSTS; do
echo $variable >> mlist.$$
done
mpirun -np 2 -machinefile mlist.$$ /usr/local/RAxML-latest raxmlHPC-MPI < -# numberOfRuns -s sequenceFileName -n outputFileName -m substitutionModel -b bootstrapRandomNumberSeed [other options]> myoutput
rm -f mlist.$$
bsub -n 2 -q queuename -o test.%J.out -e test.%J.err ./subp.sh
There are two more perl scripts toolsFrist one, it will automatically determine the best-scoring AA substitution model on a fixed starting tree.
For unpartitioned datasets execute it as following. The outfile will then contain the best-scoring AA model to use with RAxML.
perl /usr/local/RAxML-latest/ProteinModelSelection.pl alignmentFile.phylip > outfile
For partitioned datasets execute it as following. The outfile will then contain the best-scoring AA model for every partition.
perl /usr/local/RAxML-latest/ProteinModelSelection.pl alignmentFile.phylip partitionData.txt > outfile
The second perl script will convert FASTA files to RAxML’s relaxed PHYLIP format, i.e., not truncate sequence names. It will generate a file alignment.fasta.phylip that can be read by RAxML.perl /usr/local/RAxML-latest/fasta2relaxedPhylip.pl alignment.fasta