Right now this computation is very slow, which is a big problem when dealing with thousands of markers (e.g. RADseq).