diff --git a/setup.py b/setup.py
index 6987d69654bc37a341da775544429c4244bb8ee5..6beb1ecef15ecf22e517954958dfb6428e2e702f 100755
--- a/setup.py
+++ b/setup.py
@@ -5,7 +5,7 @@ with open("README.md", "r") as fh:
 
 setuptools.setup(
     name = "impetuous-gfa",
-    version = "0.21.2",
+    version = "0.21.4",
     author = "Richard Tjörnhammar",
     author_email = "richard.tjornhammar@gmail.com",
     description = "Impetuous Quantification, Alignments, Enrichments and Group Analysis",
diff --git a/src/impetuous/reducer.py b/src/impetuous/reducer.py
index aa51f5fb4f235b3cc05411026797c960e3f85f21..240693d05d9e741553ee100deb77aad3285213d9 100644
--- a/src/impetuous/reducer.py
+++ b/src/impetuous/reducer.py
@@ -39,9 +39,9 @@ def padded_rolling_average( tv , tau ) :
 		return ( tv )
 	if len(tv)<tau :
 		return ( [ np.mean(v_) for v_ in tv ] )
-	centered = lambda x:(np.min(x),np.max(x)) ; N=len(tv);
+	centered = lambda x:(x[0],x[1]) ; N=len(tv);
 	w = int(np.floor(np.abs(tau)*0.5)) ;
-	jid = lambda i,w,N:[int((i-w)>0)*(i-w)%N,i,int(i+w<N)*(i+w)%N+int(i+w>=N)*(N-1)]
+	jid = lambda i,w,N:[int((i-w)>0)*((i-w)%N),int(i+w<N)*((i+w)%N)+int(i+w>=N)*(N-1)]
 	idx = [ centered( jid(i,w,N) ) for i in range(N) ]
 	mvalues = [ np.mean(tv[i[0]:i[1]]) for i in idx ]
 	return ( mvalues )