[BioC] Restriction enzyme digestion

[hpages at fhcrc.org]

Quoting hpages at fhcrc.org:
[...]
> You could also use the getSeq() function in a vectorized way if you
> really wanted to use getSeq(). Then your code would look like:
>
>   reslt_start <- end(m)[-length(m)] + 1L
>   reslt_end <- start(m)[-1] - 1L
>   reslt <- getSeq(Hsapiens, 'chr1', start(m)[i]+2, start(m)[i+1]-2))

I meant this:

   reslt_start <- end(m)[-length(m)] + 1L
   reslt_end <- start(m)[-1] - 1L
   reslt <- getSeq(Hsapiens, 'chr1', start=reslt_start, end=reslt_end)

sorry... (but as I said in my previous post, getSeq() needs some improvements
and fixes)

Cheers,
H.