An SQL-based solution for large-scale genomic analysis
Project description
pysequila
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, 3.8
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
- support for both SQL and Dataframe/Dataset API
Setup
$ python -m pip install --user pysequila or (venv)$ python -m pip install pysequila
Usage
$ python >>> from pysequila import SequilaSession >>> ss = SequilaSession \ .builder \ .config("spark.jars.packages", "org.biodatageeks:sequila_2.12:0.7.5") \ .config("spark.driver.memory", "2g") \ .getOrCreate() >>> ss.sql( f""" CREATE TABLE IF NOT EXISTS reads USING org.biodatageeks.sequila.datasources.BAM.BAMDataSource OPTIONS(path "/features/data/NA12878.multichrom.md.bam") """ >>> ss.sql ("SELECT * FROM coverage('reads', 'NA12878','/features/data/Homo_sapiens_assembly18_chr1_chrM.small.fasta") >>> # or using DataFrame/DataSet API >>> ss.coverage("/features/data/NA12878.multichrom.md.bam", "/features/data/Homo_sapiens_assembly18_chr1_chrM.small.fasta")
ChangeLog
0.1.0 (2020-09-16)
- Initial release.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
pysequila-0.3.3.tar.gz
(18.8 kB
view hashes)
Built Distribution
Close
Hashes for pysequila-0.3.3-py2.py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 0f812b99a0bc8ada476c4dab5fd965a80e65f2e525708dd80898b67061562174 |
|
MD5 | df2f5f000f6dd56c2cd8c77ccc50efca |
|
BLAKE2-256 | c131072f7d37403b1a8ba7b49bbbbc4e7c489944e57336bc32d2e5f27d6415f4 |