joshuak94 / BAMIntervalTree

BSD 3-Clause "New" or "Revised" License
2 stars 1 forks source link

Update bamit to return vector of records. #40

Closed joshuak94 closed 2 years ago

joshuak94 commented 2 years ago

@Irallia The goal is to have a function which returns a "seqan-like" vector of sam_record objects, so that I could do something like this:

auto records_vector = bamit::get_overlap_records(input, index, start, end)
uint32_t avg_read_depth = records_vector.size() / (end - start);

At the moment I was only returning the file position, so the user would have to parse the records themself. This just makes it a bit easier.