Search this site
Embedded Files
Metagenomics
  • Metagenomics
    • Taxonomy
      • Alpha and beta diversity
      • Strain level
      • Pangenome
      • Marker genes
        • MLST
      • 16S
        • Operational taxonomic unit (OTU)
        • 16s vs shotgun
    • Phylogeny
      • Horizontal gene transfer (HGT)
      • Lowest Common Ancestor (LCA)
      • Long branch attraction (LBA)
    • Evolution
    • Ecosystem
      • Environmental gene tags (EGTs)
      • Soil metagenomics
    • Metatranscriptomics
      • RNA-seq versus microarray
    • QC
      • Coverage depth
      • Sampling depth
        • Calculate sampling depth
      • RPKM calculation
      • Phred score (Q score)
      • Giga base pairs
      • GC content
    • Wiki
      • Microbiome
      • Orthologs and paralogs
      • Short read mapping
      • Fecal microbiota transplant
      • DNA sequencing library
      • Multiplex sequencing
      • x
      • x
      • x
      • x
        • x
      • x
      • x
      • x
  • Tools
    • 16S tools
      • DADA2
        • Conda environment
          • Segmentation fault
      • Normalization
        • CSS
      • QIIME
        • Alpha & Beta diversity
        • Install QIIME 1
          • AttributeError:axisbg
          • biom error
          • Greengenes
          • h5py BIOM error
          • PackagesNotFoundError
          • qiime_config
          • R for QIIME
          • usearch61
        • OTU biom table
          • QIIME: split OTU table
          • ValueError column index exceeds matrix dimensions
        • QIIME mapfile
        • QIIME OTU clustering
          • Error running usearch61
          • Fungi 18S
          • IOError Errno 28 No space left on device
          • SILVA
        • QIIME pre-processing
          • Cannot find fastq-join
          • join-error
          • split-error
        • Taxonomy
        • Test statistics
    • Assembly
      • Contigs
      • N50 statistics
      • MEGAHIT
    • BLAST
      • BLAST error: Too many positional arguments
      • BLAST word-size
      • BLASTn output format 6
      • BLASTx
      • E-value & Bit-score
      • FastANI
        • Install FastANI
          • Install autoconf
          • Install GLS
      • Generate_database
        • Public sequences
      • Install BLAST
      • megablast
    • Bowtie2
      • Create bowtie2 index
      • Install bowtie2
        • libtbb.so.2
        • Segmentation fault
    • Genome
      • Prokka
    • Pathogen screening
    • Phylogenetic tree
      • File format
      • Sequence-Alignment
        • Alignment viewer
      • Tree-construction
        • FastTree
        • RAxML
          • Add bipartition
          • Install RAxML
      • Tree-viewer
        • Forester
          • Java error
    • SAMtools
      • BCFtools
      • Consensus sequence
      • Converting BAM to fastq
      • Error
        • Could not parse the header line
        • different line length in sequence
      • Install
        • Error curses.h
      • Number of reads in bam file
      • SAM file format
      • SAMtools: get breadth of coverage
    • Sequence data
      • Convert fastq to fasta
      • Extract sequence subset
      • Get random subset
      • Multi-FASTA format
      • NCBI ftp genome download
        • gff to ffn
    • Shotgun sequencing
      • Alignment
        • Viewer
          • Tablet
      • Data
      • NCBI SRA files
        • Install SRA-tools
        • prefetch
        • Error
          • path not found while resolving tree within virtual file system module
        • wget download
        • old fastq-dump
      • Quality control
        • Trim Galore
      • Remove host sequences
      • Remove too short reads from fastq files
      • Fastq file format
    • Ubuntu Linux
      • awk
      • bzip2
      • Extract columns from file
      • File properties
      • gzip & tar.gz
      • Loop over list of files
        • find
        • if greater than
        • String split
      • Rename multiple files
      • rsync
      • sed
      • Shell stderr redirect
      • Show disk space usage
      • Split large file
      • Data tables
        • Transpose table
        • Sorting
    • HPC SGE cluster
      • Submit job
      • Check cluster nodes
Metagenomics

RAxML

Phylogenetic tree


Create a phylogenetic tree from multiple sequence alignment

→ Install RAxML

# Documentation

http://sco.h-its.org/exelixis/web/software/raxml/hands_on.html

# Create tree from → multiple alignment sequences

raxmlHPC -m GTRGAMMA -p 12345 -s alignment.aln -n NAME

# Run RAxML in bootstrap mode

raxmlHPC -m GTRGAMMA -p 12345 -s alignment.aln -n NAME_bootstrap -f a -x 12345 -N 100 -T 12

# Results

RAxML_bestTree.NAME_bootstrap best-scoring ML tree

RAxML_bipartitions.NAME_bootstrap best-scoring ML tree with support values

RAxML_bipartitionsBranchLabels.NAME_bootstrap best-scoring ML tree with support values as branch labels

RAxML_bootstrap.NAME_bootstrap all bootstrapped trees

RAxML_info.NAME_bootstrap program info

# view tree

open RAxML_bestTree.NAME_bootstrap with → Forester


raxmlHPC -help

-p random seed

-s filename of aligned sequences

-n output filename

-T number of threads/CPUs

-f a rapid Bootstrap analysis and search for best-scoring ML tree in one program run

-m GTRGAMMA : GTR + Optimization of substitution rates + GAMMA model of rate

heterogeneity (alpha parameter will be estimated).

With the optional "X" appendix you can specify a ML estimate of base frequencies.

-N 100 Specify the number of alternative runs on distinct starting trees

In combination with the "-b" option, this will invoke a multiple boostrap analysis

-w FULL (!) path to the directory into which RAxML shall write its output files

-x Specify an integer number (random seed) and turn on rapid bootstrapping

CAUTION: unlike in version 7.0.4 RAxML will conduct rapid BS replicates under

the model of rate heterogeneity you specified via "-m" and not by default under CAT


read more

→ RAxML Manual.pdf

→ RAxML Website

→ RAxML Tutorial

www.metagenomics.wiki

Author: Matthias Scholz

Google Sites
Report abuse
Page details
Page updated
Google Sites
Report abuse