Python package for converting sequence annotations to gene feature enumerations (GFE).


Keywords
pygfe, gfe
License
LGPL-3.0
Install
pip install py-gfe==1.1.6

Documentation

py-gfe

Documentation Status Updates

Python Boilerplate contains all the boilerplate you need to create a Python package.

Docker

docker pull nmdpbioinformatics/py-gfe
docker run -it --rm -v $PWD:/opt nmdpbioinformatics/py-gfe seq2gfe \
        -f /opt/your_fastafile.fasta -l HLA-A

Example

>>> from Bio import SeqIO
>>> from BioSQL import BioSeqDatabase
>>> from seqann.sequence_annotation import BioSeqAnn
>>> import pygfe
>>> seq_file = 'test_dq.fasta'
>>> gfe = pygfe.pyGFE()
>>> server = BioSeqDatabase.open_database(driver="pymysql", user="root",
...                                       passwd="", host="localhost",
...                                      db="bioseqdb")
>>> seqann = BioSeqAnn(server=server)
>>> seq_rec = list(SeqIO.parse(seq_file, 'fasta'))[0]
>>> annotation = seqann.annotate(seq_rec, "HLA-DQB1")
>>> gfe = gfe.get_gfe(annotation, "HLA-DQB1")
>>> print(gfe)
HLA-DQB1w0-4-0-141-0-12-0-4-0-0-0-0-0

Credits

This package was created with Cookiecutter and the audreyr/cookiecutter-pypackage project template.