rsubread exactsnp ERROR "memory not mapped"
0
0
Entering edit mode
5.1 years ago
tarek.mohamed ▴ 300

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

Rsubread exactSNP • 1.8k views
ADD COMMENT
0
Entering edit mode

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)?

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2510 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6