Skip to main content
Join the official 2019 Python Developers SurveyStart the survey!

Parallel processing of sequencing reads

Project description

Pararead is a python package that simplifies parallel processing of DNA sequencing reads (BAM or SAM files), by parallelizing across chromosomes. Pararead is built for developers of python scripts that process data read-by-read. It enables you to quickly and easily parallelize your script.

Install

pararead is hosted on pypi. Install with:

pip install --user pararead

Or, within active environment:

pip install --upgrade pararead

Minimum working example

In the examples folder you can find examples/count_reads.py, which will count the number of reads in a sam/bam file in parallel.

Run this on your bam file like this:

python count_reads.py file.bam -O output.txt --cores 2

Look at the code to see how this is implemented.

Developing tools that use pararead

The main model provided is an abstract class calledParaReadProcessor, for which concrete children are created by implementing a __call__ method. This creates a callable instance that is then mapped over chromosomes.

The concept is generally described in this early blog post, which initiated the project that eventually became pararead. More details will be forthcoming.

Project details


Download files

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

Files for pararead, version 0.5.0
Filename, size File type Python version Upload date Hashes
Filename, size pararead-0.5.0.tar.gz (16.5 kB) File type Source Python version None Upload date Hashes View hashes

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN SignalFx SignalFx Supporter DigiCert DigiCert EV certificate StatusPage StatusPage Status page