FASTX-toolkit in iMAC

Hi, When I want to use FASTX-toolkit on iMAC , I get the following error "Invalid quality Score value (char '#' ord 35 quality value -29) on line 4 "

what does mean? Thank you!

rna-seq

• 1,900 views • link

updated 9.8 years ago by Michael Dondrup

• written 9.8 years ago by Arash •

1 answer

I think this means that the encoding of the quality values in the fastq file is different from what the fastx tool expects. Presence of # indicates 'Sanger' encoding (Ascii code - 33):

perl -e 'print ord("#") - 33, "\n"'
2

See : http://seqanswers.com/forums/showthread.php?t=7596 for a solution.

• 0 views • link

written 9.8 years ago by Michael Dondrup

More posts like this

Fastx quality filtering error
written 7.5 years ago by athira112athira •

Hi, I use fast-x toolkit for quality filtering of my data. I get an error when running the command fastq_quality_filter -v -Q 64 -q 25 …
filtering out low quality sequences from Sanger Sequencing
written 8.3 years ago by hassani.medamine •

Hello! I am currently facing a challenge of filtering out low quality sequences (output from Sanger sequecing, roughly I have 300 sequences, command line are …
Producing the Reverse-complement of each sequence in fastq files
written 9.4 years ago by nakanomasayuki265 •

I want to produce the reverse-complement of each sequence in fastq files. I tried fastx_reverse_complement in FASTX-Toolkit. But following an error message were obtained. fastx_reverse_complement: …
bowtie "Saw ASCII character -54 but expected 33-based Phred qual." after -Q 33 fastx reverse comple…
written 9.9 years ago by gufernandez10 •

hi, I'm working with RNA-seq data. I reverse complemented sequences in a fastq file using fastx reverse complement function adding the -Q 33 option. Without …
phred 33 from Ion S5?
written 10.3 years ago by Benn

Hello, I am using Ion S5 data from an amplicon panel, and with fastq_quality_trimmer I try to trim. However I get errors that my data …
ERROR RNA SEQ QC analysis step
written 11.7 years ago by David_emir

Hello all, I was trying to do QC with `fastx_quality_stats` to calculate out nucleotide distribution. The code is as follows [Ateeq@BIO-DT-415 NEW_NGS]$ /usr/local/bin/fastx_quality_stats -i SRR1604991.fastq …
Fastx Invalid Quality Score Value
written 12.7 years ago by sckinta

I have multiple RNAseq libraries to parse, so I wrote a pipeline using bash and submitted them in batch. Most of libraries ran well and …
What Does Windown Size Exact Mean When Trimming
written 12.9 years ago by Tonyzeng •

HI I have a question, when I use SLIDINGWINDOW option in trimmomatic, what does exact mean for the windownsize, like SLIDINGWINDOW:4:15, what is 4 mean? …
Possible Bug In Fastx Toolkit
written 13.4 years ago by Abhi

I am getting the following error with fastx_quality_stats with version 0.0.13.2. May be this has to do with the changing Illumina scores. Now they are …
Fastq Quality Check
written 14.8 years ago by toshnam

Hi all, I'm trying to check sequencing quality of FASTQ file from HiSeq2000. I used fastx_quality_stats script of FASTX-Toolkit (Version 0.0.13) for it. However I've …