Question: seqfault wiht msa
1
gravatar for jmh5ad
16 months ago by
jmh5ad20
jmh5ad20 wrote:

I am working with one of our researchers who is doing multiple sequence alignment. 

She has 200,000+ sequences in a file. Her code is basically the following:

library("msa")
mySeq <- readAAStringSet("myData.fa")

myFirstAlignment <- msa(mySeq)
msaPrettyPrint(MyFirstAlignment, output=c("pdf", "tex", "dvi", "asis"))

But, R is crashing during the msa function, with the following error message:

 *** caught segfault ***
address 0x7f5ace31d0d0, cause 'memory not mapped'

Traceback:
 1: .Call("RClustalW", inputSeqs, cluster, abs(gapOpening), abs(gapExtension),     maxiters, substitutionMatrix, type, verbose, params, PACKAGE = "msa")

 

We are running R 3.4.3.  I have uninstalled and reinstalled the msa package and its dependencies.

Is there anything else we could try, or is this a limitation due to the size of the data?

 

 

 

 

 

 

 

 

segfault msa • 252 views
ADD COMMENTlink modified 15 months ago • written 16 months ago by jmh5ad20
Answer: seqfault wiht msa
1
gravatar for jmh5ad
15 months ago by
jmh5ad20
jmh5ad20 wrote:

Thank you for your response.  We will try ClustalOmega.

 

 

ADD COMMENTlink written 15 months ago by jmh5ad20
Answer: seqfault wiht msa
0
gravatar for UBodenhofer
16 months ago by
UBodenhofer250
Johannes Kepler University, Linz, Austria
UBodenhofer250 wrote:

Thanks for reporting this issue! As far as I can see, my suspicion is indeed the vast amount of data. Actually, ClustalW is not made for such excessive amounts of data. I would rather suggest to use ClustalOmega instead. If that does not work via 'msa' either, maybe the last resort is to use stand-alone ClustalOmega. Sorry that I cannot be more helpful.

ADD COMMENTlink written 16 months ago by UBodenhofer250
Please log in to add an answer.

Help
Access

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