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