DNA sequencing and restriction mapping provide us with information on DNA sequence evolution within populations, from which the phylogenetic relationships among the sequences can be inferred. Mutations such as base substitutions, deletions, insertions and transposable element insertions can be identified in each sequence. Theoretical study of this type of sequence evolution has been initiated recently. In this paper, population genetical models for sequence evolution under multiple types of mutation are developed. Models of infinite population size with neutral mutation, infinite population size with deleterious mutation and finite population size with neutral mutation are considered.