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
|
$ gemmi align -h
Pairwise sequence alignment with scoring matrix and affine gap penalty.
Usage:
gemmi align [options] FILE[...]
Aligns sequence from the model to the full sequence (SEQRES).
Both are from the same FILE - either in the PDB or mmCIF format.
If the mmCIF format is used, option --check-mmcif can be used.
gemmi align [options] --query=CHAIN1 --target=CHAIN2 FILE1 FILE2
Aligns CHAIN1 from FILE1 to CHAIN2 from FILE2.
By default, the sequence of residues in the model is used.
To use SEQRES prepend '+' to the chain name (e.g. --query=+A).
gemmi align [options] --text-align STRING1 STRING2
Aligns two ASCII strings (used for testing).
Options:
-h, --help Print usage and exit.
-V, --version Print version and exit.
--check-mmcif checks alignment against _atom_site.label_seq_id
--query=[+]CHAIN Align CHAIN from file INPUT1.
--target=[+]CHAIN Align CHAIN from file INPUT2.
--text-align Align characters in two strings (for testing).
Scoring (absolute values):
--match=INT Match score (default: 1).
--mism=INT Mismatch penalty (default: -1).
--gapo=INT Gap opening penalty (default: -1).
--gape=INT Gap extension penalty (default: -1).
Output options:
-p Print formatted alignment with one-letter codes.
--rmsd In addition to aligning two CHAINs (--query and --target),
superpose structures and print RMSD.
-v, --verbose Verbose output.
|
