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

You should provide an unzipped VCF file to the 'SNPAnnotationFile' parameter.