Hi, I have question related to 16s metagenomics analysis. I have a pipeline including:
USEARCH (dereplicate, cluster and map to OTU) -> RDP classifier
I have a sample with 200k joined pair-end seq from illumina that I know at least 90% are valid bacterial 16s.
But after dereplicate and cluster I got only ~6k sequences which are classified with RDP at the end file from which I can calculate bacterial percentage for the taxonomy. But these 6k sequences are I think the non-redundant sequences and their copy number are not taking into account?
If I run just RDP classifier directly with the sample FASTA file I got all sequences classified which sound better..... So I wonder what is doing this USEARCH and isnit nessesary step?
•
written
by
manekineko
1 answer
I think you are asking your question wrong or you dont look good enough to your own data.
You say that you start with 200k merged reads and after dereplication and clustering you can only identify 6k read.
My question to you, what do you think dereplication and clustering does?
Log in to answer this question.
-
very few sequences classified with Kraken2
written by 10mz1 •
-
Database Recommendations for 16S-seq of the V4 Region
written by dmj6ab •
-
16S rRNA extraction of assembled genome bins
written by luyang1005 •
-
Make taxonomy bar plot using qiime2 and error in denoise procedure and training classifier
written by pple2202 •
-
Get all bacterial protein sequences from Genbank but non redundant
written by mschmid
-
Features that distinguish phyla in 16S sequences
written by grp2009 •
-
Clustering of metabarcoding reads from many environmental samples
written by lvogel •
-
16s amplicon metagenomics - no mapping file?
written by manekineko
-
Why MSA in MOTHUR reference sequences?
written by Lynxoid
-
How To Manually Check Your Putative Chimeric Sequences
written by samlambrechts299
