diff --git a/testsuite/MDAnalysisTests/analysis/.coverage.PN126116.55116.363419 b/testsuite/MDAnalysisTests/analysis/.coverage.PN126116.55116.363419
new file mode 100644
index 00000000000..6670566d664
Binary files /dev/null and b/testsuite/MDAnalysisTests/analysis/.coverage.PN126116.55116.363419 differ
diff --git a/testsuite/MDAnalysisTests/analysis/test_hydrogenbonds_analysis.py b/testsuite/MDAnalysisTests/analysis/test_hydrogenbonds_analysis.py
index eb1a87c1af1..e572f39e9c3 100644
--- a/testsuite/MDAnalysisTests/analysis/test_hydrogenbonds_analysis.py
+++ b/testsuite/MDAnalysisTests/analysis/test_hydrogenbonds_analysis.py
@@ -365,9 +365,9 @@ def test_no_bond_donor_sel(self, universe):
         u.add_TopologyAttr('mass', [15.999, 1.008, 1.008] * n_residues)
         u.add_TopologyAttr('charge', [-1.04, 0.52, 0.52] * n_residues)
         h = HydrogenBondAnalysis(u, **kwargs)
-        pairs = h._get_dh_pairs()
+        donors = u.select_atoms(h.guess_donors())
 
-        assert len(pairs) == 2
+        assert len(donors) == 2
 
     def test_first_hbond(self, hydrogen_bonds):
         assert len(hydrogen_bonds.results.hbonds) == 2