We study a “heavy” QCD axion whose coupling to the standard model is dominated by <math><mi>a</mi><mi>G</mi><mover><mi>G</mi><mo>˜</mo></mover></math> but with <math><msub><mi>m</mi><mi>a</mi></msub><mo>≫</mo><msub><mi>m</mi><mi>π</mi></msub><msub><mi>f</mi><mi>π</mi></msub><mo>/</mo><msub><mi>f</mi><mi>a</mi></msub></math>. This is well motivated as it can solve the strong <math><mi>C</mi><mi>P</mi></math> problem while evading the axion quality problem. It also poses interesting challenges for its experimental search due to its suppressed couplings to photons and leptons. Such axion with mass around a GeV is kinematically inaccessible or poorly constrained by most experimental probes except <math><mi>B</mi></math>-factories. We study <math><mi>B</mi><mo>→</mo><mi>K</mi><mi>a</mi></math> transitions as a powerful probe of the heavy QCD axion by performing necessary 2-loop calculations for the first time, together with some improvement on the existing analysis strategy. We find some of the existing limits are enhanced by at least an order of magnitude. We also demonstrate that the bounds are robust against unknown UV physics. For forthcoming data sets of the Belle II experiment, we provide a projection that <math><msub><mi>f</mi><mi>a</mi></msub></math> of a few TeV is within its future reach, which is relevant to the quality problem.