Function call_nucleosomes 

pub fn call_nucleosomes(
    scores: &[(usize, f64)],
    record_len: usize,
    params: &CallParams,
) -> Vec<NucleosomeCall>

Turn symmetry scores into overlapping nucleosome calls.

scores are (dyad, score) pairs with zero-based dyads, in ascending order. A dyad is called when its score is above zero and its whole footprint fits inside the record.

Following the reference, the bounds are strict on both sides: dyad - half_width must be greater than zero, not merely non-negative, and dyad + half_width must be less than the record length rather than within it. That drops one otherwise-callable dyad at each end, and is reproduced so the positions match.