muscle icon indicating copy to clipboard operation
muscle copied to clipboard

segmentation fault muscle -align testing.fa -output alignment.afa

Open m-leitman opened this issue 1 year ago • 1 comments

When I run the align command for two sequences, I get this error:

muscle 5.1.osxarm64 [ef62d8-dirty] 17.2Gb RAM, 8 cores Built Jan 15 2022 11:59:50 (C) Copyright 2004-2021 Robert C. Edgar. https://drive5.com

Input: 2 seqs, length avg 32604 max 33281

WARNING: Sequence length >5k may require excessive memory

00:00 12Mb 100.0% Derep 1 uniques, 0 dupes 00:00 14Mb 100.0% UCLUST 2 seqs EE<0.01, 1 centroids, 0 members 00:00 14Mb CPU has 8 cores, running 8 threads
zsh: segmentation fault muscle -super5 testing.fa -output alignment.afa

Does my device not have enough RAM to do this?

m-leitman avatar Apr 21 '24 06:04 m-leitman

Yes, as the warning message says.

rcedgar avatar Apr 21 '24 13:04 rcedgar

I am using muscle 3 and always get this error message. Tried it with only 17 protein sequences in fasta format

filtered_sequences_WP_only 17 seqs, max length 13477, avg length 2937 00:00:00 2 MB(0%) Iter 1 100.00% K-mer dist pass 1 00:00:00 2 MB(0%) Iter 1 100.00% K-mer dist pass 2 00:00:10 168 MB(1%) Iter 1 100.00% Align node 00:00:10 168 MB(1%) Iter 1 100.00% Root alignment 00:00:18 181 MB(1%) Iter 2 100.00% Refine tree 00:00:18 182 MB(1%) Iter 2 100.00% Root alignment 00:00:18 182 MB(1%) Iter 2 100.00% Root alignment zsh: segmentation fault muscle3 -in -out Refine biparts

Is it because the variation in the length of the protein sequences too large?

chingciripit avatar Jul 10 '24 21:07 chingciripit

13477 is probably too long.

rcedgar avatar Jul 12 '24 23:07 rcedgar