My apologies in advance - I am new to bioinformatics and QIIME2 (v2023.2).
I am working on a project where I am trying to compare and correlate microbiomes from highly enriched regions such as feces/intestines to regions with much low levels of the microbiome (think lungs, certain tumors, etc.). I performed MiSeq 2x250bp for the V4 region.
I performed primer trims and truncated low quality regions of the sequences, but I am getting super low frequencies.
I performed DADA2 denoising on q2galaxy server because of a slow computer, but my truncation length was F: 230 and R: 160.
What would be your suggestions for this sort of data analyses? A lot of my samples of interest (H1, H2, H3, H4) have 0 frequencies. Any help would be appreciated. Thanks!
There is a lot going on here, but the quickest bang for your buck will be to adjust your forward truncation length lower. From the interactive quality plot tab you can see that the median forward read length is 230 nts. Setting your forward truncation length to 230 thus discards 50% of your reads right off the bat. Try 225 or so.
It also looks like you are safe to extend the reverse truncation cutoff to about 220. Try changing these parameters and rerunning, and we can discuss from there.
These look like solid results to me, you're safe to move forward with these outputs.
Please note that, unfortunately, around half of your samples got sequenced with little depth. There is nothing you can do about this (except re-sequence), but don't expect great things from those analysis wise.