DADA2: loss too many reads after filtering

[Xinming_Xu]

Hi, timanix!

Thanks for your suggestions! I have two questions now after I read through some other posts.

1. So based on the below post, it explains --p-trunc-len discard sequence greater than the value it gives, which in my case, it discard sequence longer than 244 instead of discard sequence shorter than 244? If I understand correctly? I got confused by this post and your reply, which is opposite answer to me.

DADA2 truncation

Got it, thanks! So for example if my sequences are 150 bases long but I set the value for --p-trunc-len as 151, then nothing will be truncated and my sequence will remain as 150 bases?

@ek_97

Wouldn't the --p-trunc-len discard all the sequences greater than this value since this command would be truncating the right side? So if I set the value to 150, wouldn't that get rid of all the bases above the 150 position, and not below?

You're exactly correct, good catch! I'll make sure we get that documentation updated.

2. I actually run it again with lower truncation values, like

--p-trunc-len-f 230
--p-trunc-len-r 240 \

and

--p-trunc-len-f 240
--p-trunc-len-r 240 \

All get worse results, and with the lower value I give to the forward reads, I only got less than 10% reads shown as non-chimeric.

In case you also need to check the stats file, I will upload it here.
denoising-stats_230_240.qzv (1.2 MB) denoising-stats_240_240.qzv (1.2 MB) denoising-stats_no_truncation.qzv (1.2 MB)

What other parameters I can set to get rid of low-quality reads? Also, what causes the worse results after I lower truncation values since it should be better after I truncate the bad quality reads?

Again, really appreciate your help!