1

I am running the msa package to create a DNA alignment for the phangorn package and it crashes with this error

I am running RStudio with R v4.3.1 on an M1 Mac Book Pro

mult <- msa(seqs, method="Muscle", type="dna", order="input")

That results in:

*** ERROR ***  MSA::SetIdCount: cannot increase count

Fatal error, exception caught.
Error in msaFun(inputSeqs = inputSeqs, cluster = cluster, gapOpening = gapOpening,  : 
  MUSCLE finished by an unknown reason

The input file is 6440 DNA seqs with an average length 1500 bp.

When I run Muscle in a stand alone mode on a linux virtual server I have to use the super5 option, because it gives this warning

WARNING: >1k sequences, may be slow or use excessive memory, consider using -super5

Floating point exception (core dumped)
(base) ubuntu@jrmicl-clovr:~/muscle$ muscle5 -super5 seqs.fas -output seqs_aln.fas

Any ideas - is this a memory thing?

Krokomot
  • 3,208
  • 2
  • 4
  • 20

0 Answers0