Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/moshi4/pymsaviz
MSA(Multiple Sequence Alignment) visualization python package for sequence analysis
https://github.com/moshi4/pymsaviz
bioinformatics genomics matplotlib msa multiple-sequence-alignment python sequence-alignment sequence-analysis visualization
Last synced: 4 days ago
JSON representation
MSA(Multiple Sequence Alignment) visualization python package for sequence analysis
- Host: GitHub
- URL: https://github.com/moshi4/pymsaviz
- Owner: moshi4
- License: mit
- Created: 2022-11-13T09:39:51.000Z (about 2 years ago)
- Default Branch: main
- Last Pushed: 2024-05-15T14:40:51.000Z (6 months ago)
- Last Synced: 2024-07-15T09:23:52.119Z (4 months ago)
- Topics: bioinformatics, genomics, matplotlib, msa, multiple-sequence-alignment, python, sequence-alignment, sequence-analysis, visualization
- Language: Python
- Homepage: https://moshi4.github.io/pyMSAviz
- Size: 7.07 MB
- Stars: 73
- Watchers: 4
- Forks: 11
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
- Citation: CITATION.cff
Awesome Lists containing this project
README
# pyMSAviz
[](https://pypi.python.org/pypi/pymsaviz)
[](https://anaconda.org/bioconda/pymsaviz)
[](https://github.com/moshi4/pyMSAviz/actions/workflows/ci.yml)
## Table of contents
- [Overview](#overview)
- [Installation](#installation)
- [API Usage](#api-usage)
- [CLI Usage](#cli-usage)
## Overview
pyMSAviz is a MSA(Multiple Sequence Alignment) visualization python package for sequence analysis implemented based on matplotlib.
This package is developed for the purpose of easily and beautifully plotting MSA in Python.
It also implements the functionality to add markers, text annotations, highlights to specific positions and ranges in MSA.
pyMSAviz was developed inspired by [Jalview](https://www.jalview.org/) and [ggmsa](https://github.com/YuLab-SMU/ggmsa).
More detailed documentation is available [here](https://moshi4.github.io/pyMSAviz/).
**Fig.1 Simple visualization result**
**Fig.2 Customized visualization result**
## Installation
`Python 3.9 or later` is required for installation.
**Install PyPI package:**
pip install pymsaviz
**Install bioconda package:**
conda install -c conda-forge -c bioconda pymsaviz
## API Usage
Only simple example usage is described in this section.
For more details, please see [Getting Started](https://moshi4.github.io/pyMSAviz/getting_started/) and [API Docs](https://moshi4.github.io/pyMSAviz/api-docs/msaviz/).
### API Example
#### API Example 1
```python
from pymsaviz import MsaViz, get_msa_testdata
msa_file = get_msa_testdata("HIGD2A.fa")
mv = MsaViz(msa_file, wrap_length=60, show_count=True)
mv.savefig("api_example01.png")
```
#### API Example 2
```python
from pymsaviz import MsaViz, get_msa_testdata
msa_file = get_msa_testdata("MRGPRG.fa")
mv = MsaViz(msa_file, color_scheme="Taylor", wrap_length=80, show_grid=True, show_consensus=True)
mv.savefig("api_example02.png")
```
#### API Example 3
```python
from pymsaviz import MsaViz, get_msa_testdata
msa_file = get_msa_testdata("MRGPRG.fa")
mv = MsaViz(msa_file, end=180, wrap_length=60, show_consensus=True)
# Extract MSA positions less than 50% consensus identity
pos_ident_less_than_50 = []
ident_list = mv._get_consensus_identity_list()
for pos, ident in enumerate(ident_list, 1):
if ident <= 50:
pos_ident_less_than_50.append(pos)
# Add markers
mv.add_markers([1])
mv.add_markers([10, 20], color="orange", marker="o")
mv.add_markers([30, (40, 50), 55], color="green", marker="+")
mv.add_markers(pos_ident_less_than_50, marker="x", color="blue")
# Add text annotations
mv.add_text_annotation((76, 102), "Gap Region", text_color="red", range_color="red")
mv.add_text_annotation((112, 123), "Gap Region", text_color="green", range_color="green")
mv.savefig("api_example03.png")
```
## CLI Usage
pyMSAviz provides simple MSA visualization CLI.
### Basic Command
pymsaviz -i [MSA file] -o [MSA visualization file]
### Options
$ pymsaviz --help
usage: pymsaviz [options] -i msa.fa -o msa_viz.png
MSA(Multiple Sequence Alignment) visualization CLI tool
optional arguments:
-i I, --infile I Input MSA file
-o O, --outfile O Output MSA visualization file (*.png|*.jpg|*.svg|*.pdf)
--format MSA file format (Default: 'fasta')
--color_scheme Color scheme (Default: 'Zappo')
--start Start position of MSA visualization (Default: 1)
--end End position of MSA visualization (Default: 'MSA Length')
--wrap_length Wrap length (Default: None)
--wrap_space_size Space size between wrap MSA plot area (Default: 3.0)
--label_type Label type ('id'[default]|'description')
--show_grid Show grid (Default: OFF)
--show_count Show seq char count without gap on right side (Default: OFF)
--show_consensus Show consensus sequence (Default: OFF)
--consensus_color Consensus identity bar color (Default: '#1f77b4')
--consensus_size Consensus identity bar height size (Default: 2.0)
--sort Sort MSA order by NJ tree constructed from MSA distance matrix (Default: OFF)
--dpi Figure DPI (Default: 300)
-v, --version Print version information
-h, --help Show this help message and exit
Available Color Schemes:
['Clustal', 'Zappo', 'Taylor', 'Flower', 'Blossom', 'Sunset', 'Ocean', 'Hydrophobicity', 'HelixPropensity', 'StrandPropensity', 'TurnPropensity', 'BuriedIndex', 'Nucleotide', 'Purine/Pyrimidine', 'Identity', 'None']
### CLI Example
Click [here](https://github.com/moshi4/pyMSAviz/raw/main/example/example.zip) to download example MSA files.
#### CLI Example 1
pymsaviz -i ./example/HIGD2A.fa -o cli_example01.png --color_scheme Identity
#### CLI Example 2
pymsaviz -i ./example/MRGPRG.fa -o cli_example02.png --wrap_length 80 \
--color_scheme Taylor --show_consensus --show_count
#### CLI Example 3
pymsaviz -i ./example/MRGPRG.fa -o cli_example03.png --start 100 --end 160 \
--color_scheme Flower --show_grid --show_consensus --consensus_color tomato
