def get_lineage_size(df_raw, lineage_uid):
    """ Get number of sequences in a lineage """
    df = df_raw[df_raw["lineage_uid"] == lineage_uid]
    N = df.shape[0] # number of sequences in lineage
    return N