def get_called_variants(var_list, cn_prob_processed, starting_index=0):
    """
    Return called variants based on called copy number and list of variant names
    """
    total_callset = []
    if starting_index != 0:
        assert len(var_list) == len(cn_prob_processed) + starting_index
    for i, cn_called in enumerate(cn_prob_processed):
        if cn_called is not None and cn_called != 0:
            for _ in range(cn_called):
                total_callset.append(var_list[i + starting_index])
    return total_callset