Michael Gerth
  • News & Blog
  • About
  • Research
    • Experimental evolution of Spiroplasma after host shifts
    • Wolbachia phylogenomics
    • Rickettsial symbionts in green lacewings
    • The role of Wolbachia in quill mites
    • Methods in molecular phylogenetics
  • Publications
  • Resources

Two new scripts for creating simple GC-coverage plots from SPAdes assemblies and analysing PhyloBayes trace files

18/12/2017

0 Comments

 
Before writing about science again (new post is in the making), I have uploaded two scripts that I find useful for my work:
​

gc_cov.pl 

This script is basically a very simplified version of the blobplots function from the BlobTools package. It creates GC-coverage plots directly from SPAdes assembly files, without the need for mapping the reads back to the assembly. Since I use mainly SPAdes anyway, this has been quite handy. The script will also annotate the plot when a taxonomy file is provided, which can be generated, e.g., from blast outputs. Below is an example for the plots that can be generated with the script.
Picture

pbplot.R

This is a tool to plot trace parameters estimated during PhyloBayes analyses, directly from the directory in which the analyses are/were run. I use this to check if chains appear to be converging, to identify "problematic" chains, and also to estimate suitable burnin sizes. It is not meant to provide an exhaustive diagnosis of the mcmc runs, but rather as a quick look into ongoing analyses – the plots are generated almost instantly and plotted to the screen or pdf. See below for examples.
​
Picture
Trace parameter vs iteration plots.
Picture
Trace parameter density plots.
​
Find both scripts on the resources page or on github. I'd be very happy about any feedback and comments!
0 Comments



Leave a Reply.

    Welcome!

    This is the website of Michael Gerth. I am a biologist with an interest in insects and the microbes within them. Click here to learn more.

    Archives

    October 2019
    March 2019
    January 2019
    February 2018
    December 2017
    July 2017
    December 2016
    April 2016
    March 2016

    Categories

    All
    Announcements
    Bees
    General
    Gene Repertoire Analysis
    Howto
    Long Branch Attraction
    Phylogenetics
    Post-publication Review
    Publications
    Scripts
    Wolbachia

    RSS Feed


    Tweets by gerth_micha
  • News & Blog
  • About
  • Research
    • Experimental evolution of Spiroplasma after host shifts
    • Wolbachia phylogenomics
    • Rickettsial symbionts in green lacewings
    • The role of Wolbachia in quill mites
    • Methods in molecular phylogenetics
  • Publications
  • Resources