The study will be administered to 270 children with reading difficulties from 10 different schools (27 children per school). 135 children will have music lessons and 135 will not; therefore, 5 schools will be the control and 5 schools will be the intervention group. The analysis will consider the cluster structure, since the randomization was not conducted at the individual level (i.e., the school level was the randomization unit). For the inferential analysis, generalized estimation models and structural equation modeling will be used.