Hi,
First thank you for developping IQtree :)
I am running into an issue and have been trying several things.
Here is my command line, it works well for an alignment file of 2000 sequences as an input but crashes with my final file of 86k sequences. I have both been trying it on my local computer and on our cluster.
Command line: iqtree -s <path to fasta file> -m K3Pu+F+G4 -bo 10
Seed: 815072 (Using SPRNG - Scalable Parallel Random Number Generator)
Part of the log with error (for the run on my local computer):
WARNING: 86622 sequences contain more than 50% gaps/ambiguity
**** TOTAL 56.00% 18 sequences failed composition chi2 test (p-value<5%; df=3)
===> START BOOTSTRAP REPLICATE NUMBER 1
Creating bootstrap alignment (seed: 815072)...
Create initial parsimony tree by phylogenetic likelihood library (PLL)... 2693.687 seconds
NOTE: 2112 MB RAM (2 GB) is required!
Estimate model parameters (epsilon = 0.100)
1. Initial log-likelihood: -325369.792
2. Current log-likelihood: -315330.127
3. Current log-likelihood: -314971.524
4. Current log-likelihood: -314969.954
Optimal log-likelihood: -314969.918
Rate parameters: A-C: 1.00000 A-G: 1.40015 A-T: 0.52119 C-G: 0.52119 C-T: 1.40015 G-T: 1.00000
Base frequencies: A: 0.254 C: 0.222 G: 0.232 T: 0.292
Gamma shape alpha: 0.845
Parameters optimization took 4 rounds (3081.983 sec)
Computing ML distances based on estimated model parameters...
ERROR: STACK TRACE FOR DEBUGGING:
ERROR: 1 funcAbort()
ERROR: 2 ()
ERROR: 3 gsignal()
ERROR: 4 abort()
ERROR: 5 __gnu_cxx::__verbose_terminate_handler()
ERROR: 6 ()
ERROR: 7 ()
ERROR: 8 __cxa_rethrow()
ERROR: 9 operator new(unsigned long)
ERROR: 10 PhyloTree::computeDist(Params&, Alignment*, double*&, double*&)
ERROR: 11 computeMLDist(Params&, IQTree&, double, double)
ERROR: 12 runTreeReconstruction(Params&, IQTree*&)
ERROR: 13 runStandardBootstrap(Params&, Alignment*, IQTree*)
ERROR: 14 runPhyloAnalysis(Params&, Checkpoint*)
ERROR: 15 main()
ERROR: 16 __libc_start_main()
ERROR: 17 ()
ERROR:
ERROR: *** IQ-TREE CRASHES WITH SIGNAL ABORTED
ERROR: *** For bug report please send to developers:
ERROR: *** Log file: Documents/TARACommunityDistributionNAP/WORK/rDNA_V9_18S_ribotypes_diatoms/Phylogeny/Phylogeny_Sina_ExaML_Diatoms_SO/diatoms_aligned_trimmed.fasta.log
ERROR: *** Alignment files (if possible)
And here is the error part when launching it on the cluster:
WARNING: 86622 sequences contain more than 50% gaps/ambiguity
**** TOTAL 56.00% 18 sequences failed composition chi2 test (p-value<5%; df=3)
===> START BOOTSTRAP REPLICATE NUMBER 1
Creating bootstrap alignment (seed: 600398)...
Create initial parsimony tree by phylogenetic likelihood library (PLL)... 8042.695 seconds
NOTE: 1977 MB RAM (1 GB) is required!
Estimate model parameters (epsilon = 0.100)
1. Initial log-likelihood: -418723.320
2. Current log-likelihood: -329717.488
3. Current log-likelihood: -329667.347
4. Current log-likelihood: -329665.307
5. Current log-likelihood: -329664.876
6. Current log-likelihood: -329664.761
Optimal log-likelihood: -329664.721
Rate parameters: A-C: 1.00000 A-G: 1.50455 A-T: 0.65018 C-G: 0.65018 C-T: 1.50455 G-T: 1.00000
Base frequencies: A: 0.227 C: 0.148 G: 0.383 T: 0.241
Gamma shape alpha: 1.762
Parameters optimization took 6 rounds (20665.391 sec)
Computing ML distances based on estimated model parameters...ERROR: STACK TRACE FOR DEBUGGING:
ERROR:
ERROR: *** IQ-TREE CRASHES WITH SIGNAL ABORTED
ERROR: *** For bug report please send to developers:
ERROR: *** Log file: output_data/run-37981/2021-07-02-diatoms_aligned_trimmed.fasta.MFP.0.log
ERROR: *** Alignment files (if possible)
On the cluster when it crashes the memory usage is around 3G but more is allocated.
On my PC I have 30G of RAM.
Do you know where this could be coming from?
Let me know if you need additionnal info, ie full log file or alignment file.
Best,
Mathilde