Similar to Fraglength2 we should have a function that extracts genotypes from shortread objects containing a single or multiple SNPs.