A function to compute the percentage of covered random sequences for a given number of flows.

Share:

Description

Evaluation of an exact expression to obtain the percentage of all random nucleotide sequences of a given length (n) that can be sequenced with semiconductor sequencing given a certain number of flows (k).

Usage

1
flowsRandom(k = 300, n = 150)

Arguments

k

Integer. Number of flows available during semiconductor sequencing.

n

Integer. Length of the random sequences to be sequenced.

Value

P

Percentage of sequences that are completely covered

Q

Percentage of sequences that are not completely covered

delta

1 - (P + Q), should be a small number when calculation is accurate

N

Number of summands for calculation of P

M

Number of summands for calculation of Q

summands

N + M = total number of summands

Author(s)

Michael Bockmayr and Jan Budczies

Examples

1
2
3
### To compute the percentage of covered sequences of length 150 with 300 flows run:

flowsRandom(300,150)