Giter Club home page Giter Club logo

eyre's Introduction

Eyre

It is a similar workflow as Nullarbor, but it only runs isolate specific analysis. No phylogenetic tree and pangenome analysis will be generated. A snakemake workflow that process fastq files of bacterial pathogens to produce:

  1. Sequencing yield (fq),
  2. Species identification (Kraken2),
  3. de novo assemblies (shovill) and their qc (fa)
  4. Subtyping (MLST)
  5. Antimicrobial resistance profile (abricate)

Requirements

It is practically a snakemake workflow with all other softwares from Nullarbor.

  1. Nullarbor

It is best to install Nullarbor using conda.

$ conda create -n eyre nullarbor
  1. Snakemake

The snakemake software should be added to the Nullarbor conda environment

$ conda activate eyre
$ conda install snakemake
  1. Git clone the eyre repository
$ git clone https://github.com/lexleong/eyre.git
  1. Download Kraken2 database

Modification on the Snakefile is required to direct the pipeline script to the directory path containing kraken2 database files (hash.k2d, taxo.k2d, and opts.k2d).

Workflow

  1. Modify the sequencing submission sheet as per the SampleSheet.csv template

  2. Run the bcl2eyre.sh script

$ bcl2eyre.sh [dir]/SampleSheet.csv

Etymology

The Nullarbor is a huge treeless plain that spans the area between South Australia and Western Australia. If one were to travel from Adelaide to Nullarbor, one will have to go through Eyre Peninsula. So a South Australian public health microbiologist has to perform the Eyre pipeline prior to the Nullarbor pipeline for their bacterial WGS.

Disclaimer

This workflow is specific to working in slurm and conda environments. It is used mainly by SA Pathology MID PHL for downstream processing of bacterial whole genome sequencing output.

eyre's People

Contributors

lexleong avatar

Stargazers

 avatar

Watchers

 avatar

eyre's Issues

Removing short sequences

You may want to use the more portable seqtk seq -L 1000 in.fna > out.fna for this.

rule remove_small:
   input:
        "{sample}/shovill"
   output:
        "filtered_contigs/{sample}.fna"
   shell:
        "perl /data/sacgf/microbio/lleong/commands/removesmalls.pl 1000 {input}/contigs.fa > {output}"

2021-05-27

Testing issue and what the hell is this about.... :(

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.