Skip to main content

An SQL-based solution for large-scale genomic analysis

Project description

pysequila is Python entrypoint to SeQuiLa, an ANSI-SQL compliant solution for efficient sequencing reads processing and genomic intervals querying built on top of Apache Spark. Range joins, depth of coverage and pileup computations are bread and butter for NGS analysis but the high volume of data make them execute very slowly or even failing to compute.

Requirements

  • Python 3.7

Features

  • custom data sources for bioinformatics file formats (BAM, CRAM, VCF)
  • depth of coverage calculations
  • pileup calculations
  • reads filtering
  • efficient range joins
  • other utility functions

Setup

$ python -m pip install --user pysequila
or
(venv)$ python -m pip install pysequila

Usage

$ python
>>> from pyspark.sql import SparkSession
>>> spark = SparkSession \
    .builder \
    .appName(f'{app_name}') \
    .getOrCreate()
>>> from sequila import SequilaSession
>>> ss = SequilaSession(spark)
>>> ss.sql ("SELECT * FROM  coverage('reads', 'NA12878'")
>>>

ChangeLog

0.1.0 (2020-09-16)

  • Initial release.

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 pysequila, version 0.1.3
Filename, size File type Python version Upload date Hashes
Filename, size pysequila-0.1.3-py2.py3-none-any.whl (3.1 kB) File type Wheel Python version py2.py3 Upload date Hashes View
Filename, size pysequila-0.1.3.tar.gz (4.9 kB) File type Source Python version None Upload date Hashes View

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring Facebook / Instagram Facebook / Instagram PSF Sponsor Fastly Fastly CDN Google Google Object Storage and Download Analytics Huawei Huawei PSF Sponsor Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Salesforce Salesforce PSF Sponsor Sentry Sentry Error logging StatusPage StatusPage Status page