Tutorial 0: Customized FrogCap probe set

Create your own customized FrogCap probe set

This page will guide you on creating your own customized FrogCap probe set. You can draw markers from the growing FrogCap database that have been previously tested. Functionality includes: (1) The ability to select markers from certain families, subfamilies, and genera; (2) Different types of markers (UCEs, exons); (3) and marker features such as proportion of informative sites and marker length.

To begin, download the "Custom-Probe-Set" folder in the main GitHub directory to your computer. Next, open the folder and edit the R script titled: "FrogCap_Probe-Select.R", which will appear like so:

###############################################################################
########   Step 1: Create locus summary dataset               #################
###############################################################################
###############################################################################
###############################################################################

###### Directory setup ########

work.dir = "/Full/Directory/Path/to/Custom-Probe-Set"
output.file ="Custom_ProbeSet"

#The two metadata files should be placed in the working directory 
tax.data.file = "Taxonomy_dataset.csv"
marker.data.file = "Marker_Bait_dataset.csv"

#################################################
#Features you would like the probe set to have
#################################################

probe.size = 20000 #Number of probes, maximum is currently 20k for custom kits. use full kit
uce.include = TRUE #Other option FALSE, no UCEs
most.informative = TRUE #FALSE = randomly selects markers rather than by information content
min.exon.length = 200 #The minimum exon length
max.exon.length = 100000 #The maximum exon length

#Select a family or subfamily from the two lists below
clade.select = "Microhylidae"

#You can also select two or more taxonomic groups [overrides previous clade select]
#The program will select markers present in BOTH. 
#Must both be the same phylogenetic scale (Family / Subfamily)
#Cannot be from different Superfamilies (Ranoidea / Hyloidea)
clade.select = c("Mantellidae", "Rhacophoridae")

This section is the only section that needs editing. Namely, the directories should be the full path to make things go smoother. You can set various parameters here for the probe set. When finished editing, run the R script in the command line:

> Rscript FrogCap_Probe-Select.R

Which should have generated 3 files:

     /Custom-Probe-Set
      ├── FrogCap_Probe-Select.R
      ├── Taxonomy_dataset.csv
      └── /Marker_Bait_dataset.csv