I just tried to run cmaple on some 5000 mpox sequences. I got this error and don't really understand what it means.
- What is an "invalid mutation"
- How can a mutation have a "length"?
- 2^15 is around 32k, so is there a sequence that has 32k SNPs vs reference? Or more likely at least 32k non-reference - it could be Ns for example?
If it's a single sequence that's causing this issue, maybe it might be worth reporting it to the user so they can take action. Recompiling isn't really something that most users would do, so better to explain how users can fix their data.
cmaple -s /Users/cr/code/mpox/nextclade/results/prealigned.fasta
CMAPLE version 1.1.0 for Mac OS X 64-bit built Dec 10 2024
Developed by Nhan Ly-Trong, Chris Bielow, Nicola De Maio, Bui Quang Minh.
Command: cmaple -s /Users/cr/code/mpox/nextclade/results/prealigned.fasta
Seed: 238816
Reading an alignment
ERROR: Invalid mutation. Length is larger than 2^15. Recompile with larger 'LengthType' at the cost of higher memory consumption.
prealigned.fasta.zst.txt