just to follow up in case anyone else wonders about this... it seems to work fine to do essentially as described to allow indexing while writing.
On Tue, Sep 19, 2017 at 9:55 AM, Brent Pedersen <bpede...@gmail.com> wrote: > based on the API: > int hts_idx_push(hts_idx_t *idx, int tid, int beg, int end, > uint64_t offset, int is_mapped); > and a quick look at the code, it looks like I might be able to create > a .csi index while writing a bgzip output. Is that correct? > > basically, if I've just written a line/interval to fp and then do: > > hts_idx_push(idx, tid, start, stop, bgzf_tell(fp)) > > will the bgzf_tell give the position that hts_idx_push expects? or is > there some issue with waiting til the block is complete? > > thanks, > -Brent ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Samtools-help mailing list Samtools-help@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/samtools-help
