Proof of concept implementation of a highly parallelisable variant-calling algorithm based on the Li and Stephens model. The model is a simplified version of what IMPUTE2 has, works with a Gibbs Sampler, and needs ~700 iterations to converge, but is still about 3x faster in MATLAB than IMPUTE2 in C.