Convert Phred score 41 to 40
2
0
Entering edit mode
7.5 years ago

I need to convert new data (Illumina 1.8+) to 454 format (old Sanger).

According to this page (http://www.dylanstorey.com/node/54), only difference in these two formats is that new Illumina can have score 41 (letter J), whereas old Sanger goes only up to 40. Is there any tool to convert all 41 scores to 40 (letter I) in fastq?

fastq • 2.0k views
ADD COMMENT
0
Entering edit mode

I think this question has been asked before: Convert Illumina Reads To Sanger Score

ADD REPLY
0
Entering edit mode

Its not the same question as I am not interested in converting ASCII-64 to Sanger score (ASCII-33). Instead, both my formats are in ASCII-33 (Sanger), just one goes from 0..41 and second from 0..40

ADD REPLY
2
Entering edit mode
7.5 years ago
xb ▴ 420

Is this what you want?

cat <your-fastq-file> | awk ' {if(NR %4 == 0) {gsub("J","I");print} else {print}}'

If it is a .gz file, use zcat to read it and gzip to zip it back.

ADD COMMENT
1
Entering edit mode
7.5 years ago

If you don't care about converting all the 'J's to 'I's (slightly lowering the scores of some bases), you can use sed:

sed -e '4~4y/J/I/' fastqfile.fastq > fastqfile.illuminaphred.fastq

You can also add the -i option to modify the file in place, without having to copy it to fastqfile.illuminaphred.fastq.

ADD COMMENT

Login before adding your answer.

Traffic: 1615 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