Real-time quantitative polymerase chain reaction (RT-qPCR) is the most common approach to quantify changes in gene expression. Appropriate internal reference genes are essential for normalization of data of RT-qPCR. In the present study, we identified suitable reference genes for analysis of gene expression in rice seedlings subjected to different heavy metal stresses such as deficiencies of iron and zinc and toxicities of cobalt, cadmium and nickel. First, from publically available RNA-Seq data we identified 10 candidate genes having stable expression. We also included commonly used house-keeping gene OsUBQ5 (Ubiquitin 5) in our analysis. Expression stability of all the 11 genes was determined by two independent tools, NormFinder and geNorm. Our results show that selected candidate reference genes have higher stability in their expression compared to that of OsUBQ5. Genes with locus ID LOC_Os03g16690, encoding an oxysterol-binding protein (OsOBP) and LOC_Os01g56580, encoding Casein Kinase_1a.3 (OsCK1a.3) were identified to be the most stably expressed reference genes under most of the conditions tested. Finally, the study reveals that it is better to use a specific reference gene for a specific heavy metal stress condition rather than using a common reference gene for multiple heavy metal stress conditions. The reference genes identified here would be very useful for gene expression studies under heavy metal stresses in rice.