failed in modelfinder

Skip to first unread message


Jun 28, 2019, 11:20:41 PM6/28/19
It works in 3 gene matrix,but failed in 62 genes. Is there anyone who have time to give me some suggestion? Thanks!

Host:    DESKTOP-4QC0DVO (AVX2, FMA3, 63 GB RAM)
Command: bin\iqtree -s n20dal.phy -m MF -nt AUTO
Seed:    134096 (Using SPRNG - Scalable Parallel Random Number Generator)
Time:    Fri Jun 28 16:03:17 2019
Kernel:  AVX+FMA - auto-detect threads (12 CPU cores detected)
Reading alignment file n20dal.phy ... Phylip format detected
Alignment most likely contains protein sequences
WARNING: 1 sites contain only gaps or ambiguous characters.
Alignment has 20 sequences with 47119 columns, 14598 distinct patterns
25073 parsimony-informative, 7627 singleton sites, 14419 constant sites
            Gap/Ambiguity  Composition  p-value
   1  FauNEOPT     13.33%    failed      0.00%
   2  LfoCHILO     13.10%    failed      0.00%
   3  LpoXIPHOS    13.13%    failed      0.00%
   4  MgaARACH     13.21%    failed      0.00%
   5  NamDIPLO     13.51%    failed      0.00%
   6  OstOST       13.01%    failed      0.00%
   7  PaqCOLL      13.37%    failed      0.11%
   8  StuREMI      13.20%    failed      0.00%
   9  DmaBRANCH    13.10%    passed     84.89%
  10  HmaCEPHAL    12.32%    failed      0.00%
  11  ScolCHILO    13.37%    failed      0.00%
  12  SpoCHILO     13.44%    failed      0.00%
  13  AmaDIPLO     13.26%    passed     99.24%
  14  Han2SYMPH    13.42%    failed      0.00%
  15  Pge2DIPLO    12.82%    failed      0.00%
  16  Pol2DIPLO    13.04%    failed      0.00%
  17  Scu3SYMPH    13.21%    failed      0.00%
  18  Col2PYCNO    12.93%    failed      0.00%
  19  Ctas2CHILO   13.32%    failed      0.00%
  20  EuryPAURO    13.32%    failed      0.00%
****  TOTAL        13.17%  18 sequences failed composition chi2 test (p-value<5%; df=19)

Create initial parsimony tree by phylogenetic likelihood library (PLL)... 0.109 seconds
NOTE: ModelFinder requires 451 MB RAM!
ModelFinder will test 546 protein models (sample size: 47119) ...
 No. Model         -LnL         df  AIC          AICc         BIC
Measuring multi-threading efficiency up to 12 CPU cores
Increase to 10 rounds for branch lengths
21 trees examined
Threads: 1 / Time: 12.341 sec / Speedup: 1.000 / Efficiency: 100% / LogL: -935133
new_tree_lh: -865423.1692921221   tree_lh: -863818.5083998179
ERROR: C:\Users\minh\Dropbox\iqtree-git\tree\phylotree.cpp:2652: virtual double PhyloTree::optimizeAllBranches(int, double, int): Assertion `fabs(new_tree_lh-tree_lh) < max_delta_lh' failed.
ERROR: *** For bug report please send to developers:
ERROR: ***    Log file: n20dal.phy.log
ERROR: ***    Alignment files (if possible)

Minh Bui

Jun 30, 2019, 3:37:28 PM6/30/19
to IQ-TREE, YanYan
Hi YanYan,

I could not replicate this crash using the latest version 1.6.11. 

Which version did you use? Can you pls try the latest version.


You received this message because you are subscribed to the Google Groups "IQ-TREE" group.
To unsubscribe from this group and stop receiving emails from it, send an email to
To post to this group, send email to
Visit this group at
To view this discussion on the web visit
For more options, visit


Jul 1, 2019, 3:53:43 AM7/1/19

Dear Minh:

Thank you very much to give me suggestion so soon.

Yes I using latest version 1.6.11, and test in win10 and win7.

The n20dal.phy file I tried to analysis is a degenerated matrix, that is a coding method for nucleotides in which codons in a data matrix are completely degenerated, for the first and third codon positions using ambiguity coding e.g., CAT -- . CAY and TTA--. YTN.

But IQ-tree detected: Alignment most likely contains protein sequences (not DNA), and try to use protein model to analysis and then failed.

I do not know whether it is the reason that –m MF requirement failed.

I test an amino acid matrix. It works after I enter iqtree –s aamatrix.phy

Is it possible to let IQ-tree recognize degenerated DNA matrix?

在 2019年6月29日星期六 UTC+8上午11:20:41,YanYan写道:

Minh Bui

Jul 1, 2019, 6:06:18 PM7/1/19
to IQ-TREE, YanYan
Ah I see. That’s the reason why IQ-TREE mis-identify your sequences as protein. You can use -st DNA to enforce DNA sequences, instead.


You received this message because you are subscribed to the Google Groups "IQ-TREE" group.
To unsubscribe from this group and stop receiving emails from it, send an email to
To post to this group, send email to
Visit this group at


Jul 2, 2019, 3:49:14 AM7/2/19
Dear Minh:

After adding -st DNA, IQ-tree works! Thank you! I am very happy that it is really fast.
Reply all
Reply to author
0 new messages