Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

when to apply trimming parameters #52

Open
smzt opened this issue Jan 17, 2023 · 1 comment
Open

when to apply trimming parameters #52

smzt opened this issue Jan 17, 2023 · 1 comment

Comments

@smzt
Copy link

smzt commented Jan 17, 2023

Hi,
I've been going through the documentation and the paper but I have not been able to find a clear answer yet to my naive question.
When should I apply the values Figaro provides? Should I apply them before trimming the primers or after trimming the primers?

Thanks in advance.

Regards,

Sheila

@alegione
Copy link

My understanding if helpful:
The tool is written for Dada2, so the output results should be used as truncation lengths in dada2.
The tool truncates reads first, then removes the primers if you ask it to also trim the ends.

Example if using dada2 in R (all caps text to be replaced with your data)

fastqPairedFilter(fn = c(FORWARD_IN, REVERSE_IN), fout = c(FORWARD_OUT, REVERSE_OUT), truncLen = c(FIGARO TRIM POSITION 1, FIGARO TRIM POSITION 2), verbose = TRUE, maxN = 0, maxEE = 2, truncQ = 2, trimLeft = c(FORWARD PRIMER LENGTH, REVERSE PRIMER LENGTH), compress = TRUE)

Example if using dada2 in QIIME2

qiime dada2 denoise-paired \
    --i-demultiplexed-seqs demux-paired.qza \
    --p-trunc-len-f FIGARO TRIM POSITION 1 \
    --p-trunc-len-r FIGARO TRIM POSITION 2 \
    --p-trim-left-f FORWARD PRIMER LENGTH \
    --p-trim-left-r REVERSE PRIMER LENGTH \
    --o-representative-sequences representative-sequences.qza \
    --o-table table.qza \
    --o-denoising-stats denoising-stats.qza

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants