Change back to unique SNP iteration
TODO: iterate whole BSGS data, get MAFs