This report documents the binning of a Nitrospira genome from the GWW sample. See Daims et al., 2015: Complete Nitrification by Nitrospira Bacteria for further details.

Load the mmgenome package

The metagenome data is analysed using the mmgenome package.


Import data

The Rmarkdown file Load_data.Rmd describes the loading of the data and can be imported using the mmimport function. However, the preprocessed data can also be downloaded directly from figshare: Daims_GWW. Hence, here we import the prepocessed data from figshare instead.


Extract Nitrospira 7

As seen from the overview plots, there are numerous Nitrospira species. This one is located between multiple other species.

p <- mmplot(data = d,
            x = "HPD", 
            y = "HPF1", 
            log.x = T, 
            log.y = T, 
            color = "essential", 
            minlength = 1000,
            factor.shape = "solid") +
  scale_x_log10(limits=c(20, 120)) +
  scale_y_log10(limits=c(20, 100))

#sel <- mmplot_locator(p)

sel <- data.frame(HPD  =  c(52.6, 61.4, 44.5, 30.9, 31, 42),
                  HPF1  =  c(38.3, 31.6, 21.7, 22.2, 27, 36.2))

mmplot_selection(p, sel)