Skip to content
Branch: master
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Type Name Latest commit message Commit time
Failed to load latest commit information.


Various scripts to manipulate genomes and metagenomes on the SHARCNET HPC cluster
Including BLAST, BOWTIE, clustering (single linkage, MCL), aligning, treeing, etc.

These scripts are sometimes trivial, generally poorly written and often poorly documented
mea culpa

Written by Eric Collins based on an original set of scripts by Hugh Merz


A typical workflow might go something like....

1) Identify genomes of interest, e.g. 'Alteromonadales'
2) make blastlists and taxalists:
   ./ Alteromonadales AL
3) do all-vs-all blasts for Alteromonadales
   ./ blastlist-Alteromonadales
4) cluster genomes using single-linkage algorithm at several cutoffs of evalue/length fraction
   ./ AL_110001 taxalist-Alteromonadales-with-plasmids
5) make gene matrix from cluster list
   ./ taxalist-Alteromonadales-with-plasmids AL_110001.1e-10_0.7_long_on.cluster_list AL_110001
6) identify single-copy core genes
   R CMD BASH singlecopy.R
7) make concatenated alignments and ML tree from single-copy core genes
   ./ AL_110001.1e-10_0.7_long_on.single-copy.csv taxalist-Alteromonadales-with-plasmids AL_110001.1e-10_0.7_long_on.cluster_list AL_110001

You can’t perform that action at this time.
You signed in with another tab or window. Reload to refresh your session. You signed out in another tab or window. Reload to refresh your session.