Hello all!
I'm new to the community and just started my first bioinformatics project. I'm filtering out merged reads with the expected error of 1-100.
code I'm running:
vsearch -fastq_filter ./MergedFiles/${otp}_merged.fastq -fastq_maxee 1 -fastqout ./MergedFiles/${otp}_mergedfilt.fastq -fastaout ./MergedFiles/${otp}_mergedfilt.fasta
The issue is that when I run it, I get back a fatal error saying-
Reading input file 0%
Fatal error: FASTQ quality value (42) above qmax (41)
I've tried usearch, but my file is too big for the free 32-bit version, and I've tried replacing fastq_maxee 1 with fastq_qmax 42, but that doesn't filter anything out. I'm a bit new to all this and I'm stuck. Any help would be greatly appreciated!
-science panda
I'm new to the community and just started my first bioinformatics project. I'm filtering out merged reads with the expected error of 1-100.
code I'm running:
vsearch -fastq_filter ./MergedFiles/${otp}_merged.fastq -fastq_maxee 1 -fastqout ./MergedFiles/${otp}_mergedfilt.fastq -fastaout ./MergedFiles/${otp}_mergedfilt.fasta
The issue is that when I run it, I get back a fatal error saying-
Reading input file 0%
Fatal error: FASTQ quality value (42) above qmax (41)
I've tried usearch, but my file is too big for the free 32-bit version, and I've tried replacing fastq_maxee 1 with fastq_qmax 42, but that doesn't filter anything out. I'm a bit new to all this and I'm stuck. Any help would be greatly appreciated!
-science panda