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

List coverage blind spots

Project description

This software lists all positions in given bam files, that have no or only very limited coverage and stores the result in a hdf5 file.

Usage

cbs [-h] [--min-cov MIN_COV] files [files ...] output

If multiple files are given, a position is designated a blind spot, if no file contains a coverage of at least MIN_COV at this position.

Output

The output is a hdf5 file with the following structure. For every chromosome a group is created, that contains a dataset called missing_cov. The dataset contains rows of tuples of 32bit integer, representing the first blind spot and the first position having enough coverage.

Install

This software requires python3, samtools, hdf5, numpy and h5py. All python dependencies can be installed by installing cbs via pip:

pip install cbs

Project details


Release history Release notifications

Download files

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

Files for cbs, version 1.0.0
Filename, size File type Python version Upload date Hashes
Filename, size cbs-1.0.0.tar.gz (4.3 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