Skip to main content

A pipeline to deplete reads from a given reference

Project description

https://badge.fury.io/py/sequana-depletion.svg JOSS (journal of open source software) DOI https://github.com/sequana/depletion/actions/workflows/main.yml/badge.svg

This is is the depletion pipeline from the Sequana project

Overview:

select or deplete reads from input FastQ files given a reference

Input:

Fastq Files

Output:

Fastq Files

Status:

production

Citation:

Cokelaer et al, (2017), ‘Sequana’: a Set of Snakemake NGS pipelines, Journal of Open Source Software, 2(16), 352, JOSS DOI doi:10.21105/joss.00352

Installation

Install this package as follows:

pip install sequana_depletion

it requires https://sequana.readthedocs.io and https://bioconvert.readthedocs.io for the bam to fastq conversion

Usage

sequana_depletion --help
sequana_depletion --input-directory DATAPATH --reference hg38.fa --mode depletion
sequana_depletion --input-directory DATAPATH --reference covid.fa --mode selection

This creates a directory with the pipeline and configuration file. You will then need to execute the pipeline:

cd depletion
sh depletion.sh  # for a local run

This launch a snakemake pipeline. If you are familiar with snakemake, you can retrieve the pipeline itself and its configuration files and then execute the pipeline yourself with specific parameters:

snakemake -s depletion.rules -c config.yaml --cores 4 --stats stats.txt

Or use sequanix interface.

Requirements

This pipelines requires the following executable(s):

  • bwa

  • samtools

  • bamtools

#.. image:: https://raw.githubusercontent.com/sequana/depletion/master/sequana_pipelines/depletion/dag.png

Details

This pipeline runs depletion in parallel on the input fastq files (paired or not). A brief sequana summary report is also produced.

Rules and configuration details

Here is the latest documented configuration file to be used with the pipeline. Each rule used in the pipeline may have a section in the configuration file.

Changelog

Version

Description

0.3.0

  • use click and new sequana_pipetools. convert to pyproject

0.2.0

  • handle paired/unpaired data

  • refactorise to use containers/apptainers

0.1.0

First release.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

sequana_depletion-0.3.0.tar.gz (5.9 kB view details)

Uploaded Source

Built Distribution

sequana_depletion-0.3.0-py3-none-any.whl (8.0 kB view details)

Uploaded Python 3

File details

Details for the file sequana_depletion-0.3.0.tar.gz.

File metadata

  • Download URL: sequana_depletion-0.3.0.tar.gz
  • Upload date:
  • Size: 5.9 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.11.7

File hashes

Hashes for sequana_depletion-0.3.0.tar.gz
Algorithm Hash digest
SHA256 8c5cd1a2972e7517319a993731ada01ffd595e160b3c4fdfac74f8699a246232
MD5 a67cad72a22ec0bab5bdef2c57d3106b
BLAKE2b-256 f6398c30ca4d2fd8e69b942f8151b1bdb67aeafd6dfa208c7f1a1e5f6c903287

See more details on using hashes here.

File details

Details for the file sequana_depletion-0.3.0-py3-none-any.whl.

File metadata

File hashes

Hashes for sequana_depletion-0.3.0-py3-none-any.whl
Algorithm Hash digest
SHA256 474373de3fcec58676b4f42eb68c081a4bc2383a60d2a67d9a5eca61cd2d7243
MD5 bbf8128e0533d05b1b11336f116027d0
BLAKE2b-256 63648c0900ae60782194bc913ce756a87a881e8db74fa68a59395e9c63c86fa4

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page