We present a consistent treatment of heavy quarks for jet production in DIS at NLO accuracy. The method is based on the ACOT massive factorization scheme and dipole subtraction method for jets. The last had to be however extended in order to take into account initial state splittings with heavy quarks. We constructed relevant kinematics and dipole splitting functions together with their integrals. We partially implemented the method in a MC program and checked against the known inclusive result for charm structure function.