[EMBOSS] Memory problem with extractseq
michael watson (IAH-C)
michael.watson at bbsrc.ac.uk
Thu Mar 18 09:11:59 UTC 2010
Hi
I'm using EMBOSS 6.1.0 on a fairly small Linux VM which has about 3Gb of RAM.
I find it strange that extractseq reports a memory problem:
-bash-3.00# /usr/local/EMBOSS-6.1.0/bin/extractseq -sequence chr1.fasta -outseq chr1_.1.fasta -regions '34415690-34415711'
Extract regions from a sequence
Uncaught exception: Allocation failed, insufficient memory available, raised at ajstr.c:2406
Whereas if I write a Bioperl script using SeqIO and the trunk() function, it works perfectly.
I'd have thought EMBOSS would be more streamlined and memory efficient than Bioperl?
Thanks
Mick
More information about the EMBOSS
mailing list