Searching for patterns in data is important because it can lead to the discovery of sequence segments that play a functional role. The complexity of pattern statistics that are used in data analysis and the need of the sampling distribution of those statistics for inference renders efficient computation methods as paramount. This article gives an overview of the main methods used to compute distributions of statistics of overlapping pattern occurrences, specifically, generating functions, correlation functions, the Goulden-Jackson cluster method, recursive equations, and Markov chain embedding. The underlying data sequence will be assumed to be higher-order Markovian, which includes sparse Markov models and variable length Markov chains as special cases. Also considered will be recent developments for extending the computational capabilities of the Markov chain-based method through an algorithm for minimizing the size of the chain's state space, as well as improved data modeling capabilities through sparse Markov models. An application to compute a distribution used as a test statistic in sequence alignment will serve to illustrate the usefulness of the methodology. This article is categorized under: Statistical Learning and Exploratory Methods of the Data Sciences > Pattern Recognition Data: Types and Structure > Categorical Data Statistical and Graphical Methods of Data Analysis > Modeling Methods and Algorithms K E Y W O R D S Auxiliary Markov chain, distribution of a pattern statistic, Markovian sequences, sparse Markov models, VLMC