| line |
stmt |
bran |
cond |
sub |
pod |
time |
code |
|
1
|
|
|
|
|
|
|
#!/usr/bin/env perl |
|
2
|
|
|
|
|
|
|
|
|
3
|
29
|
|
|
29
|
|
148855
|
use strict; |
|
|
29
|
|
|
|
|
82
|
|
|
|
29
|
|
|
|
|
1220
|
|
|
4
|
29
|
|
|
29
|
|
135
|
use warnings; |
|
|
29
|
|
|
|
|
49
|
|
|
|
29
|
|
|
|
|
1669
|
|
|
5
|
|
|
|
|
|
|
|
|
6
|
29
|
|
|
29
|
|
386
|
use v5.10; |
|
|
29
|
|
|
|
|
115
|
|
|
7
|
29
|
|
|
29
|
|
21438
|
use Getopt::Long qw( :config bundling permute no_getopt_compat ); |
|
|
29
|
|
|
|
|
472158
|
|
|
|
29
|
|
|
|
|
236
|
|
|
8
|
29
|
|
|
29
|
|
23351
|
use Pod::Usage; |
|
|
29
|
|
|
|
|
2114880
|
|
|
|
29
|
|
|
|
|
4434
|
|
|
9
|
29
|
|
|
29
|
|
26404
|
use Path::Tiny; |
|
|
29
|
|
|
|
|
480341
|
|
|
|
29
|
|
|
|
|
2942
|
|
|
10
|
29
|
|
|
29
|
|
17101
|
use lib path($0)->absolute->parent->sibling('lib')->stringify; |
|
|
29
|
|
|
|
|
21987
|
|
|
|
29
|
|
|
|
|
177
|
|
|
11
|
|
|
|
|
|
|
|
|
12
|
29
|
|
|
29
|
|
12391
|
use Pod::Usage; |
|
|
29
|
|
|
|
|
94
|
|
|
|
29
|
|
|
|
|
3375
|
|
|
13
|
29
|
|
|
29
|
|
14188
|
use Bio::BPWrapper; |
|
|
29
|
|
|
|
|
155
|
|
|
|
29
|
|
|
|
|
978
|
|
|
14
|
29
|
|
|
29
|
|
18076
|
use Bio::BPWrapper::AlnManipulations; |
|
|
29
|
|
|
|
|
161
|
|
|
|
29
|
|
|
|
|
12675
|
|
|
15
|
|
|
|
|
|
|
|
|
16
|
|
|
|
|
|
|
################################################################################ |
|
17
|
|
|
|
|
|
|
# Option parsing |
|
18
|
|
|
|
|
|
|
################################################################################ |
|
19
|
29
|
|
|
|
|
2221668
|
my %opts; |
|
20
|
29
|
50
|
|
|
|
371
|
GetOptions( |
|
21
|
|
|
|
|
|
|
\%opts, |
|
22
|
|
|
|
|
|
|
"help|h", |
|
23
|
|
|
|
|
|
|
"man", |
|
24
|
|
|
|
|
|
|
"version|V", |
|
25
|
|
|
|
|
|
|
"aln-index|I=s", |
|
26
|
|
|
|
|
|
|
"avg-pid|a", |
|
27
|
|
|
|
|
|
|
"binary", # 0/1 for binary sites (for phylip legacy programs) |
|
28
|
|
|
|
|
|
|
"bin-inform", # extract binary & informative sites (for clique): discard constant, 3/4-states, non-informative |
|
29
|
|
|
|
|
|
|
"bin-ref=s", |
|
30
|
|
|
|
|
|
|
"boot|b", |
|
31
|
|
|
|
|
|
|
"codon-view|c:20", # Optional value, default 20 |
|
32
|
|
|
|
|
|
|
"con-blocks|B:6", |
|
33
|
|
|
|
|
|
|
"concat|A", |
|
34
|
|
|
|
|
|
|
"consensus|C:50", |
|
35
|
|
|
|
|
|
|
"delete|d=s", |
|
36
|
|
|
|
|
|
|
"dna2pep|D", |
|
37
|
|
|
|
|
|
|
"gap-char=s", |
|
38
|
|
|
|
|
|
|
"gap-states", |
|
39
|
|
|
|
|
|
|
"gap-states2", |
|
40
|
|
|
|
|
|
|
"input|i=s", |
|
41
|
|
|
|
|
|
|
"length|l", |
|
42
|
|
|
|
|
|
|
"list-ids|L", |
|
43
|
|
|
|
|
|
|
"match|m", |
|
44
|
|
|
|
|
|
|
"no-flat|F", |
|
45
|
|
|
|
|
|
|
"no-gaps|g", |
|
46
|
|
|
|
|
|
|
"num-seq|n", |
|
47
|
|
|
|
|
|
|
"output|o=s", |
|
48
|
|
|
|
|
|
|
"pair-diff", # pairwise sequence diff |
|
49
|
|
|
|
|
|
|
"pair-diff-ref=s", # pairwise sequence diff to a ref seq |
|
50
|
|
|
|
|
|
|
"pep2dna|P=s", |
|
51
|
|
|
|
|
|
|
"permute-states|M", |
|
52
|
|
|
|
|
|
|
"phy-nonint", # non-interleaved phylip (for e.g.,clique) |
|
53
|
|
|
|
|
|
|
"pick|p=s", |
|
54
|
|
|
|
|
|
|
"random-slice=i", |
|
55
|
|
|
|
|
|
|
"ref-seq|r=s", |
|
56
|
|
|
|
|
|
|
"remove-third", |
|
57
|
|
|
|
|
|
|
"resample|R:i", # Optional value, default is floor of num_sequences/2 |
|
58
|
|
|
|
|
|
|
"rm-col|E=s", |
|
59
|
|
|
|
|
|
|
"select-third", |
|
60
|
|
|
|
|
|
|
"shuffle-sites|S", |
|
61
|
|
|
|
|
|
|
"slice|s=s", |
|
62
|
|
|
|
|
|
|
"split-cdhit=s", |
|
63
|
|
|
|
|
|
|
"trim-ends", |
|
64
|
|
|
|
|
|
|
"uniq|u", |
|
65
|
|
|
|
|
|
|
"upper", # make upper case (for DNAStatistics) |
|
66
|
|
|
|
|
|
|
"var-sites|v", |
|
67
|
|
|
|
|
|
|
"window|w:30", |
|
68
|
|
|
|
|
|
|
# "dnadist|D=s", # Needs fixing |
|
69
|
|
|
|
|
|
|
# "inform|Z", |
|
70
|
|
|
|
|
|
|
) or pod2usage(2); |
|
71
|
|
|
|
|
|
|
|
|
72
|
29
|
|
|
29
|
|
292
|
use constant PROGRAM => File::Basename::basename(__FILE__); |
|
|
29
|
|
|
|
|
58
|
|
|
|
29
|
|
|
|
|
3416867
|
|
|
73
|
29
|
50
|
|
|
|
96844
|
Bio::BPWrapper::print_version(PROGRAM) if $opts{"version"}; |
|
74
|
|
|
|
|
|
|
|
|
75
|
29
|
|
|
|
|
231
|
initialize(\%opts); |
|
76
|
29
|
|
|
|
|
17743
|
write_out(\%opts); |
|
77
|
|
|
|
|
|
|
|
|
78
|
|
|
|
|
|
|
################# POD Documentation ############ |
|
79
|
|
|
|
|
|
|
|
|
80
|
|
|
|
|
|
|
__END__ |