Fastx_reverse_complement error. Invalid quality score even with -Q33
1
0
Entering edit mode
4.9 years ago
roblogan6 ▴ 30

I am using the fastx toolkit to get the reverse complement of fastq files that are in Sanger phred 33 format. I need to add the -Q33 flag, as discussed in this previous BioStars post Fastx Invalid Quality Score Value . However, I am getting the following error:

fastx_reverse_complement: Error: invalid quality score data on line 1464 (quality_tok = "+"

The "+" symbol on this line is a low scoring phred 33 symbol, not the "+" that separates sequence lines from symbol lines. The fastq files were generated from PacBio, which are in Sanger 33 format. I don't know why it is not being recognized or accepted. Any insights on how to overcome this? Thanks

fastq reverse-complement fastx • 2.1k views
ADD COMMENT
2
Entering edit mode
4.9 years ago

Use BBMap's reformat tool:

reformat.sh in=reads.fq out=rcomped.fq rcomp qin=33

That gets used daily with PacBio reads without problems. Fastx is a bit obsolete.

ADD COMMENT
0
Entering edit mode

Thanks Brian for pointing me in the right direction. BBMap worked very well. I noticed you are the author of BBMap too. Very cool! Thanks again.

ADD REPLY
0
Entering edit mode

You're welcome, glad I could help!

ADD REPLY

Login before adding your answer.

Traffic: 2303 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6