Fig. 1

Schematic representation of the procedure used to estimate the number of cutting sites present in the Antarctic fur seal genome for the enzyme SbfI. a Two DNA fragments were generated for every restriction site (denoted by vertical red lines) present in the fur seal genome (denoted by a continuous black horizontal line) and the first ~ 200 bp of each of them was sequenced in multiple copies (horizontal green lines); b The resulting sequence data were delivered as a collection of raw sequence reads; c The Stacks pipeline was then employed to assemble the raw reads into RAD loci (represented by blue horizontal lines); d The consensus sequence of every RAD locus was mapped to the reference genome (denoted by the continuous grey horizontal line) allowing us to distinguish among five different scenarios: (I) a single pair of RAD loci map around a restriction enzyme recognition sequence; (II) multiple RAD loci map around a recognition sequence; (III) a single RAD locus maps to a recognition sequence; (IV and V) a pair or a single RAD locus maps to a genomic region not containing the recognition sequence. Scenarios I, II and III are indicative of the presence of a restriction enzyme cutting site that is represented by the assembled RADseq data