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