Question: Expected error filtering on Ion Torrent data using fastq_filter USEARCH command
0
gravatar for sukhu.sakhare
9 months ago by
sukhu.sakhare0 wrote:

Hello, I am doing 16s rRNA analysis on Ion Torrent data using Usearch pipeline. I got error at quality filtering step. Command used for the quality filtering is:- usearch -fastq_filter Reads_trunc.fq -fastq_maxee 1.0 -fastaout Reads_filtered.fa

Error got is:- "CharToIntQual('L') Phred score 43 out of range 0..41"

I have changed the -fastq_maxee value to 0.5 but still got the same error. I can't figure out the reason behind the error. Need help

Thanks in advance.

metagenomics • 490 views
ADD COMMENTlink modified 5 weeks ago by Simply Bioinformatics80 • written 9 months ago by sukhu.sakhare0
0
gravatar for Brian Bushnell
9 months ago by
Walnut Creek, USA
Brian Bushnell15k wrote:

The problem is that the quality scores are out of the expected range. Before you can proceed further, you need to figure out why that is the case. For example, they may have been incorrectly processed, the file may be corrupt, or it may be that they are actually not Ion Torrent reads but very old Illumina reads. I am not familiar with Ion Torrent data so it could be that 43 is a normal quality score for it, in which case you might need to restrict the range to 0-41, which you can do with BBMap's reformat.sh. But first find the reason why you have a base with Q43.

ADD COMMENTlink written 9 months ago by Brian Bushnell15k

Thank you for your suggestion. I have used -fastq_qmax parameter and it worked. Thank you

ADD REPLYlink written 9 months ago by sukhu.sakhare0
0
gravatar for Simply Bioinformatics
5 weeks ago by
WashingtonDC
Simply Bioinformatics80 wrote:

USEARCH Fastq-Filter for Ion-Torrent 16S
usearch -fastq_filter in.fastq -fastqout out.fastq -fastq_ascii 33 -fastq_qmin 0 -fastq_qmax 47

ADD COMMENTlink modified 5 weeks ago • written 5 weeks ago by Simply Bioinformatics80
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 617 users visited in the last hour