Microbiome - QC and Contamination Filtering

Pathogens of all samples report generation and visualization

Annotation of an assembled bacterial genomes to detect genes, potential plasmids, integrons and Insertion sequence (IS) elements.

This workflow

• Reconstruct phylogeny (insert fragments in a reference)
• Alpha rarefaction analysis
• Taxonomic analysis

Create Meryl Database used for the estimation of assembly parameters and quality control with Merqury. Part of the VGP pipeline.

VGP Workflow #1

This workflow produces a Meryl database and Genomescope outputs that will be used to determine parameters for following workflows, and assess the quality of genome assemblies. Specifically, it provides information about the genomic complexity, such as the genome size and levels of heterozygosity and repeat content, as well about the data quality.

Inputs

• A collection of Hifi long reads in FASTQ format
• k-mer length
• Ploidy

Outputs

• Meryl Database of kmer counts

...

Contiging Solo w/HiC:

Generate phased assembly based on PacBio Hifi Reads using HiC data from the same individual for phasing.

Inputs

1. Hifi long reads [fastq]
2. HiC forward reads (if multiple input files, concatenated in same order as reverse reads) [fastq]
3. HiC reverse reads (if multiple input files, concatenated in same order as forward reads) [fastq]
4. K-mer database [meryldb]
5. Genome profile summary generated by Genomescope [txt]
6. Name of first assembly
7. Name of second assembly ...

Purge contigs marked as duplicates by purge_dups (could be haplotypic duplication or overlap duplication). This workflow is the 6th workflow of the VGP pipeline. It is meant to be run after one of the contigging steps (Workflow 3, 4, or 5)

Scaffolding with Bionano

Scaffolding using Bionano optical map data

Inputs

1. Bionano data [cmap]
2. Estimated genome size [txt]
3. Phased assembly generated by Hifiasm [gfa1]

Outputs

1. Scaffolds
2. Non-scaffolded contigs
3. QC: Assembly statistics
4. QC: Nx plot
5. QC: Size plot

Purge Duplicate Contigs

Purge contigs marked as duplicates by purge_dups in a single haplotype(could be haplotypic duplication or overlap duplication) This workflow is the 6th workflow of the VGP pipeline. It is meant to be run after one of the contigging steps (Workflow 3, 4, or 5)

Inputs

1. Genomescope model parameters [txt] (Generated by the k-mer profiling workflow)
2. Hifi long reads - trimmed [fastq] (Generated by Cutadapt in the contigging workflow)
3. Assembly to purge (e.g. hap1) ...