The Vertebrate Genomes Pipelines in Galaxy are intended to allow a user to generate high-quality near error-free assemblies of species from a user's own data or from the GenomeArk database

IWC - Intergalactic Workflow Commission

Scaffolding using HiC data with YAHS.

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 ...

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 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)

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) ...

COVID-19: variation analysis on ARTIC PE data

The workflow for Illumina-sequenced ampliconic data builds on the RNASeq workflow for paired-end data using the same steps for mapping and variant calling, but adds extra logic for trimming amplicon primer sequences off reads with the ivar package. In addition, this workflow uses ivar also to identify amplicons affected by primer-binding site mutations and, if possible, excludes reads derived from such ...

Assembly with Hifi reads and Trio Data

Generate phased assembly based on PacBio Hifi Reads using parental Illumina data for phasing

Inputs

1. Hifi long reads [fastq]
2. Concatenated Illumina reads : Paternal [fastq]
3. Concatenated Illumina reads : Maternal [fastq]
4. K-mer database [meryldb]
5. Paternal hapmer database [meryldb]
6. Maternal hapmer database [meryldb]
7. Genome profile summary generated by Genomescope [txt]
8. Bloom Filter
9. Name of first haplotype
10. Name of second haplotype ...

dada2 amplicon analysis for paired end data

The workflow has three main outputs:

• the sequence table (output of makeSequenceTable)
• the taxonomy (output of assignTaxonomy)
• the counts which allow to track the number of sequences in the samples through the steps (output of sequence counts)

This workflow performs segmentation and counting of cell nuclei using fluorescence microscopy images. The segmentation step is performed using Otsu thresholding (Otsu, 1979). The workflow is based on the tutorial: https://training.galaxyproject.org/training-material/topics/imaging/tutorials/imaging-introduction/tutorial.html

Assemble long reads with Flye, then view assembly statistics and assembly graph

Run velocyto to get loom with counts of spliced and unspliced. It will extract the 'barcodes' from the bundled outputs.

This workflow processes the CMO fastqs with CITE-seq-Count and include the translation step required for cellPlex processing. In parallel it processes the Gene Expresion fastqs with STARsolo, filter cells with DropletUtils and reformat all outputs to be easily used by the function 'Read10X' from Seurat.

This workflow is composed with the XCMS tool R package (Smith, C.A. 2006) able to extract and the metaMS R package (Wehrens, R 2014) for the field of untargeted metabolomics.

https://training.galaxyproject.org/training-material/topics/metabolomics/tutorials/gcms/tutorial.html

MMGBSA simulation and calculation

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

...

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

This workflow is composed with the XCMS tool R package (Smith, C.A. 2006) able to extract, filter, align and fill gapand the possibility to annotate isotopes, adducts and fragments using the CAMERA R package (Kuhl, C 2012).

https://training.galaxyproject.org/training-material/topics/metabolomics/tutorials/lcms-preprocessing/tutorial.html

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 ...