Samples_Per_Lane: Estimate the number of samples per single Illumina lane
In angelgr2/radseq_tools: Tools for the Exploratory Analysis of RADseq Data

Description Usage Arguments Value Author(s) See Also Examples

This function calculates the estimated maximum number of samples a user can multiplex in one Illumina lane given the number of cutsites, while maintaining a desired user-defined coverage.

1
2
3

Samples_Per_Lane(num_cutsites,
                 desired_coverage = NULL,
                 sequencing_machine = NULL)

`num_cutsites`	Number of cutsites present in the genome. This number can be the output from `number_cutsites()`
`desired_coverage`	Number of the coverage the user wishes to keep. If no value is provided, the default value is `30`.
`sequencing_machine`	Illumina sequencing platform that the user wants to use for sequencing. Argument has two parameters, `hiseq2500` and `hiseq4000`, each has different read throughputs. If no machine is provided, the function will default to `hiseq2500`.

Returns a vector with an estimate of low, medium, and high estimates of the number of samples a user can fit in a single Illumina lane.

Kira Long

number_cutsites

#Define parameters
numberCuts <- 39000
coverage <- 30

#Run arguments in Function
Samples_Per_Lane(numberCuts,
                 coverage,
                 "hiseq4000")