Hello! We are interested in processing sequences f...
# ampliseq
l
Hello! We are interested in processing sequences from the 515FY/926R PCR primer pair for both 18S and 16S amplicons. This requires a splitting step during the pipeline (if not, the 18S reads are lost as they do not overlap). Is this possible to do in ampliseq? Here is a pipeline describing it, and here is the paper. Thank you!
d
Afaik this is the intended application for --mergepairs_strategy consensus. A discussion about that addition is in https://github.com/nf-core/ampliseq/pull/803. Let me know if that serves your purpose.
Because I had a minute, I think indeed that approach should aid your approach:
This approach is particularly beneficial for datasets containing both prokaryotic and eukaryotic sequences, where lengths vary, leading to differing overlap extents.
n
Hi
--mergepairs_strategy consensus
has been indeed implemented to deal with this primer pair, we've detailed the approach in the supplementary of this preprint: https://doi.org/10.1101/2025.08.11.669668 , the paper is currently under review
l
Thank you so much both, that's great! I will try it out 🙏