def extract_bcf_stats(path):
    """
    Extract relevant information from BCF Stats file for a single sample

    Specifically extract SNPs, MNPs, indels, "others", multiallelic sites,
    and multiallelic SNPsites.

    No effort is made to convert strings to numbers for the stat values.

    Args:
        path (str): path to the BCF stats file

    Returns:
        dict: keys as stat names and the values as stat values

    """

    # Not ideal to hardcode here nor below, but gets the job done
    stats_of_interest = {"number of SNPs:",
                         "number of MNPs:",
                         "number of indels:",
                         "number of others:",
                         "number of multiallelic sites:",
                         "number of multiallelic SNP sites:"}

    stats = {}

    with open(path) as statsf:
        for line in statsf:

            if line.startswith("SN"):
                parts = line.strip().split("\t")
                stat = parts[-2]
                num = parts[-1]

                if stat in stats_of_interest:
                    stats[stat.strip(":")] = num

    return stats