pairToEvents {amplican} | R Documentation |
Parse EMBOSS needle (or needleall) "pair" format into GRanges representation with events of deletions, insertions and mismatches. Make sure that each file corresponds to single subject (single amplicon). Assumes that bottom sequence "-bsequence" corresponds to the "subject" and full sequence alignment is returned.
pairToEvents(file, ID = "NA", strand_info = "+")
file |
(character) File path. |
ID |
(character) ID of the experiment, will be used as seqnames of the reutner ranges. |
strand_info |
(character) Strand to assign. |
(GRanges
) Same as events.