Hi everybody, I have 12 fasta files of DNA sequences that I downloaded from NCBI. Every file contains ~ 30000 rows and ~2MLN of characters.
I have Ubuntu 18.04 VPS with 32 GB of RAM, 8 CPUs.
I want multiple alignment of sequences with "msaClustalOmega" function of the Bioconductor package msa. I executed the following code:
mySeqs <- readDNAStringSet(file.choose())
When I executed the istruction
align<-msaClustalOmega(mySeqs, dealign=FALSE) the R process is killed due to lack of RAM.
How to resolve the RAM problem?
How to obtain non-aligned sequences as output? I tried
dealign=TRUE, but I got the same output.
Thank you so much and sorry for my bad English