if ( !stats->is_sorted ) error("The BAM must be sorted in order for -t to work.\n");
regions_t *reg = &stats->regions[bam_line->core.tid];
if ( !stats->is_sorted ) error("The BAM must be sorted in order for -t to work.\n");
regions_t *reg = &stats->regions[bam_line->core.tid];