Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

N<4, iqtree: ERROR: It makes no sense to perform bootstrap with less than 4 sequences. #85

Closed
cmkobel opened this issue Feb 2, 2024 · 1 comment
Labels
wontfix This will not be worked on

Comments

@cmkobel
Copy link
Owner

cmkobel commented Feb 2, 2024

Either there should be an if statement that shuts bootstrapping off when N<4, or this error should be simply ignored. I guess FastTree is good enough when N<4.

Activating conda environment: ../../../../../../conda-envs/108c93bf8519b05838cdca347b9bf6fa
IQ-TREE multicore version 2.2.5 COVID-edition for Linux 64-bit built Sep 15 2023
Developed by Bui Quang Minh, James Barbetti, Nguyen Lam Tung,
Olga Chernomor, Heiko Schmidt, Dominik Schrempf, Michael Woodhams, Ly Trong Nhan.

Host:    threadripper (AVX2, FMA3, 503 GB RAM)
Command: iqtree -s results_ac2/roary/core_gene_alignment.aln -m GTR --boot 100 --prefix results_ac2/iqtree/core_genome_iqtree -redo
Seed:    488167 (Using SPRNG - Scalable Parallel Random Number Generator)
Time:    Fri Feb  2 12:58:02 2024
Kernel:  AVX+FMA - 1 threads (128 CPU cores detected)

HINT: Use -nt option to specify number of threads because your CPU has 128 cores!
HINT: -nt AUTO will automatically determine the best number of threads to use.

Reading alignment file results_ac2/roary/core_gene_alignment.aln ... Fasta format detected
Reading fasta file: done in 0.100144 secs using 95.78% CPU
Alignment most likely contains DNA/RNA sequences
Constructing alignment: done in 0.505062 secs using 31.51% CPU
Alignment has 3 sequences with 1882304 columns, 105 distinct patterns
0 parsimony-informative, 10079 singleton sites, 1872225 constant sites
       Gap/Ambiguity  Composition  p-value
Analyzing sequences: done in 7.54371e-06 secs
   1  116_2    0.65%    passed     98.28%
   2  E8202    0.84%    passed     99.89%
   3  SRR24    0.78%    passed     99.22%
****  TOTAL    0.76%  0 sequences failed composition chi2 test (p-value<5%; df=3)
NOTE: minimal branch length is reduced to 0.000000053126 for long alignment
ERROR: It makes no sense to perform bootstrap with less than 4 sequences.

@cmkobel
Copy link
Owner Author

cmkobel commented Mar 28, 2024

A possible solution would be to make a checkpoint that checks if there is enough samples, but it seems a bit overengineered.

@cmkobel cmkobel added the wontfix This will not be worked on label May 14, 2024
@cmkobel cmkobel closed this as completed May 14, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
wontfix This will not be worked on
Projects
None yet
Development

No branches or pull requests

1 participant