Skip to content

Remove chimeric reads

Gavin Douglas edited this page Feb 17, 2016 · 24 revisions

Our script chimeraFilter.pl wraps usearch (v6.1), specifically the uchime algorithm, to remove chimeric reads.

Here is an example command:

chimeraFilter.pl -type 1 -db /usr/local/db/single_strand/Bacteria_RDP_trainset15_092015.udb fasta_files/*

Where "-type 1" means that any reads clearly called as chimeric AND reads that are ambiguous are filtered out. Note that a DB file needs to be input as well and if you'd like to use the UDB format rather than FASTA then you'll need to use the "-makeudb_usearch" function of usearch v6.1 (the same usearch version as used for chimera checking).

Clone this wiki locally