Build a core-genome phylogeny

From www.wiki.chlamdb.ch
Revision as of 12:07, 24 October 2018 by Tpillone (talk | contribs)
Jump to navigation Jump to search


Overview

This tutorial describes step by step the reconstruction of a core genome phylogeny for a given clade of bacteria.

Setup

  • clone eutils repository
1# download repo from metagenlab
2git clone https://github.com/metagenlab/eutils.git
3# move to eutils directory
4cd eutils
5# add path to pythonpath and path
6echo "export PATH=:$PWD"':$PATH' >> ~/.bashrc
7echo "export PYTHONPATH=:$PWD"':$PYTHONPATH' >> ~/.bashrc
8# reload bashrc
9source ~/.bashrc
  • install miniconda
  • create a dedicated environment and install orthofinder
1# create a new environment
2conda create -n orthofinder
3source activate orthofinder
4# install orthofinder
5conda install orthofinder
6# install biopython
7conda install biopython


Get table of genome assemblies

1df -h

Download gbff files

1df -h

Prepare fasta files for orthofinder

Download gbff files