diff --git a/pyemma/coordinates/data/featurizer.py b/pyemma/coordinates/data/featurizer.py
index <HASH>..<HASH> 100644
--- a/pyemma/coordinates/data/featurizer.py
+++ b/pyemma/coordinates/data/featurizer.py
@@ -447,7 +447,7 @@ class BackboneTorsionFeature(DihedralFeature):
             self._psi_inds = indices[np.in1d(indices[:, 1],
                                              topology.select(selstr), assume_unique=True)]
 
-        dih_indexes = np.vstack((self._psi_inds, self._psi_inds))
+        dih_indexes = np.vstack((self._phi_inds, self._psi_inds))
         super(BackboneTorsionFeature, self).__init__(topology, dih_indexes,
                                                      deg=deg, cossin=cossin)
[coor/featurizer] fix bug in BackboneTorsions (returned (psi, psi) for each residue).
Luckily this has not been released in the prior version.