Question: (Closed) Multiple Sequence Alignment using muscle
gravatar for atitparajuli2018
7 months ago by
atitparajuli20180 wrote:

I tried to align 112626 sequences of alfalfa using muscle command, but I get segmentation fault 11 error which is mainly due to low running memory. Then I split into 8 smaller files with 14100 sequences in each file. The muscle command ran smoothly in those files and generated output. However, when I used the -profile option in muscle command to join those 8 aligned files (2 at a time) it gives me error message (* WARNING Invalid byte hex 00 in FASTA sequence data, ignored ERROR * Internal error MSA::ExpandCache, ColCount changed). Why am I getting this error message? and how can I join those 8 small files into single aligned file?

ADD COMMENTlink written 7 months ago by atitparajuli20180

It seems you want to align a transcriptome or an EST set using muscle, and this makes no sense. Muscle is designed to align sequences with moderate or high similarity (either because they are homologous, or share structural features). What is your goal? Maybe you want first to group the 112626 sequences into clusters with a certain level of similarity, then align these clusters?

ADD REPLYlink written 7 months ago by h.mon29k

While you posted a specific error in this new thread, we have tried to help you with this line of analysis in a past thread: Multiple Sequence Alignment using muscle

What you are trying to do does not make scientific sense.

ADD REPLYlink modified 7 months ago • written 7 months ago by genomax78k

Hello atitparajuli2018!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.


ADD REPLYlink written 7 months ago by h.mon29k
Please log in to add an answer.
The thread is closed. No new answers may be added.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1597 users visited in the last hour