R scripts to estimate locus, allele and SNP from RADseq data processed with Stacks