R version: /usr/local/share/applications/R-3.4.3/bin/R homo: Aesc_assa_HG013_c98.Aesc_assa_HG013_c98.homo dir: Aesc_assa_HG013_c98.Aesc_assa_HG013_c98 cds1: Aesc_assa_HG013_c98.cds cds2: Aesc_assa_HG013_c98.cds sp1: Aesc_assa_HG013_c98 sp2: Aesc_assa_HG013_c98 2020-05-13 06:55:56: INFO 2020-05-13 06:55:56: INFO codeml found 2020-05-13 06:55:56: INFO 2020-05-13 06:55:56: INFO 2020-05-13 06:55:57: INFO Translating CDS file N/A% (0 of 42195) | | Elapsed Time: 0:00:00 ETA: --:--:-- 0% (263 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 1% (504 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:16 1% (785 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 2% (1043 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 3% (1354 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 3% (1603 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 4% (1856 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 5% (2137 of 42195) | | Elapsed Time: 0:00:00 ETA: 0:00:15 6% (2634 of 42195) |# | Elapsed Time: 0:00:00 ETA: 0:00:13 6% (2862 of 42195) |# | Elapsed Time: 0:00:00 ETA: 0:00:13 7% (3125 of 42195) |# | Elapsed Time: 0:00:01 ETA: 0:00:13 8% (3419 of 42195) |# | Elapsed Time: 0:00:01 ETA: 0:00:13 8% (3739 of 42195) |# | Elapsed Time: 0:00:01 ETA: 0:00:13 9% (3999 of 42195) |# | Elapsed Time: 0:00:01 ETA: 0:00:13 10% (4273 of 42195) |# | Elapsed Time: 0:00:01 ETA: 0:00:13 10% (4534 of 42195) |## | Elapsed Time: 0:00:01 ETA: 0:00:12 11% (4808 of 42195) |## | Elapsed Time: 0:00:01 ETA: 0:00:12 11% (5048 of 42195) |## | Elapsed Time: 0:00:01 ETA: 0:00:12 12% (5282 of 42195) |## | Elapsed Time: 0:00:01 ETA: 0:00:12 13% (5587 of 42195) |## | Elapsed Time: 0:00:01 ETA: 0:00:12 13% (5876 of 42195) |## | Elapsed Time: 0:00:02 ETA: 0:00:12 14% (6160 of 42195) |## | Elapsed Time: 0:00:02 ETA: 0:00:12 15% (6410 of 42195) |## | Elapsed Time: 0:00:02 ETA: 0:00:12 15% (6654 of 42195) |## | Elapsed Time: 0:00:02 ETA: 0:00:12 16% (6934 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 17% (7177 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 17% (7447 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 18% (7671 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 18% (7889 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 18% (8012 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 19% (8358 of 42195) |### | Elapsed Time: 0:00:02 ETA: 0:00:12 20% (8546 of 42195) |### | Elapsed Time: 0:00:03 ETA: 0:00:12 20% (8810 of 42195) |### | Elapsed Time: 0:00:03 ETA: 0:00:12 21% (9065 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:12 22% (9484 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 23% (9777 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 23% (9976 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 24% (10149 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 24% (10366 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 25% (10603 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 25% (10892 of 42195) |#### | Elapsed Time: 0:00:03 ETA: 0:00:11 26% (11166 of 42195) |#### | Elapsed Time: 0:00:04 ETA: 0:00:11 27% (11470 of 42195) |#### | Elapsed Time: 0:00:04 ETA: 0:00:11 27% (11747 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:11 28% (12040 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:11 29% (12279 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:11 29% (12566 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:10 30% (12819 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:10 31% (13137 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:10 31% (13353 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:10 32% (13591 of 42195) |##### | Elapsed Time: 0:00:04 ETA: 0:00:10 32% (13854 of 42195) |##### | Elapsed Time: 0:00:05 ETA: 0:00:10 33% (14162 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:10 34% (14422 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:10 34% (14698 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:10 35% (14956 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:09 35% (15167 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:10 36% (15410 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:10 37% (15662 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:10 37% (15972 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:09 38% (16309 of 42195) |###### | Elapsed Time: 0:00:05 ETA: 0:00:09 39% (16558 of 42195) |####### | Elapsed Time: 0:00:05 ETA: 0:00:09 39% (16756 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 40% (16947 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 40% (17092 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 41% (17362 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 41% (17570 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 42% (17820 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 42% (18068 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 43% (18309 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 44% (18585 of 42195) |####### | Elapsed Time: 0:00:06 ETA: 0:00:09 44% (18796 of 42195) |######## | Elapsed Time: 0:00:06 ETA: 0:00:09 45% (19036 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:09 45% (19229 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 46% (19510 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 46% (19757 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 47% (19986 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 47% (20217 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 48% (20485 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 49% (20752 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 49% (21042 of 42195) |######## | Elapsed Time: 0:00:07 ETA: 0:00:08 50% (21298 of 42195) |######### | Elapsed Time: 0:00:07 ETA: 0:00:08 50% (21496 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:08 51% (21747 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:08 51% (21899 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:08 52% (22158 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:08 53% (22418 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:07 53% (22637 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:07 54% (22872 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:07 54% (23127 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:07 55% (23330 of 42195) |######### | Elapsed Time: 0:00:08 ETA: 0:00:07 55% (23502 of 42195) |########## | Elapsed Time: 0:00:08 ETA: 0:00:07 56% (23755 of 42195) |########## | Elapsed Time: 0:00:08 ETA: 0:00:07 56% (24015 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:07 57% (24236 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:07 57% (24467 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:07 58% (24716 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:07 59% (24993 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:07 59% (25207 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:06 60% (25412 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:06 60% (25593 of 42195) |########## | Elapsed Time: 0:00:09 ETA: 0:00:07 61% (25811 of 42195) |########### | Elapsed Time: 0:00:09 ETA: 0:00:07 61% (26068 of 42195) |########### | Elapsed Time: 0:00:09 ETA: 0:00:06 62% (26393 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 62% (26563 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 63% (26706 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 63% (26957 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 64% (27156 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 64% (27240 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 65% (27512 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 65% (27774 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 66% (27966 of 42195) |########### | Elapsed Time: 0:00:10 ETA: 0:00:06 66% (28132 of 42195) |############ | Elapsed Time: 0:00:10 ETA: 0:00:06 67% (28309 of 42195) |############ | Elapsed Time: 0:00:10 ETA: 0:00:06 67% (28563 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 68% (28746 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 68% (28992 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 69% (29226 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 69% (29377 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 70% (29592 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 70% (29862 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 71% (30071 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 71% (30249 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 72% (30427 of 42195) |############ | Elapsed Time: 0:00:11 ETA: 0:00:05 72% (30595 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:05 73% (30852 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:05 73% (31112 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:05 74% (31331 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 74% (31513 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 75% (31669 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 75% (31863 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 75% (32037 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 76% (32263 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 76% (32423 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 77% (32581 of 42195) |############# | Elapsed Time: 0:00:12 ETA: 0:00:04 77% (32778 of 42195) |############# | Elapsed Time: 0:00:13 ETA: 0:00:04 78% (32987 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 78% (33116 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 79% (33362 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 79% (33518 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 79% (33650 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 80% (33796 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 80% (33972 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 80% (34142 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 81% (34329 of 42195) |############## | Elapsed Time: 0:00:13 ETA: 0:00:04 81% (34469 of 42195) |############## | Elapsed Time: 0:00:14 ETA: 0:00:04 82% (34634 of 42195) |############## | Elapsed Time: 0:00:14 ETA: 0:00:04 82% (34795 of 42195) |############## | Elapsed Time: 0:00:14 ETA: 0:00:04 82% (34924 of 42195) |############## | Elapsed Time: 0:00:14 ETA: 0:00:04 83% (35077 of 42195) |############## | Elapsed Time: 0:00:14 ETA: 0:00:04 83% (35214 of 42195) |############### | Elapsed Time: 0:00:14 ETA: 0:00:04 83% (35378 of 42195) |############### | Elapsed Time: 0:00:14 ETA: 0:00:03 84% (35532 of 42195) |############### | Elapsed Time: 0:00:14 ETA: 0:00:03 84% (35653 of 42195) |############### | Elapsed Time: 0:00:14 ETA: 0:00:03 84% (35786 of 42195) |############### | Elapsed Time: 0:00:14 ETA: 0:00:03 85% (35884 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 85% (36039 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 85% (36181 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 86% (36320 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 86% (36496 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 86% (36671 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 87% (36840 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 88% (37214 of 42195) |############### | Elapsed Time: 0:00:15 ETA: 0:00:03 89% (37567 of 42195) |################ | Elapsed Time: 0:00:15 ETA: 0:00:02 89% (37863 of 42195) |################ | Elapsed Time: 0:00:15 ETA: 0:00:02 90% (38323 of 42195) |################ | Elapsed Time: 0:00:16 ETA: 0:00:02 92% (38926 of 42195) |################ | Elapsed Time: 0:00:16 ETA: 0:00:01 93% (39525 of 42195) |################ | Elapsed Time: 0:00:16 ETA: 0:00:01 94% (40059 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 96% (40593 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 97% (41067 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 98% (41401 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 98% (41644 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 99% (41903 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 99% (42172 of 42195) |################# | Elapsed Time: 0:00:16 ETA: 0:00:00 100% (42195 of 42195) |##################| Elapsed Time: 0:00:16 Time: 0:00:16 2020-05-13 06:56:14: WARNING There were 0 warnings during translation 2020-05-13 06:56:14: INFO Started whole paranome Ks analysis 2020-05-13 06:56:14: WARNING Filtered out the 2 largest gene families because n*(n-1)/2 > `max_pairwise` 2020-05-13 06:56:14: WARNING If you want to analyse these large families anyhow, please raise the `max_pairwise` parameter. 2020-05-13 06:56:14: INFO Started analysis in parallel (n_threads = 1) 2020-05-13 06:56:14: INFO Performing analysis on gene family GF_000003 2020-05-13 06:56:36: INFO Performing analysis on gene family GF_000004 2020-05-13 06:56:43: INFO Performing analysis on gene family GF_000005 2020-05-13 06:57:50: INFO Performing analysis on gene family GF_000006 2020-05-13 06:58:07: INFO Performing analysis on gene family GF_000007 2020-05-13 06:58:34: INFO Performing analysis on gene family GF_000008 2020-05-13 06:59:13: INFO Performing analysis on gene family GF_000009 2020-05-13 06:59:18: INFO Performing analysis on gene family GF_000010 2020-05-13 06:59:40: INFO Performing analysis on gene family GF_000011 2020-05-13 06:59:55: INFO Performing analysis on gene family GF_000012 2020-05-13 06:59:57: INFO Performing analysis on gene family GF_000013 2020-05-13 07:00:03: INFO Performing analysis on gene family GF_000014 2020-05-13 07:00:10: INFO Performing analysis on gene family GF_000015 2020-05-13 07:00:12: INFO Performing analysis on gene family GF_000016 2020-05-13 07:00:29: INFO Performing analysis on gene family GF_000017 2020-05-13 07:00:36: INFO Performing analysis on gene family GF_000018 2020-05-13 07:00:38: INFO Performing analysis on gene family GF_000019 2020-05-13 07:00:42: INFO Performing analysis on gene family GF_000020 2020-05-13 07:00:47: INFO Performing analysis on gene family GF_000021 2020-05-13 07:00:53: INFO Performing analysis on gene family GF_000022 2020-05-13 07:01:04: WARNING No Ks value for Aesc_assa_c98.83532 - Aesc_assa_c98.10806! 2020-05-13 07:01:04: WARNING No Ka value for Aesc_assa_c98.83532 - Aesc_assa_c98.10806! 2020-05-13 07:01:04: WARNING No ω value for Aesc_assa_c98.83532 - Aesc_assa_c98.10806! Traceback (most recent call last): File "/home/kuokuo/anaconda3/envs/python36/bin/wgd", line 8, in sys.exit(cli()) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/click/core.py", line 764, in __call__ return self.main(*args, **kwargs) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/click/core.py", line 717, in main rv = self.invoke(ctx) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/click/core.py", line 1137, in invoke return _process_result(sub_ctx.command.invoke(sub_ctx)) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/click/core.py", line 956, in invoke return ctx.invoke(self.callback, **ctx.params) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/click/core.py", line 555, in invoke return callback(*args, **kwargs) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/wgd_cli.py", line 632, in ksd max_pairwise=max_pairwise File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/wgd_cli.py", line 773, in ksd_ max_pairwise=max_pairwise, File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/wgd/ks_distribution.py", line 645, in ks_analysis_paranome ) for family in sorted_families) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/parallel.py", line 779, in __call__ while self.dispatch_one_batch(iterator): File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/parallel.py", line 625, in dispatch_one_batch self._dispatch(tasks) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/parallel.py", line 588, in _dispatch job = self._backend.apply_async(batch, callback=cb) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/_parallel_backends.py", line 111, in apply_async result = ImmediateResult(func) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/_parallel_backends.py", line 332, in __init__ self.results = batch() File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/parallel.py", line 131, in __call__ return [func(*args, **kwargs) for func, args, kwargs in self.items] File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/joblib/parallel.py", line 131, in return [func(*args, **kwargs) for func, args, kwargs in self.items] File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/wgd/ks_distribution.py", line 305, in analyse_family results_dict, msa=msa_path_protein, method=method) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/wgd/ks_distribution.py", line 98, in _weighting tree_path, pairwise_estimates['Ks']) File "/home/kuokuo/anaconda3/envs/python36/lib/python3.6/site-packages/wgd/phy.py", line 138, in phylogenetic_tree_to_cluster_format node.name = id_map[node.name] KeyError: 'Aesc_assa_c98.54211'