c_malossi
(Camila Dantas Malossi)
December 7, 2023, 2:36pm
1
Hi,
I have a Illumina run (150x2 v2) and amplicon size is 291, with a overlap of 9nt. I cutoff my primers/adaptors and I use dada2 (denoise_pair) to denoising sequences. In january, I can use this comand with --p-min-overlap 4 to change the overlap default into 4nt.
qiime dada2 denoise-paired
--i-demultiplexed-seqs demux-paired-end.qza
--p-trunc-len-f 150
--p-trunc-len-r 150
--p-min-overlap 4
--o-table table.qza
--o-representative-sequences rep-seqs.qza
--o-denoising-stats denoising-stats.qza
Now, I can´t use again, I have this error message:
There was a problem with the command:
(1/1?) no such option: --p-min-overlap
Why can I use this before, and not now? Is there a new tool to perform this analysis?
jwdebelius
(Justine Debelius)
December 7, 2023, 2:38pm
2
Hi @c_malossi ,
I'm assuming you're doing 515F-806R 2x150 sequencing?
It's pretty widely accepted to just use forward reads for this purpose, so no need to modify the trim length. Here's a more detailed forum post on the issue:
Hi all,
Running into a weird problem I’m helping to trouble shoot. We have a run which may have some technical issues, which we are re-running, but wanted to see if we could use only the forward reads to de-multiplex and annotate.
Essentially, the forward and reverse runs are in fastq files and the barcode is in another fastq files. I realized that since the barcode is not in the sequence itself it is not useful if I use cutadapt. The forward read/reverse read fastq/barcode fastq have been per…
Best,
Justine
2 Likes
system
(system)
Closed
January 7, 2024, 8:39pm
3
This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.