[BioC] narrow on GappedAlignments fails with error

Hubert Rehrauer Hubert.Rehrauer at fgcz.ethz.ch
Fri Jan 20 14:14:17 CET 2012


Hi

There seems to be an overvflow problem if the GappedAlignments object 
holds many reads:

With foo being a gapped alignments having 17836148 reads I get:
>  bar = narrow( foo[1:10000000],  di[1:10000000]) ## fine
> bar = narrow( foo[10000000:17836148],  di[10000000:17836148]) ## fine
> bar = narrow( foo,  di)
> Error during wrapup: in 'cigar' element 11380048: CIGAR is empty after 
> narrowing
> Warning messages:
> 1: In successiveIRanges(cigarToWidth(cigar)) :
>   Integer overflow in 'cumsum'; use 'cumsum(as.numeric(.))'
> 2: In start(x) + width(x) : NAs produced by integer overflow
> 3: In start(x) + width(x) : NAs produced by integer overflow
> 4: In start(x) + width(x) : NAs produced by integer overflow
can you help???

hubert



More information about the Bioconductor mailing list