Thank you for the reply, Erik. The "count breakpoints in each bin" part doesn't seem to work. This is the error I got:
Error in hist.default(synteny[2, 1][[1]][, c("start2", "end2")], breaks = seq(from = 0, :
some 'x' not counted; maybe 'breaks' do not span range of 'x'
Although R is not my weapon of choice, this package seems useful. Though, I'm worried that you (and the package) say synteny, which is not really what I'm looking for. I'm interested in normal sequence alignments, like what BLAST gives. Finding where the endpoints of the scaffolds map to the reference is not what I'm looking for either, since those aren't necessarily breakpoints.
What I have right now is a script that parses BLAST's output and reports the endpoints of all the alignments, except the endpoints of the entire scaffold. I'm considering these coordinates the "breakpoints".
I appreciate your taking the time to answer so here's an upvote!