The aim of this study was to detect the position of genomic regions associated with milk production and milk components in Pinzgau cattle. The dataset consisted of milk yield records of in total 7729 cows of 35 sires representing active bulls (19 animals) and DNA material stored in gene bank (16 animals) of Pinzgau cattle in Slovakia. In total 130087 test-day records of milk, protein and fat yield were used for the association analysis. The Illumina BovineSNP50 BeadChip V2 was used and after quality control final dataset consisted of 41,487 autosomal loci covering overall length 2,500,315 kb of the genome. Identification of genomic regions associated with milk production and composition were performed using GEMMA software with use of linear mixed model approach based on genetic-relationship matrix estimated from SNP genotypes to model correlation between the phenotypes. Based on this were found SNPs in the regions of important QTLs significantly associated with milk yield, dressing percentage, protein yield, SCS score, marbling score and fat yield.