DNA_reads_Per_Lane: Estimate per-lane reads

In angelgr2/radseq_tools: Tools for the Exploratory Analysis of RADseq Data

Description

This function calculates how many reads are required from the sequencing process in order to obtain a desired coverage for a given number of samples and cutsites. This function is to assist in determining the necessary lanes of sequencing required, and/or the Illumina platform to use for a given experiment.

Usage

DNA_Reads_Per_Lane(num_cutsites,
                   num_samples,
                   desired_coverage = NULL)

Arguments

num_cutsites	Number of cutsites present in the genome. This number can be the output from number_cutsites().
num_samples	Number of samples that the user wants to sequence in a single Illumina lane.
desired_coverage	Number of the coverage the user wishes to keep. If no value is provided, the default value is 30.

Value

Returns a single numberic value.

Author(s)

Kira Long

See Also

number_cutsites

Examples

#Define parameters
numberCuts <- 49000
samples <- 96
coverage <- 40

#Run arguments in Function
DNA_Reads_Per_Lane(numberCuts, samples, coverage)

angelgr2/radseq_tools documentation built on May 15, 2019, 3:59 a.m.