Data are converted to phyloseq format.

phyloseq()
# phyloseq object for experiment samples
metadata %>%
  filter(SourceSink == "sink")

# phyloseq object to compare calculus samples
LS0tCnRpdGxlOiAiRGlmZmVyZW50aWFsIGFidW5kYW5jZSIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CiN8IGxhYmVsOiBzZXR1cAojfCBpbmNsdWRlOiBmYWxzZQpsaWJyYXJ5KHBoeWxvc2VxKQpsaWJyYXJ5KEFOQ09NQkMpCm90dV90YWJsZSA8LSByZWFkcjo6cmVhZF9jc3YoKQpgYGAKCkRhdGEgYXJlIGNvbnZlcnRlZCB0byBgcGh5bG9zZXFgIGZvcm1hdC4KCmBgYHtyfQpwaHlsb3NlcSgpCmBgYAoKYGBge3J9CiMgcGh5bG9zZXEgb2JqZWN0IGZvciBleHBlcmltZW50IHNhbXBsZXMKbWV0YWRhdGEgJT4lCiAgZmlsdGVyKFNvdXJjZVNpbmsgPT0gInNpbmsiKQoKIyBwaHlsb3NlcSBvYmplY3QgdG8gY29tcGFyZSBjYWxjdWx1cyBzYW1wbGVzCmBgYAoK