diff --git a/reacnetgenerator/tools.py b/reacnetgenerator/tools.py
index b9621f1eb..d47c13abe 100644
--- a/reacnetgenerator/tools.py
+++ b/reacnetgenerator/tools.py
@@ -36,7 +36,7 @@ def read_species(
     >>> from reacnetgenerator.tools import read_species
     >>> import matplotlib.pyplot as plt
     >>> step_idx, n_species = read_species('methane.species')
-    >>> plt.plot(step_idx, n_species['[H]C([H])([H])[H]'])
+    >>> plt.plot(step_idx, n_species['[H][C]([H])([H])[H]'])
     >>> plt.savefig("methane.svg")
     """
     step_idx = []