convControl: Conversion control

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Function to identify all potential conversion controls in a given input sequence, for a given list of fragments

Usage

1
convControl(sequence, fragmentation, multiple = FALSE)

Arguments

sequence

Nucleotide sequence input as a character string

fragmentation

List of MassArrayFragment objects corresponding to the fragmentation pattern of the sequence input

multiple

Logical value specifying whether or not to include multiple CGs on the same conversion control fragment where possible (default is FALSE).

Details

Potential conversion controls are identified from the nucleotide sequence input through pattern recognition of fragments that contain non-CG cytosines. Any conversion controls that contain CG dinucleotides or have a molecular weight outside of the usable mass window are screened out. Additionally, conversion controls that are located in identified primer sequence or have molecular weight overlap with other fragments are removed from consideration. Lastly, if the consideration of the fragment as a conversion control will cause new molecular weight overlap(s) with one or more other fragments, the control is also removed from consideration.

Value

Returns a list of MassArrayFragment objects identical to the input, with the exception that conversion controls are labeled and updated accordingly.

Author(s)

Reid F. Thompson (rthompso@aecom.yu.edu), John M. Greally (jgreally@aecom.yu.edu)

See Also

See Also as MassArrayFragment

Examples

1
2
3
4
data(MassArray.example.data)
MassArray.example.data$fragments.T[[54]]
conversion.data <- convControl(MassArray.example.data$sequence, MassArray.example.data$fragments.T)
conversion.data[[54]]

MassArray documentation built on Nov. 8, 2020, 5:16 p.m.