Skip to content

Peter-J-Freeman/NextFlow_RD_Genomic

Folders and files

NameName
Last commit message
Last commit date
Feb 28, 2024
Dec 14, 2023
Mar 6, 2024
Dec 14, 2023
Feb 14, 2024
Sep 28, 2023
Jan 18, 2024
Feb 16, 2024
May 9, 2024
Feb 14, 2024
Dec 12, 2023
Feb 28, 2024
Feb 1, 2024

Repository files navigation

NextFlow_RD_Genomic

Description

A simple base Rare disease and germline genomics pipeline to test the effects of down-sampling on variant calling

Basic Overview

Using the NextFlow workflow software to run the following pipeline

Pipeline

Index genome > Fastqc analysis > Align reads > Downsample bam files > Sort bam > Mark duplicates > Index bam > Call variants > Hard filter

Setup

To run the pipeline, we need to obtain

  • A genome build (GRCh38) - provided by the Broad institute
$ cd data/genome
$ wget https://storage.googleapis.com/gcp-public-data--broad-references/hg38/v0/Homo_sapiens_assembly38.fasta
  • FastQ sample (for workflow development)
$ cd ../samples
$ wget https://ftp.sra.ebi.ac.uk/vol1/fastq/SRR151/003/SRR1518253/SRR1518253_1.fastq.gz && \
wget https://ftp.sra.ebi.ac.uk/vol1/fastq/SRR151/003/SRR1518253/SRR1518253_2.fastq.gz && \
gunzip *.gz
  • When scaling up, FastQ samples
$ wget https://genomics.viapath.co.uk/benchmark/files/FASTQ/NA12878_WES.zip

Running the pipeline

# Using Docker
$ nextflow run -profile docker main.nf

# Using docker in singularity
$ nextflow run -profile singularity main.nf

Validating the pipeline

See https://genomics.viapath.co.uk/benchmark

DNANexus applet setup (A local applet for basic testing)

pip install -r requirements.txt
  • Routine maintenance Periodically update dxpy
$ pip install --upgrade dxpy

DNANexus Tutorial

$ dx select <your-project-name>
$ dx build --nextflow

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published