stringCoverage-extractedReadsListReverseComplement-method: Get string coverage STR identified objects.

Description Usage Arguments Value Examples

Description

stringCoverage takes an extractedReadsList-object and finds the coverage of every unique string for every marker in the provided list.

Usage

1
2
3
  ## S4 method for signature 'extractedReadsListReverseComplement'
stringCoverage(extractedReadsListObject,
  control = stringCoverage.control())

Arguments

extractedReadsListObject

an extractedReadsList-object, created using the identifySTRRegions-function.

control

an stringCoverage.control-object.

Value

Returns a list, with an element for every marker in extractedReadsList-object, each element contains the string coverage of all unique strings of a given marker.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# Regions identified using 'identifySTRs()'
data("identifiedSTRs")

# Limiting and restructuring
sortedIncludedMarkers <- sapply(names(identifiedSTRs$identifiedMarkersSequencesUniquelyAssigned),
                                function(m) which(m == flankingRegions$Marker))

# Aggregate the strings
stringCoverage(extractedReadsListObject = identifiedSTRs,
               control = stringCoverage.control(
                   motifLength = flankingRegions$MotifLength[sortedIncludedMarkers],
                   Type = flankingRegions$Type[sortedIncludedMarkers],
                   numberOfThreads = 1,
                   trace = FALSE,
                   simpleReturn = TRUE))

svilsen/STRMPS documentation built on May 31, 2019, 7:55 a.m.