def get_sample_region_coverage(sample_tabix_object, region_dict):
    """
    get_sample_region_coverage
    ==========================
    Get the coverage values for a specific sample at a specific region interval.

    A list containing per base coverage values will be returned.

    Parameters:
    -----------
    1) sample_tabix_object: (psyam Object) A pysam tabix object representing the bed file of a sample
    2) region_dict: (dict) A dictionary with region information. Required keys include "chrom", "start", and "end"

    Returns:
    +++++++
    1) (list) A list of coverage values representing a per base coverage score
    """

    coverage = []
    ## Iterate over the overlapping coverage intervals
    for i, coverage_line in enumerate(
        sample_tabix_object.fetch(
            region_dict["chrom"], int(region_dict["start"]), int(region_dict["end"])
        )
    ):

        ## Create a dict for the current line
        cov_dict = dict(
            zip(
                ["chrom", "start", "end", "coverage_value"],
                coverage_line.strip().split("\t"),
            )
        )

        ## Add the coverage value for each position defined by the current interval
        coverage.extend(
            [
                int(cov_dict["coverage_value"])
                for x in range(int(cov_dict["end"]) - int(cov_dict["start"]))
            ]
        )

    return coverage