def get_vcf_column_headers(line, sample_ids, sample_indices):
    """Get VCF column headers, including only relevant samples
    
    Parameters
    ----------
    line : str
        A line from a vcf file (e.g. from an open file descriptor)
    sample_ids : list, tuple
        All sample IDs in the input VCF file
    sample_indices : container
        Indices indicating which samples from sample_ids should be kept
    
    Returns
    -------
    list
        Column headers for the output VCF file
    """
    
    return line.split()[:9] + [sample_ids[i] for i in sample_indices]