def METHOD_NAME(self, somatic_variant, alt, transcript):
    (phased_somatic_variant, potential_proximal_variants) = self.find_phased_somatic_variant_and_potential_proximal_variants(somatic_variant, alt, transcript)
    if phased_somatic_variant is None:
        print("Warning: Main somatic variant not found in phased variants file: {}, {}".format(somatic_variant, alt))
        return []
    if len(potential_proximal_variants) == 0:
        return []
    proximal_variants = []
    sample = self.proximal_variants_vcf.header.samples.names[0]
    phased_somatic_variant_genotype = phased_somatic_variant.call_for_sample[sample]
    if 'HP' in phased_somatic_variant.FORMAT:
        somatic_phasing = phased_somatic_variant_genotype.data['HP']
        for (entry, csq_entry, proximal_alt) in potential_proximal_variants:
            proximal_variant_genotype = entry.call_for_sample[sample]
            #identify variants that are in phase with the phased_somatic_variant
            if 'HP' in entry.FORMAT:
                proximal_variant_phasing = proximal_variant_genotype.data['HP']
                if proximal_variant_phasing == somatic_phasing:
                    proximal_variants.append([entry, csq_entry, proximal_alt])
            #proximal variant is hom var
            elif proximal_variant_genotype.is_variant and not proximal_variant_genotype.is_het:
                #main somatic variant is het
                if phased_somatic_variant_genotype.is_het:
                    proximal_variants.append([entry, csq_entry, proximal_alt])
                #main somatic variant is hom var
                if phased_somatic_variant_genotype.is_variant and not phased_somatic_variant_genotype.is_het:
                    proximal_variants.append([entry, csq_entry, proximal_alt])
    else:
        for (entry, csq_entry, proximal_alt) in potential_proximal_variants:
            proximal_variant_genotype = entry.call_for_sample[sample]
            #proximal variant is hom var
            if proximal_variant_genotype.is_variant and not proximal_variant_genotype.is_het:
                #main somatic variant is het
                if phased_somatic_variant_genotype.is_het:
                    proximal_variants.append([entry, csq_entry, proximal_alt])
                #main somatic variant is hom var
                if phased_somatic_variant_genotype.is_variant and not phased_somatic_variant_genotype.is_het:
                    proximal_variants.append([entry, csq_entry, proximal_alt])
    return proximal_variants