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