diff --git a/Histone_TFSEE/cluster1_enriched_tfs.png b/Histone_TFSEE/cluster1_enriched_tfs.png
deleted file mode 100644
index ef6c247344729c52152cd4e91ae5edba694aea69..0000000000000000000000000000000000000000
Binary files a/Histone_TFSEE/cluster1_enriched_tfs.png and /dev/null differ
diff --git a/Histone_TFSEE/cluster2_enriched_tfs.png b/Histone_TFSEE/cluster2_enriched_tfs.png
new file mode 100644
index 0000000000000000000000000000000000000000..1aded192feb99d33515b97068ee5669fa791abd2
Binary files /dev/null and b/Histone_TFSEE/cluster2_enriched_tfs.png differ
diff --git a/Histone_TFSEE/cluster3_enriched_tfs.png b/Histone_TFSEE/cluster3_enriched_tfs.png
index acc812e4c3471d743c4c25a5e71157be6a7b3634..31178fdfdfa3d33829fc47eb15f4684aca482e2c 100644
Binary files a/Histone_TFSEE/cluster3_enriched_tfs.png and b/Histone_TFSEE/cluster3_enriched_tfs.png differ
diff --git a/Histone_TFSEE/rank_order.py b/Histone_TFSEE/rank_order.py
index 2e032f4d0784cf873d9b786d517b2f036cb4b20e..fe840bcc0435335db9599b8dbeba79c62d2e6c43 100644
--- a/Histone_TFSEE/rank_order.py
+++ b/Histone_TFSEE/rank_order.py
@@ -8,30 +8,30 @@ import matplotlib.pyplot as plt
 
 tfsee = pd.read_csv('clustering_tfs.csv')
 
-tfsee_cluster1 = tfsee[tfsee['cluster'] == 0]
+tfsee_cluster2 = tfsee[tfsee['cluster'] == 1]
 
 
-tfsee_cluster1['early'] = tfsee_cluster1[['ES_D0','ES_D2','ES_D5']].mean(axis=1)
-tfsee_cluster1['late'] = tfsee_cluster1[['ES_D7','ES_D10']].mean(axis=1)
-tfsee_cluster1['diff'] = tfsee_cluster1['early'] - tfsee_cluster1['late']
-tfsee_cluster1['rank'] = tfsee_cluster1['diff'].rank()
+tfsee_cluster2['early'] = tfsee_cluster2[['ES_D0','ES_D2','ES_D5']].mean(axis=1)
+tfsee_cluster2['late'] = tfsee_cluster2[['ES_D7','ES_D10']].mean(axis=1)
+tfsee_cluster2['diff'] = tfsee_cluster2['early'] - tfsee_cluster2['late']
+tfsee_cluster2['rank'] = tfsee_cluster2['diff'].rank()
 
-x = list(tfsee_cluster1['rank'])
+x = list(tfsee_cluster2['rank'])
 
-z = np.polyfit(tfsee_cluster1['rank'], tfsee_cluster1['diff'], 3)
+z = np.polyfit(tfsee_cluster2['rank'], tfsee_cluster2['diff'], 3)
 f = np.poly1d(z)
-x_new = np.linspace(1, 24, num=len(x)*10)
+x_new = np.linspace(1, 19, num=len(x)*10)
 
 
-plt.figure(figsize=(25,20))316A45
+plt.figure(figsize=(25,20))
 plt.plot(x_new, f(x_new), color = 'k', linewidth=4.0)
-plt.scatter(x=tfsee_cluster1['rank'], y=tfsee_cluster1['diff'], color='#316A45', s=600)
+plt.scatter(x=tfsee_cluster2['rank'], y=tfsee_cluster2['diff'], color='#316A45', s=600)
 plt.ylim([-1.5,4.5])
-plt.xlim([0,25])
+plt.xlim([0,20])
 plt.suptitle('Early Enriched TFS', fontsize=8, fontweight='bold')
 plt.xlabel('Rank Order',fontsize=8, fontweight='bold')
 plt.ylabel('delta z',fontsize=8, fontweight='bold')
-plt.savefig('cluster1_enriched_tfs.png')
+plt.savefig('cluster2_enriched_tfs.png')
 plt.clf()
 
 
@@ -48,14 +48,14 @@ x = list(tfsee_cluster3['rank'])
 
 z = np.polyfit(tfsee_cluster3['rank'], tfsee_cluster3['diff'], 3)
 f = np.poly1d(z)
-x_new = np.linspace(1, 17, num=len(x)*10)
+x_new = np.linspace(1, 24, num=len(x)*10)
 
 
 plt.figure(figsize=(25,20))
 plt.plot(x_new, f(x_new), color = 'k', linewidth=4.0)
 plt.scatter(x=tfsee_cluster3['rank'], y=tfsee_cluster3['diff'], color='#3999CA', s=600)
 plt.ylim([-1.5,4.5])
-plt.xlim([0,20])
+plt.xlim([0,25])
 plt.suptitle('Late Enriched TFS', fontsize=8, fontweight='bold')
 plt.xlabel('Rank Order',fontsize=8, fontweight='bold')
 plt.ylabel('delta z',fontsize=8, fontweight='bold')