## ----setup, include = FALSE----------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----message = FALSE-----------------------------------------------------
library(metagenomeFeatures)
gg85 <- get_gg13.8_85MgDb()

## ------------------------------------------------------------------------
gg85

## ------------------------------------------------------------------------
data_dir <- system.file("extdata", package = "metagenomeFeatures")
soil_gg_ids <-  readRDS(file.path(data_dir, "qiita_study_94_gg_ids.RDS"))

## ------------------------------------------------------------------------
soil_mgF <- annotateFeatures(gg85, soil_gg_ids)

## ------------------------------------------------------------------------
soil_mgF

## ------------------------------------------------------------------------
mgF_seq(soil_mgF)

## ------------------------------------------------------------------------
mgF_tree(soil_mgF)