Currently requires "aa" notation.
Usage
coverage_at_aa(coverage, aa)
Arguments
- coverage
A data frame with columns labelled "position" (0-indexed) and "coverage".
- aa
A vector of mutations. aa:orf1a:I300V, ins:28215:3, del:27378:25, or C703T
Details
Finds the maximum coverage for the three positions of the amino acid. This ensures that the number of observed mutations can never be larger than the coverage.