amplicanOverlap {amplican} | R Documentation |
To determine which deletions, insertions and mismatches (events) are probably created by CRISPR we check whether they overlap expected cut sites. Expected cut sites should be specified in UPPER CASE letters in the amplicon sequences.
amplicanOverlap(aln, cfgT, cut_buffer = 5, relative = FALSE)
aln |
(data.frame) Contains relevant events in GRanges style. |
cfgT |
(data.frame) Contains amplicon sequences. |
cut_buffer |
(numeric) Number of bases that should expand 5' and 3' of the specified expected cut sites. |
relative |
(boolean) Sets whether events are relative to the position of the target site. |
(bolean vector) Where TRUE means that given event overlaps cut site.
Other analysis steps: amplicanAlign
,
amplicanConsensus
,
amplicanFilter
, amplicanMap
,
amplicanNormalize
,
amplicanPipeline
,
amplicanReport
,
amplicanSummarize
file_path <- system.file("test_data", "test_aln.csv", package = "amplican") aln <- data.table::fread(file_path) cfgT <- data.table::fread( system.file("test_data", "test_cfg.csv", package = "amplican")) all(aln$overlaps == amplicanOverlap(aln, cfgT))