Ketil Malde <[email protected]> writes:

> I've run into something...troublesome when using the samtools library.
> I've written a function, left_links, that takes an IdxHandle, reads all
> matches in the region (0,350), and filters out just the ones that point
> "left".  Or so I thought.

My own fault, really.  My guess is that IdxHandle contains a bunch of
state, and that it doesn't really go down well with lazy reading from a
region of the Bam file.  It /seems/ that strictly evaluating the list of
Bam1's solves the problem.

-k
-- 
If I haven't seen further, it is by standing in the footprints of giants
_______________________________________________
Biohaskell mailing list
[email protected]
http://malde.org/cgi-bin/mailman/listinfo/biohaskell

Reply via email to