def get_gene_capacity(genes_df, database, col='GeneID'):
    """Get gene capcaity from the stored metadata"""
    capacity = (database.groupby('geneid').capacity.mean()
                .to_frame(name='GeneCapacity'))
    genes_df = genes_df.merge(capacity, how='left', left_on='GeneID', right_index=True)
    return genes_df