Description Usage Arguments Value
It uses formula p(D=g|G_1=g)=1-e,p(D=other\ call|G_1=g)=e/3 to generate the read for a given genotype G_1 and error rate e.
It is called from the generate_seqdata
functions.
1 | seq_call(genotype, error, ndepth)
|
genotype |
vector of two alleles for a given genotype, ie. genotype CT is given as ('C', 'T') |
error |
vector of sequencing errors, length of this vector is the same as read depth (ndepth). |
ndepth |
the read depth (integer) |
a row vector of sequence reads consisting of 'A','G','T' or 'C' for a single variant with length the same as read depth (ndepth).
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.