Hi, I have a question regarding how I should configure the
nf-core Ampliseq pipeline for
16S metabarcoding analysis targeting the V3–V4 region. I used the
Zymo Quick-16S Plus NGS Library Prep Kit (V3–V4) (
https://zymoresearch.eu/products/quick-16s-plus-ngs-library-prep-kit-v3-v4-udi?srsltid=AfmBOooKXo-LX85v4ocNVcKSIs2Z9jSEPgg7fvqERGKS1Mv_kLrgZrtz) and sequenced the libraries on an
Illumina MiniSeq.
The sequencing was performed in *paired-end mode*; however, one limitation of this platform is that it produces reads of
2 × 150 bp. This detail appears to strongly affect the
read merging step during the DADA2 stage. I used the
default quality-based trimming parameters, and the resulting DADA2 statistics are shown in the file
DADA2_stats.tsv. I am analyzing
butter cheese (queso mantecoso) samples, and I am also attaching the script that I used. After observing the merging issue, I would like to know how I should proceed with the pipeline—whether it would be more appropriate to analyze the data as
single-end (SE) or
paired-end (PE) reads.
This is particularly important because the analysis will be part of a
manuscript that we plan to publish, and this is the first time I have encountered this situation. I would greatly appreciate your guidance on how to proceed.
Thank you in advance for your support.