def ConSpliceML_score(rf, var_df, feature_col_names, new_col_name="ConSpliceML"):

    """
    ConSpliceML_score
    =================
    Score a set of variants using the ConSpliceML Random Forest trained model.

    Parameters:
    -----------
    1) rf:          (sklearn RF) The ConSpliceML trained RF model
    2) var_df:       (Pandas DF) A pandas Dataframe with variant info. Must contain the feature columns used to train the model
    3) feature_col_names: (list) A sorted list of features that match the scored ConSpliceML model to used to score the variant df.
    4) new_col_name:       (str) The name of the new column to add to the dataframe

    Returns:
    ++++++++
    1) (Pandas DataFrame) An updated pandas dataframe with the ConSpliceML score for the variants in the file
    2)              (str) The name of the ConSpliceML column added to the df.
    """

    var_df[new_col_name] = rf.predict_proba(var_df[feature_col_names])[:, 1]
    scored_df = var_df

    return (scored_df, new_col_name)