This paper presents a study of early epidemiological assessment of COVID-19 transmission dynamics in Indonesia. The aim is to quantify heterogeneity in the numbers of secondary infections. To this end, we estimate the basic reproduction number [Formula: see text] and the overdispersion parameter [Formula: see text] at two regions in Indonesia: Jakarta-Depok and Batam. The method to estimate [Formula: see text] is based on a sequential Bayesian method, while the parameter [Formula: see text] is estimated by fitting the secondary case data with a negative binomial distribution. Based on the first 1288 confirmed cases collected from both regions, we find a high degree of individual-level variation in the transmission. The basic reproduction number [Formula: see text] is estimated at 6.79 and 2.47, while the overdispersion parameter [Formula: see text] of a negative-binomial distribution is estimated at 0.06 and 0.2 for Jakarta-Depok and Batam, respectively. This suggests that superspreading events played a key role in the early stage of the outbreak, i.e., a small number of infected individuals are responsible for large numbers of COVID-19 transmission. This finding can be used to determine effective public measures, such as rapid isolation and identification, which are critical since delay of diagnosis is the most common cause of superspreading events.