|
insertions[[length(new.seq)+1]] <- c( |
|
substr(seq, left, left+len-1), |
|
substr(qual, left, left+len-1) |
|
) |
Given that new.seq is a single string which grows through a series of paste0() calls, is the use of length(new.seq) supposed to be nchar(new.seq)? Currently it looks like the Insertions attribute always generates a NULL element followed by a single pair of insertions.