Limit of sequences for SMS - Smart model selection
Entering edit mode
2.2 years ago
flogin ▴ 280

Hello guys, I'm testing the SMS to select best substitution models on alignment (, but, due the large number of sequences in my alignments (6413) I'm testing in a local server.

My problem is that this tool appear accept only datasets with 4000 or less sequences, and when I tried to run with my dataset the error of 4k sequences occurred.

Starting SMS v1.8.4
Input alignment: teste.phy
Data type: DNA
Number of taxa: 6413
Number of sites: 228 Number of branches : 12823
Criterion : AIC

Step 1 : Set a fixed topology Unexpected error while computing Log-likelihood

The number of taxa should not exceed 4000phyml: io.c:1293: Get_Seq_Phylip: Assertion `FALSE' failed.

Someone here know the best way to avoid this situation? (with SMS or with other tool to select the best model of nucleotide or aminoacids in large scale?).


Ps: I run this tool with a subset of 3k sequences, and everything runs ok...

alignment sms • 426 views
Entering edit mode
2.2 years ago
flogin ▴ 280

After an e-mail to author of SMS, this is how can I've solved this problem:

I download the source code of phyml on and I modify the value of variable N_MAX_OTU from /src/utilities.h. After this, I compiled the phyml, and copy all files of phyml directory to sms/phyml-src/phyml-3.3.20190909/.

So these steps allow the use of SMS with alignments with more than 4,000 sequences....


Login before adding your answer.

Traffic: 1425 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6