Thank you for the suggestions. I will consider changing my strategy.
• 0 views
•
link
written
by
atitparajuli2018 •
Multiple Sequence Alignment using muscle
Hi,
I keep getting an error segmentation fault 11 while running muscle command to align multiple sequences from a fasta file containing 112626 ESTs. Is it because of the large number of sequences that is causing the problem?
I used the command
$ ./muscle -in input.fa -out output.afa -maxiters 1 -diags1 -sv
• 3,399 views
•
link
updated
by
h.mon
•
written
by
atitparajuli2018 •
•
written
by
atitparajuli2018 •
1 answer
You have been told repeatedly your approach is flawed, but you insist in following it. Please reconsider, and try to take a different approach: there are several papers describing SNP prospection from ESTs, why don't you follow one of them? The general workflow is to assemble the ESTs (generally with CAP3), then map the ESTs to the assembly and detect polymorphisms. Here are some links with references:
Log in to answer this question.
More posts like this
-
Segmentation fault (core dumped) error when running MUSCLE
written by Zhang •
-
Segmentation fault while running muscle
written by alicecol •
-
Muscle alignment error
written by hypeanut •
-
Samtools faidx not printing reverse complement
written by jamie.pike
-
Multiple Sequence Alignment using muscle
written by atitparajuli2018 •
-
Error characters in fasta file
written by atitparajuli2018 •
-
Trouble running Muscle on Mac OS
written by ricfoz •
-
Using Biopython to access MUSCLE
written by jtbioinfo •
-
muscle v3.8.31 large sequence Segmentation fault
written by singh.punit1990 •
-
RNAstar alignment segmentation fault
written by sam •
Probably. Why in the world would you like to MSA 112626 ESTs?
Most likely yes.
Are you trying to deduplicate this dataset? If that is the case then there are programs like CD-HIT which are more appropriate.
I am preparing an exome capture library from these 112K ESTs.
Can you clarify what that means? What are you hoping to do by aligning 112K sequences? Are they all related to each other? Otherwise there is no point in trying to do MSA with them.
I have collected 320k ests from NCBI database, 47209 Ests from DDBJ and 54K Ests from Alfalfa Genome index (AGED). I am aligning these Ests together to generate a sequence which I will use as a reference sequence for Identifying the SNPs to the alfalfa lines that I have. I am working on identifying genetic variants associated with self incompatibility in alfalfa.
You are unlikely to get the results you want by following the strategy you describe here. I suggest taking a step back and thinking about why doing an MSA is not going to work in this case. Think about how you could use CD-HIT to get started.
So you design capturing beads? What does a MSA help? There is software provided by the companies that make the beads for this task. Better use them instead of custom approaches that are completely untested.