diff --git a/config.yaml b/config.yaml
index 4bfe9cfc6ec68aa4fe9e522d4f3c077cb5ee03e8..a165debec251c77e18fa8d54fd233333745de79f 100644
--- a/config.yaml
+++ b/config.yaml
@@ -48,3 +48,10 @@ preprocess:
     # http://seqanswers.com/forums/archive/index.php/t-42552.html
     bbmap_reference: '/home/mticlla/Documents/Git_repositories/metagenomicsnake/ref/hg19_main_mask_ribo_animal_allplant_allfungus.fa.gz'
     bbmap_ref_dir: '/home/mticlla/Documents/Git_repositories/metagenomicsnake/ref/ref_bbsplit'
+    
+phlanprof:
+    strphlan_clade_profiling:
+        metadatas:
+            - 'SubjectID'
+            - 'Group'
+        colorized_metadata: 'Group'
diff --git a/rules/phlanprof.smk b/rules/phlanprof.smk
index e95430237902b6ed7a6200d35777b16ac3dbd0d2..ebad71b027e189649cc948811adbba120151b836 100644
--- a/rules/phlanprof.smk
+++ b/rules/phlanprof.smk
@@ -220,7 +220,9 @@ rule strphlan_clade_profiling:
         clade_results = OUT_DIR + '/{dataset}/PhlAnProf/logs/strphlan/clades_profiled/{clade}.log'
     params:
         mpa_pkl = METASNK_DB_DIR + '/metaphlan_databases/mpa_v20_m200.pkl',
-        strpa_dir = METASNK_DB_DIR+'/strainphlan_markers'
+        strpa_dir = METASNK_DB_DIR+'/strainphlan_markers',
+        metadatas = config['phlanprof']['strphlan_clade_profiling']['metadatas'],
+        colorized_metadata = config['phlanprof']['strphlan_clade_profiling']['colorized_metadata']
     threads: cpus_avail
     singularity:singularity_metaprof
     group: 'strphlan_clade_analysis'
@@ -239,12 +241,20 @@ rule strphlan_clade_profiling:
         add_metadata_tree.py \
         --ifn_trees {output.clade_result_dir}/RAxML_bestTree.{wildcards.clade}.tree \
         --ifn_metadatas {input.samples_metadata} \
-        --metadatas SubjectID; \
+        --metadatas {params.metadatas}; \
         plot_tree_graphlan.py \
         --ifn_tree {output.clade_result_dir}/RAxML_bestTree.{wildcards.clade}.tree.metadata \
-        --colorized_metadata SubjectID \
+        --colorized_metadata {params.colorized_metadata} \
         --leaf_marker_size 60 \
-        --legend_marker_size 60; fi') &>{log}
+        --legend_marker_size 60; \
+        deactivate; \
+        source activate py27breadcrumbs; \ 
+        strainphlan_ggtree.R \
+        {output.clade_result_dir}/RAxML_bestTree.{wildcards.clade}.tree \
+        {input.samples_metadata} \
+        {output.clade_result_dir}/{wildcards.clade}.fasta \
+        {output.clade_result_dir}/RAxML_bestTree.{wildcards.clade}.tree_1.png \
+        {output.clade_result_dir}/RAxML_bestTree.{wildcards.clade}.tree_2.png; fi') &>{log}
         '''
 
 def aggregate_clade_profiles(wildcards):