Generation of biomolecule sequence coevolution data structures, matrices, scores, and sectors are described. Generally, the generated coevolution data removes covariant noise due to phylogenetic drift and can reveal coevolution of residue positions in multiple phylogenetic distances. Scores can be built upon the data structures and matrices to reveal sectors of residue positions that function and evolve together. Furthermore, the coevolution data structures, matrices, scores, and sectors can be used to predict structure or function of residue variants.