added --subset argument: now only a subset of snps can only be analyzed