rsubread exactsnp ERROR "memory not mapped"

Question

[tarek.mohamed]

Hi All,

I am trying to call SNPs from RNAseq data. I am using exactSNP function from Rsubread package.

I got an error ( * caught segfault * address (nil), cause 'memory not mapped')

The script is as follows;

library(Rsubread)
exactSNP(SNPAnnotationFile="All_20151104.vcf.gz",readFile="PB01_1_sv.sam",isBAM=FALSE,refGenomeFile="BSgenome.Hsapiens.NCBI.GRCh38.fasta",outputFile="PB01_1_sv.sam.vcf",nthreads=8)

 *** caught segfault ***

address (nil), cause 'memory not mapped'

Traceback:
 1: .C("R_SNPcalling_wrapper", as.integer(n), as.character(cmd),     PACKAGE = "Rsubread")
 2: exactSNP(SNPAnnotationFile = "All_20151104.vcf.gz", readFile = "PB01_1_sv.sam",     isBAM = FALSE, refGenomeFile = "BSgenome.Hsapiens.NCBI.GRCh38.fasta",     outputFile = "PB01_1_sv.sam.vcf", nthreads = 8)
aborting ...

Thanks,

Tarek

[GenoMax]

With memory errors the first thing to suspect is that you don't have enough available. Have you tried to run this with a smaller number of threads (instead of 8 above)?

[tarek.mohamed]

Hi, No I have not, but I am using the university server and I know that these threads are available. Thanks Tarek

[tarek.mohamed]

Hi, I tried it with only 2, and I got the same error

[GenoMax]

Are you allocating a specific amount of memory for the job (is there a default value, especially if you are using a job scheduler)?