We present a Caputo fractional order mathematical model that describes the cellular infection of the Hepatitis B virus and the immune response of the body with Holling type II functional response. We study the existence of unique positive solutions and the local and global stability of virus-free and endemic equilibria. Finally, we present numerical results using the Adam-type predictor–corrector iterative scheme.