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

Reply via email to