Quality-filtering for 16S Microbiome Taxonomic Profiling

For 16S microbiome taxonomic profiling, the following criteria are used to filter out sequencing reads with low quality:

  • Sequences with the lengths of <80 bp or >2,000 bp
  • Averaged Q value is <25.
  • Not predicted as a 16S gene by the Hidden Markov Model (HMM) based search.
  • Sequences are first assigned to the reference 16S database. All sequences that do not match any of reference sequences with at least 97% similarity cutoff are clustered using UCLUST method using 97% the cutoff. If a sequence is found to be a singleton, we assume that it is an erroneous one that should be excluded in the subsequent analyses. This algorithm is widely used, especially for Illumina short read sequencing [See QIIME manual’s step 5].

 

 


The BIOiPLUG team / Last edited on Feb. 19, 2018