Parameter

Parameter

Name

Variable

Default

Value

Parameter

Range

Description
PCR_DISTRIBUTION
 default{default, none, file}parameter specifies the name of a PCR distribution file, or 'default' to use a distribution with 15 rounds and 20 bins. A value of 'none' disables PCR amplification.
PCR_PROBABILITY
0.7

duplication probability in each step of the simulated PCR. The parameter value is only taken into account if GC_MEAN is 'NaN'.

GC_MEAN
meanGC0.5NaN, 0 < meanGC < 1

mean of duplication probability distribuiton with respect to GC content; the distribution is approximated by a normal distribution. A value of 'NaN' disables the GC-dependency of PCR and employs the (constant) probability PCR_PROBABILITY

GC_SD
SDGC 0.1 0 < SDGC < 1standard deviation of duplication probability distribuiton with respect to GC content

Algorithm

The efficiency of the polymerase chain reaction (PCR) amplification is either specified by an universal success rate , or, by a normal distribution  parameterized to capture GC preferential biases (defaultmeanGC=0.5 and SDGC=0.1). Given p, the number of copies produced from a certain fragment is determined by random sampling under the cumulative binomial:

 

with S denoting the PCR cycle and N the number of molecules. As default, we assume 15 PCR cycles (S=15), and sample randomly the number of duplicates yielded by PCR amplification under the corresponding probability distribution P15(N) for all possible values of N=[1;215]. The recursion terminates by

  • No labels