diff --git a/scripts/functions/stability_count_bp.R b/scripts/functions/stability_count_bp.R
index 8f7d9ed..e5ed684 100644
--- a/scripts/functions/stability_count_bp.R
+++ b/scripts/functions/stability_count_bp.R
@@ -15,39 +15,42 @@ theme_set(theme_grey())
   ## ...opt args
 #==========================================================
 stability_count_bp <- function(plotdf
-         , df_colname
-         , leg_title = "Legend title"
-         , axis_text_size = 25
-         , axis_label_size = 22
-         , leg_text_size = 20
-         , leg_title_size = 22
+         , df_colname = ""
+         , leg_title = "Legend Title"
+         , ats = 25     # axis text size
+         , als = 22     # axis label size
+         , lts = 20     # legend text size
+         , ltis = 22    # label title size
+         , geom_ls = 10 # geom_label size
          , yaxis_title = "Number of nsSNPs"
          , bp_plot_title = ""
          , label_categories = c("Destabilising", "Stabilising")
          , title_colour = "chocolate4"
          , subtitle_text = NULL
-         , subtitle_size = 20
+         , sts = 20
          , subtitle_colour = "pink"
          #, leg_position = c(0.73,0.8) # within plot area
          , leg_position = "top"){
  
-  OutPlot_count = ggplot(plotdf, aes(x = eval(parse(text = df_colname)))) + 
+#  OutPlot_count = ggplot(plotdf, aes(x = eval(parse(text = df_colname)))) +
+  OutPlot_count = ggplot(plotdf, aes_string(x = df_colname)) +     
     geom_bar(aes(fill = eval(parse(text = df_colname))), show.legend = TRUE) +
     geom_label(stat = "count"
                , aes(label = ..count..)
                , color = "black"
                , show.legend = FALSE
-               , size = 10) +
+               , size = geom_ls) +
     theme(axis.text.x = element_blank()
           , axis.title.x = element_blank()
-          , axis.title.y = element_text(size =  axis_label_size)
-          , axis.text.y = element_text(size = axis_text_size)
+          , axis.title.y = element_text(size =  als)
+          , axis.text.y = element_text(size = ats)
           , legend.position = leg_position
-          , legend.text = element_text(size = leg_text_size)
-          , legend.title = element_text(size =  leg_title_size)
-          , plot.title = element_text(size =  axis_label_size
-                                      , colour = title_colour)
-          , plot.subtitle = element_text(size = subtitle_size
+          , legend.text = element_text(size = lts)
+          , legend.title = element_text(size =  ltis)
+          , plot.title = element_text(size =  als
+                                      , colour = title_colour
+                                      , hjust = 0.5)
+          , plot.subtitle = element_text(size = sts
                                          , hjust = 0.5
                                          , colour = subtitle_colour)) + 
     labs(title      = bp_plot_title