CFT consists for most mathematicians - who are interested in this topic - currently of the study of vertex operator algebras, see this question on math overflow:
You can find a little bit more about the topic and the work of several mathematicians here:
As you can see from the answers on mathoverflow, vertex algebras were not invented for the study of CFT, and that they form an axiomatic abstraction of operator algebra products was noted only later.
A personal and very subjective note: One should not underestimate the amount of theoretical physics that is necessary to understand what a QFT is to physicists. Most mathematicians that encounter physics for the first time since highschool through some QFT framework seem to be quite taken aback by the high intrance fee they'd have to pay to understand this. This is my own, personal explanation for the observtion that most mathematicians study the formal machinery only, in order to use it to prove some new mathematical theorems, but only very rarely in order to better understand what physicists do. Although you'll find quite a lot of work by quite a lot of rather famous mathematicians if you follow the links above, AFAIK there is none doing the work you describe.
Your statement
working with and subtracting infinities ... which in general have no mathematical meaning
is not really correct, and seems to have a common misunderstanding in it. The technical difficulties from QFT do not come from infinities. In fact, ideas basically equivalent to renormalization and regularization have been used since the beginning of math--see, e.g., many papers by Cauchy, Euler, Riemann, etc. In fact, G.H. Hardy has a book published on the topic of divergent series:
http://www.amazon.com/Divergent-AMS-Chelsea-Publishing-Hardy/dp/0821826492
There is even a whole branch of math called "integration theory" (of which things like Lebesgue integration is a subset) that generalizes these types of issues. So having infinities show up is not an issue at all, in a sense, they show up out of convenience.
So the idea that infinities have anything to do with making QFT axiomatic is not correct.
The real issue, from a more formal point of view, is that you "want" to construct QFTs via some kind of path integral. But the path integral, formally (i.e., to mathematicians) is an integral (in the general sense that appears in topics like "integration theory") over a pretty pathological looking infinite dimensional LCSC function space.
Trying to define a reasonable measure on an infinite dimensional function space is problematic (and the general properties of these spaces doesn't seem to be particularly well understood). You run into problems like having all reasonable sets being "too small" to have a measure, worrying about measures of pathological sets, and worrying about what properties your measure should have, worrying if the "$\mathcal{D}\phi$" term is even a measure at all, etc...
At best, trying to fix this problem, you'd run into an issue like you have in the Lebesgue integral's definition, where it defines the integral and you construct some mathematically interesting properties, but most of its utility is in letting you abuse the Riemann integral in the way you wanted to. Actually calculating integrals from the definition of the Lebesgue integral is not generally easy. This isn't really enough to attract the attention of too many physicists, since we already have a definition that works, and knowing all of its formal properties would be nice, and would certainly tell us some surprising things, but it's not clear that it would be all that useful generally.
From an algebraic point of view, I believe you run into trouble with trying to define divergent products of operators that depend on renormalization scheme, so you need to have some family of $C^*$-algebras that respects renormalization group flow in the right way, but it doesn't seem like people have tried to do this in a reasonable way.
From a physics point of view, we don't care about any of this, because we can talk about renormalization, and demand that our answers have "physically reasonable" properties. You can do this mathematically, too, but the mathematicians are not interested in getting a reasonable answer; what they want is a set of "reasonable axioms" that the reasonable answers follow from, so they're doomed to run into technical difficulties like I mentioned above.
Formally, though, one can define non-interacting QFTs, and quantum mechanical path integrals. It's probably the case that formally defining a QFT is within the reach of what we could do if we really wanted, but it's just not a compelling topic to the people who understand how renormalization fixes the solutions to physically reasonable ones (physicists), and the formal aspects aren't well-understood enough that it's something one could get the formalism for "for free."
So my impression is that neither physicists or mathematicians generally care enough to work together to solve this problem, and it won't be solved until it can be done "for free" as a consequence of understanding other stuff.
Edit:
I should also add briefly that CFTs and SCFTs are mathematically much more carefully defined, and so a reasonable alternative to the classic ideas I mentioned above might be to start with a SCFT, and define a general field theory as some kind of "small" modification of it, done in such a way to keep just the right things well-defined.
Best Answer
Lack of convergence does not mean there is nothing mathematically rigorous one can extract from perturbation theory. One can use Borel summation. In fact, Borel summability of perturbation theory has been proved for some QFTs:
In fact these articles obtain such results by using an alternative to ordinary perturbation theory called a multiscale (or phase cell or phase space) cluster expansion. The latter is based on combinatorial structures which mimic Feynman diagrams. However, these expansions converge at small coupling.
Edit as per Timur's comment: Glimm and Jaffe's book is what you want to read in order to understand why one needs cluster expansions. It is excellent at giving the big picture: how axiomatic, Euclidean, constructive QFTs fit together, as well as with scattering theory. But for learning how to do a cluster expansion the book is outdated. The cluster expansion explained in GJ is the early one invented by Glimm, Jaffe and Spencer in their Annals of Math article. It was the first in the QFT context and as such quite a mathematical feat. However there has been many improvements and simplifications since then (around 1973). If you want to learn about cluster expansions in 2011, here is a more efficient path: