FindMaxRepeatDel: Return the number of repeat units in which a deletion is...
In ICAMS: In-Depth Characterization and Analysis of Mutational Signatures ('ICAMS')

Description Usage Arguments Details Value ID classification Examples

Return the number of repeat units in which a deletion is embedded

1	FindMaxRepeatDel(context, rep.unit.seq, pos)

`context`	A string that embeds `rep.unit.seq` at position `pos`
`rep.unit.seq`	A substring of `context` at `pos` to `pos + nchar(rep.unit.seq) - 1`, which is the repeat unit sequence.
`pos`	The position of `rep.unit.seq` in `context`.

This function is primarily for internal use, but we export it to document the underlying logic.

For example FindMaxRepeatDel("xyaczt", "ac", 3) returns 0.

If substr(context, pos, pos + nchar(rep.unit.seq) - 1) != rep.unit.seq then stop.

If this functions returns 0, then it is necessary to look for microhomology using the function FindDelMH.

Warning
This function depends on the variant caller having "aligned" the deletion within the context of the repeat.

For example, a deletion of CAG in the repeat

1	GTCAGCAGCATGT

can have 3 "aligned" representations as follows:

1
2
3

CT---CAGCAGGT
CTCAG---CAGGT
CTCAGCAG---GT

In these cases this function will return 2. (Please not that the return value does not include the rep.uni.seq in the count.)

However, the same deletion can also have an "unaligned" representation, such as

1	CTCAGC---AGGT

(a deletion of AGC).

In this case this function will return 1 (a deletion of AGC in a 2-element repeat of AGC).

The number of repeat units in which rep.unit.seq is embedded, not including the input rep.unit.seq in the count.

See https://github.com/steverozen/ICAMS/raw/master/data-raw/PCAWG7_indel_classification_2017_12_08.xlsx for additional information on ID (small insertion and deletion) mutation classification.

See the documentation for Canonicalize1Del which first handles deletions in homopolymers, then handles deletions in simple repeats with longer repeat units, (e.g. CACACACA, see FindMaxRepeatDel), and if the deletion is not in a simple repeat, looks for microhomology (see FindDelMH).

See the code for unexported function CanonicalizeID and the functions it calls for handling of insertions.