Welcome! Here you can find a variety of tools and pipelines tailored to using PacBio HiFi Reads for metagenomics. In addition to the resources currently available, we will continue to add new tools as they are developed.
-
Genome-Binning-Pipeline
: Perform genome binning and MAG assessment on HiFi metagenomic assemblies with a streamlined workflow. Includes steps with minimap2, MetaBAT2, CheckM, and GTDB-Tk. Outputs high-quality MAG sequences and associated metadata. -
Taxonomic-Functional-Profiling-Protein
: Align HiFi reads to a protein database using DIAMOND and prepare inputs for MEGAN6, for the purpose of taxonomic and functional profiling. Provides access to taxonomic analyses (NCBI and GTDB), and allows functional analyses based on annotations from several databases (EC, SEED, InterPro2GO, eggNOG). -
Taxonomic-Profiling-Nucleotide
: Align HiFi reads to a nucleotide database using minimap2 and prepare inputs for MEGAN6, for the purpose of taxonomic profiling. This only provides access to NCBI taxonomic analysis, but may give higher resolution than with protein alignments. The GTDB taxonomy and functional annotations are not accessible downstream. -
MEGAN-RMA-Summary
: Obtain functional and taxonomic class count summaries for a set of RMA files created from the above MEGAN6 pipelines. Outputs absolute read counts and normalized read counts across samples for classes of each database, and summarizes the proportion of HiFi reads receiving functional and taxonomic annotations. Workflows are available for both protein-based RMA files and for nucleotide-based RMA files.
Each of these pipelines can be found in their respective folders. They are made available as Snakemake workflows. Snakemake is a python-based workflow manager. Snakemake workflows are highly portable because dependencies and environments are automatically setup using Anaconda/Conda. Snakemake also allows reproducibility, checkpointing, and the ability to scale workflows using HPC and cloud environments. Snakemake v5+ is required, and the workflows have been tested using v5.19+.
All available documentation can be found in the docs/
folder above.
Currently available documentation:
- Tutorial: Genome-Binning-Pipeline
- Tutorial: Taxonomic-Functional-Profiling-Protein
- Tutorial: Taxonomic-Profiling-Nucleotide
- Tutorial: MEGAN-RMA-Summary
THIS WEBSITE AND CONTENT AND ALL SITE-RELATED SERVICES, INCLUDING ANY DATA, ARE PROVIDED "AS IS," WITH ALL FAULTS, WITH NO REPRESENTATIONS OR WARRANTIES OF ANY KIND, EITHER EXPRESS OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, ANY WARRANTIES OF MERCHANTABILITY, SATISFACTORY QUALITY, NON-INFRINGEMENT OR FITNESS FOR A PARTICULAR PURPOSE. YOU ASSUME TOTAL RESPONSIBILITY AND RISK FOR YOUR USE OF THIS SITE, ALL SITE-RELATED SERVICES, AND ANY THIRD PARTY WEBSITES OR APPLICATIONS. NO ORAL OR WRITTEN INFORMATION OR ADVICE SHALL CREATE A WARRANTY OF ANY KIND. ANY REFERENCES TO SPECIFIC PRODUCTS OR SERVICES ON THE WEBSITES DO NOT CONSTITUTE OR IMPLY A RECOMMENDATION OR ENDORSEMENT BY PACIFIC BIOSCIENCES.