dada2 icon indicating copy to clipboard operation
dada2 copied to clipboard

removePrimers misses valid primer matches in certain circumstances

Open gmagoon opened this issue 4 years ago • 5 comments

Currently, using removePrimers with orient=TRUE appears to miss valid primer matches in certain circumstances. Specifically, it seems that in the current removePrimers logic, if there is a fwd-only primer match on the original (non-reverse-complement) sequence, a full (fwd+rev) match on the reverse complement will not be considered. This can lead to unexpected behavior, e.g. "loosening" of matching thresholds can result in fewer passing read sequences. I believe #956 should fix this.

gmagoon avatar Apr 19 '20 22:04 gmagoon