usearch -filter_phix command 信息搬

信息来源：https://www.drive5.com/usearch/manual/cmd_filter_phix.html

Filter reads for hits to the PhiX genome which is used as a spike in Illumina amplicon sequencing.

Input is a FASTA or FASTQ file. A reverse (R2) read file can be specified using -reverse, in which case a read pair is discarded if the forward (R1) and/or reverse (R2) read has a hit to PhiX. If -reverse is specified, -output2 must be specified for the reverse read output file. Output is written in the same format (FASTA or FASTQ) as the input.

Standard output files are supported for reporting PhiX hits.

Multithreading is supported.

Examples

usearch -filter_phix reads.fq -output filtered_reads.fq -alnout hits.txt

usearch -filter_phix reads_fwd.fq -reverse reads_rev.fq -output fil_fwd.fq -output2 fil_rev.fq

search_phix command

Search sequences for hits to the PhiX genome which is used as a spike in Illumina amplicon sequencing. This command works exactly like ublast except that no database file is needed (no -db option) because the database is compiled into usearch.

Use the filter_phix command to filter paired reads.

Example

usearch -search_phix reads.fq -notmatchedfq filtered_reads.fq -alnout hits.txt