Predict splicing variant effect from VCF
Project description
# mmsplice
[![pypi](https://img.shields.io/pypi/v/mmsplice.svg)](https://pypi-hypernode.com/pypi/mmsplice)
[![travis](https://img.shields.io/travis/s6juncheng/mmsplice.svg)](https://travis-ci.org/s6juncheng/mmsplice)
Predict splicing variant effect from VCF
* Free software: MIT license
## Usage example
------
NOTE: make sure you split and left-normalize the input VCF file.
Check notebooks/example.ipynb
```python
# Import
from mmsplice.vcf_dataloader import SplicingVCFDataloader
from mmsplice import MMSplice, predict_all_table
from mmsplice.utils import max_varEff
# example files
gtf = 'tests/data/test.gtf'
vcf = 'tests/data/test.vcf.gz'
fasta = 'tests/data/hg19.nochr.chr17.fa'
gtfIntervalTree = '../tests/data/test.pkl' # pickle exon interval Tree
# dataloader to load variants from vcf
dl = SplicingVCFDataloader(gtf,
fasta,
vcf,
out_file=gtfIntervalTree,
split_seq=False)
# Specify model
model = MMSplice(
exon_cut_l=0,
exon_cut_r=0,
acceptor_intron_cut=6,
donor_intron_cut=6,
acceptor_intron_len=50,
acceptor_exon_len=3,
donor_exon_len=5,
donor_intron_len=13)
# Do prediction
predictions = predict_all_table(model, dl, batch_size=1024, split_seq=False, assembly=False)
# Summerize with maximum effect size
predictionsMax = max_varEff(predictions)
```
=======
History
=======
0.1.0 (2018-07-17)
------------------
* First release on PyPI.
[![pypi](https://img.shields.io/pypi/v/mmsplice.svg)](https://pypi-hypernode.com/pypi/mmsplice)
[![travis](https://img.shields.io/travis/s6juncheng/mmsplice.svg)](https://travis-ci.org/s6juncheng/mmsplice)
Predict splicing variant effect from VCF
* Free software: MIT license
## Usage example
------
NOTE: make sure you split and left-normalize the input VCF file.
Check notebooks/example.ipynb
```python
# Import
from mmsplice.vcf_dataloader import SplicingVCFDataloader
from mmsplice import MMSplice, predict_all_table
from mmsplice.utils import max_varEff
# example files
gtf = 'tests/data/test.gtf'
vcf = 'tests/data/test.vcf.gz'
fasta = 'tests/data/hg19.nochr.chr17.fa'
gtfIntervalTree = '../tests/data/test.pkl' # pickle exon interval Tree
# dataloader to load variants from vcf
dl = SplicingVCFDataloader(gtf,
fasta,
vcf,
out_file=gtfIntervalTree,
split_seq=False)
# Specify model
model = MMSplice(
exon_cut_l=0,
exon_cut_r=0,
acceptor_intron_cut=6,
donor_intron_cut=6,
acceptor_intron_len=50,
acceptor_exon_len=3,
donor_exon_len=5,
donor_intron_len=13)
# Do prediction
predictions = predict_all_table(model, dl, batch_size=1024, split_seq=False, assembly=False)
# Summerize with maximum effect size
predictionsMax = max_varEff(predictions)
```
=======
History
=======
0.1.0 (2018-07-17)
------------------
* First release on PyPI.
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
mmsplice-0.2.1.tar.gz
(448.5 kB
view details)
Built Distribution
mmsplice-0.2.1-py2.py3-none-any.whl
(446.1 kB
view details)
File details
Details for the file mmsplice-0.2.1.tar.gz
.
File metadata
- Download URL: mmsplice-0.2.1.tar.gz
- Upload date:
- Size: 448.5 kB
- Tags: Source
- Uploaded using Trusted Publishing? No
- Uploaded via: twine/1.12.1 pkginfo/1.4.2 requests/2.18.4 setuptools/39.0.1 requests-toolbelt/0.8.0 tqdm/4.19.7 CPython/3.6.3
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | 86ef86e3b223acc35590165d507a852a7236237b2c9f7a8ded1620298db893dd |
|
MD5 | 052f74a0c986ad0838a6e7995331a76c |
|
BLAKE2b-256 | 04dc75b8214cea697276c1f1a512a3cafc8f135ccf8e85252b16a639cf9621da |
File details
Details for the file mmsplice-0.2.1-py2.py3-none-any.whl
.
File metadata
- Download URL: mmsplice-0.2.1-py2.py3-none-any.whl
- Upload date:
- Size: 446.1 kB
- Tags: Python 2, Python 3
- Uploaded using Trusted Publishing? No
- Uploaded via: twine/1.11.0 pkginfo/1.4.2 requests/2.19.1 setuptools/39.1.0 requests-toolbelt/0.8.0 tqdm/4.23.4 CPython/3.6.0
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | 02ba51c7ea161739240f3b5d9ab3750d3ab6fafb81669592114a983391fa0f1b |
|
MD5 | 129869919c99b62d6706c5b31acd3f09 |
|
BLAKE2b-256 | b1457fdb9377f6f1f73481b61cce298356fce875deea396c4f3a8793da05ab1c |