pymugen 1.1

Management about FASTA and VCF files.

pymugen

This module allows to manage data about genomic sequences, specificaly getting sequences from a FASTA file.

The module’s other functionality is the ability to get sequences with a mutation from a VCF file and transform it to another representation.

Installation

Install pymugen with python

`pip3 install pymugen`

Usage/Examples

Fasta example

from pymugen.fasta import FastaReader

fasta = FastaReader("fasta.fa")

chromosome = fasta['chr1']
sequence = chromosome.sequence(10, 5, 5)

sequence[1]
>>> A

sequence[:]
>>> CATGCATGCAT

str(sequence)
>>> C-A-T-G-C,A,T-G-C-A-T

chromosome[5: 16]
>>> CATGCATGCAT

VCF example

from pymugen.transformers import ExtendedTransformer
from pymugen.fasta import FastaReader

vcf = ExtendedTransformer("vcf.vcf", "fasta.fa")

fasta = FastaReader("fasta.fa")

chromosome = fasta['chr1']
sequence = chromosome.sequence(10, 5, 5)

mutated = vcf.method(sequence, "CGT")

str(mutated)
>>> 'A|w-q-r-e-w,s-d-f,v-c-x-z-v'

mutated[:]
>>> 'wqrewsdfvcxzv'