diff --git a/ketrface/ketrface/dbscan.py b/ketrface/ketrface/dbscan.py
index a7f1a30..8b103f7 100644
--- a/ketrface/ketrface/dbscan.py
+++ b/ketrface/ketrface/dbscan.py
@@ -54,7 +54,7 @@ def DBSCAN(points, eps = MAX_DISTANCE, minPts = MIN_PTS, verbose = True):
 
     sub_perc = -1
     sub_total = len(S)
-    T = S
+    T = S.copy()
     for j, Q in enumerate(S):                       # Process every seed point
       if verbose == True:
         sub_new_perc = int(100 * (j+1) / sub_total)
@@ -73,7 +73,7 @@ def DBSCAN(points, eps = MAX_DISTANCE, minPts = MIN_PTS, verbose = True):
       N = RangeQuery(points, Q, eps)  # Find neighbors
       if len(N) >= minPts:            # Density check (if Q is a core point)
         T = Union(T, N)               # Add new neighbors to seed set
- 
+    S = T
   return clusters
         
 def RangeQuery(points, Q, eps):