Workflows

This workflow performs the scaffolding of a genome assembly using HiC data with YAHS. Can be used on any assembly with Hi-C data, and the assembly in the gfa format. You can generate a gfa from a fasta using the gfastat tool. Part of the VGP set of workflows, it is meant to be run after the contigging (workflows 3,4, or 5), optional purging step (Workflow 6 or 6b), and an optionnal scaffolding with Bionano data (Workflow 7). This workflow includes QC with Assembly statistics, Busco, and Hi-C maps. ...

This workflow constructs Metagenome-Assembled Genomes (MAGs) using SPAdes or MEGAHIT as assemblers, followed by binning with four different tools and refinement using Binette. The resulting MAGs are dereplicated across the entire input sample set, then annotated and evaluated for quality. You can provide pooled reads (for co-assembly/binning), individual read sets, or a combination of both. The input samples must consist of the original reads, which are used for abundance estimation. In all cases, ...

Microbiome - Variant calling and Consensus Building

This workflow performs subtyping and consensus sequence generation for batches of Illumina PE sequenced Influenza A isolates.

Workflow for variant analysis against a reference genome in GenBank format

Find and annotate variants in ampliconic SARS-CoV-2 Illumina sequencing data and classify samples with pangolin and nextclade

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

Microbiome - QC and Contamination Filtering

This workflow takes as input a collection of paired fastqs. Remove adapters with cutadapt, map pairs with bowtie2. Keep MAPQ30 and concordant pairs. MACS2 for paired bam.

This workflow takes as input a collection of fastqs (single reads). Remove adapters with cutadapt, map with bowtie2. Keep MAPQ30. MACS2 for bam with fixed extension or model.