Multifractal detrended moving average analysis (MFDFA) has been performed on time series of three different epidemiological indices, namely daily number of new cases per million population, reproduction rate and positive rate, recorded during the Covid-19 pandemic. We consider India, Japan, USA and UK for this study. For all the series the scaling behavior of MFDMA fluctuation functions has been examined, and the multifractal variables like the generalised Hurst exponent spectra and singularity spectra have been calculated. It is found that the time series of the Covid-19 indices are of multifractal nature. The degree of multifractality of the reproduction rate series is significantly weaker than the other two series for all the countries. To find out the source(s) of multifractality, we generate a set of one hundred shuffled series and Amplitude-Adjusted Fourier Transform surrogate series for each of the original series. Comparing the results obtained from the original series with those from the shuffled and surrogate series we understand that the main source of multifractality in Covid-19 data is long-range temporal correlation among the series values.