-
Notifications
You must be signed in to change notification settings - Fork 1.7k
/
_ClustalOmega.py
184 lines (165 loc) · 8.26 KB
/
_ClustalOmega.py
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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
# -*- coding: utf-8 -*-
# Copyright 2011 by Andreas Wilm. All rights reserved.
# Based on ClustalW wrapper copyright 2009 by Cymon J. Cox.
#
# Wrapper for Clustal Omega by Andreas Wilm (2011). Used _Clustalw.py
# as template.
#
# This code is part of the Biopython distribution and governed by its
# license. Please see the LICENSE file that should have been included
# as part of this package.
"""Command line wrapper for the multiple alignment program Clustal Omega.
"""
from Bio.Application import _Option, _Switch, AbstractCommandline
class ClustalOmegaCommandline(AbstractCommandline):
"""Command line wrapper for clustal omega
http://www.clustal.org/omega
Example:
>>> from Bio.Align.Applications import ClustalOmegaCommandline
>>> in_file = "unaligned.fasta"
>>> out_file = "aligned.fasta"
>>> clustalomega_cline = ClustalOmegaCommandline(infile=in_file, outfile=out_file, verbose=True, auto=True)
>>> print clustalomega_cline
clustalo -i unaligned.fasta -o aligned.fasta --auto -v
You would typically run the command line with clustalomega_cline() or via
the Python subprocess module, as described in the Biopython tutorial.
Citation:
Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, Lopez R,
McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG (2011).
Fast, scalable generation of high-quality protein multiple
sequence alignments using Clustal Omega.
Molecular Systems Biology 7:539 doi:10.1038/msb.2011.75
Last checked against versions: 1.1.0
"""
def __init__(self, cmd="clustalo", **kwargs):
# order parameters in the same order as clustalo --help
self.parameters = \
[
# Sequence Input
_Option(["-i", "--in", "--infile", "infile"],
"Multiple sequence input file",
filename=True,
equate=False),
_Option(["--hmm-in", "HMM input", "hmm_input"],
"HMM input files",
filename=True,
equate=False),
_Switch(["--dealign", "dealign"],
"Dealign input sequences"),
_Option(["--profile1", "--p1", "profile1"],
"Pre-aligned multiple sequence file (aligned columns will be kept fix).",
filename=True,
equate=False),
_Option(["--profile2", "--p2", "profile2"],
"Pre-aligned multiple sequence file (aligned columns will be kept fix).",
filename=True,
equate=False),
_Option(["-t", "--seqtype", "seqtype"],
"{Protein, RNA, DNA} Force a sequence type (default: auto).",
equate=False,
checker_function=lambda x: x in ["protein", "rna", "dna",
"Protein", "RNA", "DNA",
"PROTEIN"]),
_Option(["--infmt", "infmt"],
"""Forced sequence input file format (default: auto)
Allowed values: a2m, fa[sta], clu[stal], msf, phy[lip], selex, st[ockholm], vie[nna]
""",
equate=False,
checker_function=lambda x: x in ["a2m", "fa", "fasta",
"clu", "clustal",
"msf",
"phy", "phylip",
"selex",
"st", "stockholm",
"vie", "vienna"]),
# Clustering
_Option(["--distmat-in", "distmat_in"],
"Pairwise distance matrix input file (skips distance computation).",
filename=True,
equate=False),
_Option(["--distmat-out", "distmat_out"],
"Pairwise distance matrix output file.",
filename=True,
equate=False),
_Option(["--guidetree-in", "guidetree_in"],
"Guide tree input file (skips distance computation and guide-tree clustering step).",
filename=True,
equate=False),
_Option(["--guidetree-out", "guidetree_out"],
"Guide tree output file.",
filename=True,
equate=False),
_Switch(["--full", "distmat_full"],
"Use full distance matrix for guide-tree calculation (might be slow; mBed is default)"),
_Switch(["--full-iter", "distmat_full_iter"],
"Use full distance matrix for guide-tree calculation during iteration (might be slowish; mBed is default)"),
# Alignment Output
_Option(["-o", "--out", "--outfile", "outfile"],
"Multiple sequence alignment output file (default: stdout).",
filename=True,
equate=False),
_Option(["--outfmt", "outfmt"],
"MSA output file format:"
" a2m=fa[sta],clu[stal],msf,phy[lip],selex,st[ockholm],vie[nna]"
" (default: fasta).",
equate=False,
checker_function=lambda x: x in ["a2m", "fa", "fasta",
"clu", "clustal",
"msf",
"phy", "phylip",
"selex",
"st", "stockholm",
"vie", "vienna"]),
# Iteration
_Option(["--iterations", "--iter", "iterations"],
"Number of (combined guide-tree/HMM) iterations",
equate=False,
checker_function=lambda x: isinstance(x, int)),
_Option(["--max-guidetree-iterations", "max_guidetree_iterations"],
"Maximum number of guidetree iterations",
equate=False,
checker_function=lambda x: isinstance(x, int)),
_Option(["--max-hmm-iterations", "max_hmm_iterations"],
"Maximum number of HMM iterations",
equate=False,
checker_function=lambda x: isinstance(x, int)),
# Limits (will exit early, if exceeded):
_Option(["--maxnumseq", "maxnumseq"],
"Maximum allowed number of sequences",
equate=False,
checker_function=lambda x: isinstance(x, int)),
_Option(["--maxseqlen", "maxseqlen"],
"Maximum allowed sequence length",
equate=False,
checker_function=lambda x: isinstance(x, int)),
# Miscellaneous:
_Switch(["--auto", "auto"],
"Set options automatically (might overwrite some of your options)"),
_Option(["--threads", "threads"],
"Number of processors to use",
equate=False,
checker_function=lambda x: isinstance(x, int)),
_Option(["-l", "--log", "log"],
"Log all non-essential output to this file.",
filename=True,
equate=False),
_Switch(["-h", "--help", "help"],
"Outline the command line params."),
_Switch(["-v", "--verbose", "verbose"],
"Verbose output"),
_Switch(["--version", "version"],
"Print version information and exit"),
_Switch(["--long-version", "long_version"],
"Print long version information and exit"),
_Switch(["--force", "force"],
"Force file overwriting."),
]
AbstractCommandline.__init__(self, cmd, **kwargs)
def _test():
"""Run the module's doctests (PRIVATE)."""
print "Running ClustalOmega doctests..."
import doctest
doctest.testmod()
print "Done"
if __name__ == "__main__":
_test()