findLQR: Find Off-targets and Fragmented alignments from reads.
In amplican: Automated analysis of CRISPR experiments

Description Usage Arguments Value See Also Examples

Will try to detect off-targets and low quality alignments (outliers). It tries k-means clustering on normalized number of events per read and read alignment score. If there are 3 clusters (decided based on silhouette criterion) cluster with high event count and low alignment score will be marked for filtering. When there is less than 1000 scores in aln it will filter nothing.

1	findLQR(aln)

aln

(data.frame) Should contain events from alignments in GRanges style with columns eg. seqnames, width, start, end, score.

(logical vector) where TRUE indicates events that are potential off-targets or low quality alignments.

findPD findEOP

Other filters: findEOP, findPD

file_path <- system.file("extdata", "results", "alignments",
                         "raw_events.csv", package = "amplican")
aln <- data.table::fread(file_path)
aln <- aln[seqnames == "ID_1"] # for first experiment
findLQR(aln)