diff --git a/openmc/data/neutron.py b/openmc/data/neutron.py
index 670144ec1..bcb830aa0 100644
--- a/openmc/data/neutron.py
+++ b/openmc/data/neutron.py
@@ -622,7 +622,7 @@ class IncidentNeutron(EqualityMixin):
         data.energy[strT] = energy
         total_xs = ace.xss[i + n_energy : i + 2*n_energy]
         absorption_xs = ace.xss[i + 2*n_energy : i + 3*n_energy]
-        heating_number = ace.xss[i + 3*n_energy : i + 4*n_energy]*EV_PER_MEV
+        heating_number = ace.xss[i + 4*n_energy : i + 5*n_energy]*EV_PER_MEV
 
         # Create redundant reactions (total, absorption, and heating)
         total = Reaction(1)
diff --git a/tests/regression_tests/tallies/results_true.dat b/tests/regression_tests/tallies/results_true.dat
index 9830ff872..945e20fc6 100644
--- a/tests/regression_tests/tallies/results_true.dat
+++ b/tests/regression_tests/tallies/results_true.dat
@@ -1 +1 @@
-de773a799f84348241ebd65bf7beae8114861d8674b652bfd443d578c146a7d37ca38f2efc5d05124fdbb1cf12829907768351e6b2d6b5f4bd2c121417757507
\ No newline at end of file
+cc3993048a54cb162a3224d69c5734df14e17a4c8784a9b8781537a5a5ecb3d9b8145c3b0e95957e67a68dc098e8ce1bc6c61d39e9125d0d0327a659f4b49808
\ No newline at end of file