Considering the ever growing world population as well as climate change, it has become essential to improve crops in order to reach higher yields, more tolerance to pests and droughts, and better adaptation to different environments. The advancement of next-generation sequencing technology enables us to perform comparative genome analyses for related genomes, specifically soybean, other Glycine species, and other related legumes. Although plants and their chemistry have been under study for a very long time, the understanding of the value of their secondary metabolites and their practical application is rather a newer process. One of the most studied pathways in legumes is the flavonoid biosynthetic pathway. Flavonoids play many important roles in plants and their interactions with their environment. Despite all the previous studies on the flavonoid biosynthetic pathway, our knowledge of the connection between the genotype and phenotype of the pathway, and the gene families elaborating the different characteristics of different species, particularly the gene families determining the flower colors and seed coat colors of different species, is still very inadequate. As part of this study I have identified candidate gene families and their putative orthologs for twelve enzymes from the flavonoid pathway among eight legumes and Arabidopsis thaliana. The study provided data for looking for candidate genes for flower colors and seed coat colors in sixteen soybean varieties. Many of these pigments can be used as markers to study different biological and evolutionary processes within the legumes. Isoflavone synthase, one of the enzymes of this pathway, was targeted for re-sequencing of ~150,000 bpregions in seven perennial legumes, where I have annotated these targeted regions, identified orthologous genes and performed evolutionary analyses to understand the genomic dynamics between these species.