small change to script + added --subset argument: now only a subset of snps can only be analyzed