...

Exome Alignment Workflow

Bootstrapping-for-BQSR @ NCI-Gadi is a pipeline for bootstrapping a variant resource to enable GATK base quality score recalibration (BQSR) for non-model organisms that lack a publicly available variant resource. This implementation is optimised for the National Compute Infrastucture's Gadi HPC. Multiple rounds of bootstrapping can be performed. Users can use Fastq-to-bam @ NCI-Gadi and Germline-ShortV @ NCI-Gadi to ...

Local Cromwell implementation of GATK4 germline variant calling pipeline

See the GATK website for more information on this toolset

Assumptions

• Using hg38 human reference genome build
• Running 'locally' i.e. not using HPC/SLURM scheduling, or containers. This repo was specifically tested on Pawsey Nimbus 16 CPU, 64GB RAM virtual machine, primarily running in the /data volume storage partition.
• Starting from short-read Illumina paired-end fastq ...

Fastq-to-BAM @ NCI-Gadi is a genome alignment workflow that takes raw FASTQ files, aligns them to a reference genome and outputs analysis ready BAM files. This workflow is designed for the National Computational Infrastructure's (NCI) Gadi supercompter, leveraging multiple nodes on NCI Gadi to run all stages of the workflow in parallel, either massively parallel using the scatter-gather approach or parallel by sample. It consists of a number of stages and follows the BROAD Institute's best practice ...

SLURM HPC Cromwell implementation of GATK4 germline variant calling pipeline

See the GATK website for more information on this toolset

Assumptions

• Using hg38 human reference genome build
• Running using HPC/SLURM scheduling. This repo was specifically tested on Pawsey Zeus machine, primarily running in the /scratch partition.
• Starting from short-read Illumina paired-end fastq files as input

Dependencies

The following versions have been ...

A pipeline for mapping, calling, and annotation of SARS-CoV2 variants.

A workflow for mapping and consensus generation of SARS-CoV2 whole genome amplicon nanopore data implemented in the Nextflow framework. Reads are mapped to a reference genome using Minimap2 after trimming the amplicon primers with a fixed length at both ends of the amplicons using Cutadapt. The consensus is called using Pysam based on a majority read support threshold per position of the Minimap2 alignment and positions with less than 30x coverage are masked using ‘N’.

Workflow for NonSpliced RNAseq data with multiple aligners.

Steps:

• workflow_quality.cwl:
• FastQC (control)
• fastp (trimming)
• bowtie2 (read mapping)
• sam_to_sorted-bam
• featurecounts (transcript read counts)
• kallisto (transcript [pseudo]counts)

Workflow to build different indices for different tools from a genome and transcriptome.

This workflow expects an (annotated) genome in GBOL ttl format.

Steps:

• SAPP: rdf2gtf (genome fasta)
• SAPP: rdf2fasta (transcripts fasta)
• STAR index (Optional for Eukaryotic origin)
• bowtie2 index
• kallisto index

Alignment, assembly and annotation of generated transcripts from RNASEQ reads.

Alignment, assembly and annotation of RNASEQ reads as well as annotation of generated transcripts.

Alignment, assembly RNASEQ reads and annotation of generated transcripts.

Alignment, assembly and annotation of RNQSEQ reads using TOPHAT (without filtering out host reads).