[BioC] 'Rsamtools' produces segfaults with XS tag

Misha Kapushesky ostolop at ebi.ac.uk
Fri Aug 13 18:36:44 CEST 2010


Hi,

I can add more detail to this - the same problem:

> library( "Rsamtools" )
...
> tags <- c("NM", "MD", "XA", "XS","NS")
> and_file = "SRR034792/tophat_out/accepted_hits.sortednames.bam"
> and_bam = scanBam(and_file, param = ScanBamParam(flag = scanBamFlag(isUnmappedQuery = FALSE), tag = tags))

*** caught segfault ***
address 0x2b3449beb8cd, cause 'invalid permissions'

Traceback:
1: .Call(func, file, index, "rb", NULL, flag, simpleCigar, ...)
2: .io_bam(.scan_bam, file, index, reverseComplement, tmpl, param = param)
3: scanBam(and_file, param = ScanBamParam(flag = scanBamFlag(isUnmappedQuery = FALSE),     tag = tags))
4: scanBam(and_file, param = ScanBamParam(flag = scanBamFlag(isUnmappedQuery = FALSE),     tag = tags))

This error goes away if we don't include "XS" in tags, also if we use 
version 1.0.0 of Rsamtools. Any ideas?

The sessionInfo() output is below. Thanks in advance!

--Misha

> sessionInfo()
R version 2.11.1 (2010-05-31)
x86_64-unknown-linux-gnu

locale:
[1] C

attached base packages:
[1] tools     stats     graphics  grDevices utils     datasets  methods
[8] base

other attached packages:
[1] Rsamtools_1.1.11    Biostrings_2.16.9   GenomicRanges_1.0.7
[4] IRanges_1.6.11      svMisc_0.9-57       JavaGD_0.5-2
[7] rJava_0.8-4

loaded via a namespace (and not attached):
[1] Biobase_2.8.0


On Thu, 12 Aug 2010, Andrew Tikhonov wrote:

> Dear Biocore Team,
>
> 	We've got a problem where 'Rsamtools' produces segfaults and 
> core-dumps when reading a .bam file using the "XS" tag. We have a faulty .bam 
> and the list or R commands to reproduce it. Could you please point us to 
> someone who we may discuss this matter with and, if possible, fix the 
> problem.
>
> 	Thanks a lot in advance.
>
> Regards,
> Andrew
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: 
> http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list