forked from briney/nwalign3
-
Notifications
You must be signed in to change notification settings - Fork 0
/
PKG-INFO
90 lines (70 loc) · 3 KB
/
PKG-INFO
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
Metadata-Version: 1.0
Name: nwalign3
Version: 0.1.0
Summary: Needleman-Wunsch global sequence alignment
Home-page: https://gitub.com/briney/nwalign3
Author: Bryan Briney
Author-email: [email protected]
License: BSD
Description:
------------------------------------------------------------------------------
nwalign3: fast `cython`_ - `Needleman-Wunsch`_ alignment
------------------------------------------------------------------------------
.. _`Needleman-Wunsch`: http://en.wikipedia.org/wiki/Needleman-Wunsch_algorithm
.. _`scoring matrix`: http://en.wikipedia.org/wiki/Substitution_matrix
.. _`cython`: http://cython.org
This module provides a python module and a command-line interface to do global-
sequence alignment using the `Needleman-Wunsch` algorithm. It uses `cython`_
and numpy for speed.
Command-Line Usage
==================
the nwalign3 executable is installed to the PATH by setuptools
::
$ nwalign3 alphabet alpet
alphabet
alp---et
specify an alignment `scoring matrix`_
::
$ nwalign3 --matrix /usr/share/ncbi/data/BLOSUM62 EEAEE EEEEG
EEAEE-
EE-EEG
with specified penalties
::
$ nwalign3 --gap_open -10 --gap_extend -4 --match 12 ASDFF ASFF
ASDFF
AS-FF
Python Usage
============
Alignment
---------
::
>>> import nwalign3 as nw
>>> nw.global_align("CEELECANTH", "PELICAN", matrix='PAM250')
('CEELE-CANTH', '-PEL-ICAN--')
# with a specified penalty for open and extend.
>>> nw.global_align("CEELECANTH", "PELICAN", gap_open=-10, gap_extend=-4, matrix='PAM250')
('CEELECANTH', '-PELICAN--')
the `matrix` is specified as the full path to an `scoring matrix`_ as
is distributed with the NCBI toolset.
Scoring
-------
get the score of an alignment. (the first 2 args are from an alignment
and must have the same length.
::
>>> nw.score_alignment('CEELECANTH', '-PELICAN--', gap_open=-5,
... gap_extend=-2, matrix='PAM250')
11
>>> nw.score_alignment('CEELE-CANTH', '-PEL-ICAN--', gap_open=-5,
... gap_extend=-2, matrix='PAM250')
6
Keywords: sequence bioinformatics alignment needleman-wunsch
Platform: UNKNOWN
Classifier: Development Status :: 3 - Alpha
Classifier: Intended Audience :: Developers
Classifier: Intended Audience :: Science/Research
Classifier: License :: OSI Approved :: BSD License
Classifier: Operating System :: OS Independent
Classifier: Programming Language :: Python
Classifier: Topic :: Scientific/Engineering :: Bio-Informatics
Classifier: Topic :: Scientific/Engineering
Classifier: Topic :: Text Processing