17Low or uneven read depth is a common limitation of genotyping-by-sequencing (GBS) 18 and restriction site-associated DNA sequencing (RAD-seq), resulting in high missing data rates, 19 heterozygotes miscalled as homozygotes, and uncertainty of allele copy number in heterozygous 20 polyploids. Bayesian genotype calling can mitigate these issues, but previously has only been 21 implemented in software that requires a reference genome or uses priors that may be 22 inappropriate for the population. Here we present several novel Bayesian algorithms that 23 estimate genotype posterior probabilities, all of which are implemented in a new R package, 24 polyRAD. Appropriate priors can be specified for mapping populations, populations in Hardy-25Weinberg equilibrium, or structured populations, and in each case can be informed by genotypes 26 at linked markers. The polyRAD software imports read depth from several existing pipelines, 27 and outputs continuous or discrete numerical genotypes suitable for analyses such as genome-28 wide association and genomic prediction. 29 30 Approximately 70% of vascular plant species are recent polyploids, yet genomic 31 resources and bioinformatics tools for polyploids typically lag behind those for diploids (Moghe 32 and Shiu 2014; Renny-Byfield and Wendel 2014; Bourke et al. 2018b). Reduced representation 33 DNA sequencing methods, such as genotyping-by-sequencing (GBS) and restriction site-34 associated DNA sequencing (RAD-seq), have made high-density genotyping considerably more 35 accessible and affordable (Poland and Rife 2012; Davey et al. 2013). However, the two most 36 popular pipelines for processing GBS and RAD-seq data, Stacks (Catchen et al. 2013) and 37 TASSEL (Glaubitz et al. 2014), do not output polyploid genotypes. Though pipelines for 38polyploids are available, each have limitations that prevent their general application. For 39 example, the UNEAK pipeline is designed for diploidized polyploids only (Lu et al. 2013). 40 HaploTag is specialized for self-fertilizing polyploids (Tinker et al. 2016). FreeBayes and 41 GATK can output polyploid genotypes, but require a reference genome (McKenna et al. 2010; 42 Garrison and Marth 2012). The software EBG imports read depth from other pipelines to 43 estimate auto-or allopolyploid genotypes (Blischak et al. 2018) but requires allele frequency 44 estimations from the parent species for allopolyploids. The R package updog estimates 45 polyploid genotypes from read depth, modeling preferential pairing and accounting for multiple 46 technical issues that can arise with sequencing data, and can output posterior mean genotypes 47