Mutational spectra analysis package
Project description
PyMutSpec
Python library for mutational spectra analysis
Requirements
- python 3.8+
Installation
pip3 install pymutspec
https://pypi.org/project/PyMutSpec/
Example code
from Bio import SeqIO
from pymutspec.annotation import calculate_mutspec, CodonAnnotation
from pymutspec.draw import plot_mutspec12, plot_mutspec192
coda = CodonAnnotation(gencode=2) # mitochondrial genetic code
path_to_observed_mutations = ...
path_to_reference_seq = ...
# load data (mutations and sequence)
gene = SeqIO.parse(path_to_reference_seq, format='fasta')
observed_mutations = pd.read_csv(path_to_observed_mutations, sep='\t')
for col in ['Mut', 'MutType']:
assert col in observed_mutations.columns
# sample only syn mutations
mut_syn = observed_mutations[observed_mutations.MutType >= 1] # 0 for all mutations, 1 for syn, 2 for fourfold syn (syn4f)
# derive expected mutations from reference gene
sbs12_freqs, sbs192_freqs = coda.collect_exp_mut_freqs(gene, labels['all', 'syn', 'syn4f'])
sbs12_freqs_syn = sbs12_freqs['syn']
sbs192_freqs_syn = sbs192_freqs['syn']
# calculate mutation spectra
spectra12 = calculate_mutspec(mut_syn, sbs12_freqs_syn, use_context=False)
spectra192 = calculate_mutspec(mut_syn, sbs192_freqs_syn, use_context=True)
# plot mutation spectra
plot_mutspec12(spectra12)
plot_mutspec192(spectra192)
Spectra barplots
Links
- Iqtree - efficient software for phylogenomic inference
- Genetic codes
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
PyMutSpec-0.0.8.tar.gz
(41.2 kB
view hashes)
Built Distribution
PyMutSpec-0.0.8-py3-none-any.whl
(49.2 kB
view hashes)
Close
Hashes for PyMutSpec-0.0.8-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 71ca95665c7cf042b32720778b664a0888c20462ba00e27ca8e89e0f106b1d1f |
|
MD5 | b0083c262028fdaa658c9d115c3141d6 |
|
BLAKE2b-256 | 4895c3a0ebc73ce4211474f376ed8c6a9b1753d975ba5d8729ad8dfaaa270b72 |