On Wed, 2015-05-13 at 21:03 -0600, Brent Pedersen wrote:
> I finally got around to this. For now, for simplicity, I went with
> option 3: just Read() the requested chunk into a buffer.
>
https://github.com/brentp/faidx
> I'm happy to contribute this to biogo if you think it's appropriate.
>
I don't think this gets us anything that current biogo seq io doesn't
already provide. If you implement the mmap approach then it would be a
more useful addition.
The code that you have here makes a basis for that though.
Before it would be merged, indexing would need to be changed to 0-based
half open and I think providing an faidx association rather than using
an implied file extension would be better (I think
base.{bai,csi,fai,tbi} is a horrible pattern and should not ever work
it's way into library routines - app-level code may want to do that
though, maybe).
Dan