This seems to be an absurd question, but bear with me.

In quantum field theory, the Dirac fermion mass Lagrangian term reads
$$
m\bar\psi \psi = m(\bar\psi\_L \psi\_R + \bar\psi\_R \psi\_L) = m(\psi\_L^\dagger \gamma\_0\psi\_R + \psi\_R^\dagger \gamma\_0\psi\_L)
$$
Assuming this fermion mass Lagrangian term is used as the integrand (in the exponential) in functional integral formalism of Dirac-fermion-related QFT, in the following we will regard $\psi$ components as anti-commuting Grassmann variables given the anti-symmetric nature of the fermions, rather than classical commuting complex variables.

Note that $m$ is supposed to be real throughout this post. The question here is about whether the mass Lagrangian term $m\bar\psi \psi$ is real or imaginary, not about $m$ parameter itself. (Note that there could be legitimate pseudoscalar mass term as in $m\bar{\psi} i\gamma\_5\psi$. But we will not discuss about pseudoscalar mass in the current post. The interested readers can see [here](https://physics.stackexchange.com/questions/611255/is-there-a-square-root-of-klein-gordon-equation-which-is-different-from-dirac#611255) for more details.)

However, if you look under the hood of the familiar mass term $m\bar\psi \psi$, there are some surprises lurking around. Let's look at a simple example in the Weyl basis
$$
\psi = (\xi, 0, \chi, 0)^T
$$
where $\xi$ and $\chi$ are *real Grassmann numbers* ($\xi^\*= \xi$, $\chi^\*= \chi$, they are not 2 component columns, we will get to the general case of **complex** Grassmann numbers towards the end of the post). In the Weyl basis, $(\xi, 0, 0, 0)^T$ and $(0, 0, \chi, 0)^T$ represent left ($\psi\_L$) and right ($\psi\_R$) handed part of the Dirac spinor, respectively.

Let's calculate the mass term:
$$
m\bar\psi \psi = m(\psi\_L^\dagger \gamma\_0\psi\_R + \psi\_R^\dagger \gamma\_0\psi\_L)
= m(\xi^\*\chi + \chi^\*\xi) = m(\xi\chi + \chi\xi) = 0
$$
Oops, it's identical to zero since $\xi$ and $\chi$ are anti-commuting Grassmann numbers! Note that the $\gamma\_0$ matrix just flip the left-handed components to the right-handed components and *vice versa* in the Weyl basis.

Now, let's tweak the the trial spinor to make one of its component *imaginary Grassmann* (multiplying $\chi$ by an $i$)
$$
\psi = (\xi, 0, i\chi, 0)^T
$$
Low and behold (I will omit the $m$ parameter hereafter):
$$
\bar\psi \psi = \psi\_L^\dagger \gamma\_0\psi\_R + \psi\_R^\dagger \gamma\_0\psi\_L
= \xi^\*(i\chi) + (i\chi)^\*\xi = i\xi\chi - i\chi\xi = 2i\xi\chi \neq 0
$$
The mass term is non-zero.

**The remarkable and weird feature is that the mass term is imaginary with an $i$!**

The interested reader can try all sorts of $\psi$ configurations in any representation (Weyl basis or not), and you will end up with the same result of imaginary mass. Any efforts of building a real mass term will be a wild goose chase, since the $\psi\_L^\dagger \gamma\_0\psi\_R$ portion will always cancel out the $\psi\_R^\dagger \gamma\_0\psi\_L$ portion.

Let's double check whether the imaginary mass term is Hermitian:
$$
(2i\xi\chi)^\dagger = -2 i\chi^\*\xi^\* = -2i\chi\xi = 2i\xi\chi
$$
So the Hermitian property is secured.

On the other hand, a real mass term (if it existed)
$$
(2\xi\chi)^\dagger = 2 \chi^\*\xi^\* = 2\chi\xi = -2\xi\chi
$$
would be non-Hermitian.

The key point here is that the Hermitian operator is by definition
$$
(AB)^\dagger = B^\dagger A^\dagger
$$
Note that there is NO minus sign, even if both $A$ and $B$ are Grassmann odd. (As a side note, transpose is defined as: $(AB)^T = -B^T A^T$, if $A$ and $B$ are Grassmann valued. Note that there is a minus sign! See [here](https://physics.stackexchange.com/questions/388079/majorana-flip-relations?noredirect=1&lq=1).)

At the end of the day, physicists seem not troubled by the imaginary nature of the mass term as long as it's Hermitian. **I have to underscore (in response to @octonion's comments) that being Hermitian and being real are two disparate notions.**

You might wonder why the imaginary mass is not mentioned in the usual text books. It's because when we deal fermions, the common practice is to use *complex* Grassmann numbers
$$
\xi = \xi\_1 + i\xi\_2 \\
\chi = \chi\_1 + i\chi\_2
$$
where $\xi\_1$, $\xi\_2$, $\chi\_1$, and $\chi\_2$ are *real* Grassmann numbers

And thus the mass term of $\psi = (\xi, 0, \chi, 0)^T$ is
$$
\bar\psi \psi = \psi\_L^\dagger \gamma\_0\psi\_R + \psi\_R^\dagger \gamma\_0\psi\_L
= \xi^\*\chi + \chi^\*\xi
$$
the imaginary nature is hidden in plain sight. Only when we write out the explicit terms
$$
\xi^\*\chi + \chi^\*\xi = (\xi\_1 + i\xi\_2)^\*(\chi\_1 + i\chi\_2) + (\chi\_1 + i\chi\_2)^\*(\xi\_1 + i\xi\_2) = 2i(\xi\_1\chi\_2 + \chi\_1\xi\_2)
$$
the imaginary mass is manifest.

---

Added note in response to the wrong answer below (by @alexarvanitakis) saying that "The signs and/or presence of factors of i in fermion lagrangians is somewhat superfluous and convention dependent".

Of course, one can have a Dirac Lagrangian with no $i$. For example, one can simply changing the metric from (+, -, -, -) to (-, +, +, +), see more details [here](https://physics.stackexchange.com/questions/513768/dirac-equation-without-i/513772#513772). However, using what ever convention you choose, you still end up with a imaginary mass term! This is because that each column element of the Dirac wave function is valued in **complex** Grassmann space. Whereas a **real** Grassmann valued Dirac wave function will imply zero mass.

The wrong answer further says that "For example you can work with a convention where complex conjugation does not reverse the order of a product of fermions which changes radically the appearances of factors of i".

Please note that in deriving the imaginary mass above, there is no reverse invoked for any product of fermions. So the mass term is still imaginary. The only place the complex conjugation convention in reverse of a product of fermions is relevant is the proof that the mass term is Hermitian, albeit imaginary. If one takes the wrong answer's convention, the mass term will be both non-Hermitian and imaginary!

Additionally, the wrong answer says that "You instead want to look at the Klein-Gordon equation satisfied by the fermion field...You need to correlate the sign of $m^2$ with the convention for the gamma matrices so that the above operator does not admit tachyonic-type solutions."

I am talking about $m\bar\psi \psi$ being imaginary, not $m$ being imaginary. The wrong answer's proof of $m$ being real (or $m^2$ being positive) is totally irrelevant to the question here!

---

More added note:

I am not talking about whether the expectation value $\langle m\bar\psi \psi\rangle$ being real. That is not the problem. What I am talking about is the Lagrangian/action in integrand before path integral, while the expectation value is after path integral. Actually, even-though the Lagrangian mass term is imaginary, the path integral expectation value $\langle m\bar\psi \psi\rangle$ IS real. I am not contesting about the expectation value.