def METHOD_NAME(self, contig, pos):
    """
    Given a contig name and a position, estimate the copy number of that
    position in the genome.
    Return the number of instances expected (1 for non-repetitive sequence).
    """
    # TODO: use contig
    # Get the set of overlapping things
    overlaps = self.tree[pos]
    # Keep track of the number of copies of the most numerous repeat
    # observed.
    max_copies = 1
    for interval in overlaps:
        # For each repeat we are in
        # Max in how many copies of it exist
        max_copies = max(max_copies, self.counts[interval.data])
    return max_copies