ERROR RNA SEQ QC analysis step
1
0
Entering edit mode
8.2 years ago
David_emir ▴ 460

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 -o SRR1604991.txt-Q33  But I got the following Error: fastx_quality_stats: Invalid quality score value (char '#' ord 35 quality value -29) on line 4  I used -Q33 as well but failed to get any resolution!!! Guys, please help me out, I am so new to this. Thanks a ton ... A few lines of my FASTQ file for reference: [Ateeq@BIO-DT-415 NEW_NGS]$ head SRR1604991.fastq
NTCCTACTCTTCTTAGCGCCTACCCTCATACCTATCTCCCTCCTCCCATCTCCTAGGGGACTGGCGCCAAATGGTCTCTCCCTGCCAATTTTGGTATNTTCCGGGACCAAAATAAAGAGCAAGCAGGCCCCCTTCACTGAGGTGCTGGGTAGGGCTCAGTGCCACATTACTGTGCTTTGAGAAAGAGGAAGGGGATTNGT
#1=DFFFFHHHHHJJJJJJJJJJJJJJJIJJJJJJJIIJJIJJHIJJJIJGIJIJJJJJIIHHHHFFDDDDDDDCCDDDDDDDDDDDDDDDDDD<C####CCCFFFFFHHHHHJJJJIJJJJJJJJJJJJJJHIJIJIJIII=FGGIIJ=FHIIJEHHH;@DDDDEECEEEE>CDDDDDDDDDDDCCDDDDDDDBDD###
NTGAGTAGCACTCTCTGAGAGCTCCAATTTCATCCGTCTGCCATCGGCGCCATCCTGCAATCTAAGCCACAATGGTGCGCATGAATGTCCTGGCAGATGCCTCTTTTCGGCATTGTTGATACTCTTGAGAGCATCTGCCAGGACATTCATGCGCACCATTGTGGCTTAGATTGCAGGATGGCGCCGATGGCAGACGGATG
#4=DDDDFHHHHHHJJJJJJJJJJJJJJJJJJJJJJHHHIIJJJJJJIGIJIIIJHHHHHGFDFFFFEEE;65;>;@CDDDDDDDEDDEEDDDDDDDDDDCCCFFFFFHHHHHJJJJJJJJJIJJJJJJJJJJJJHIIJJJJJJJJJJJIIHIJGIJIJJJEHGHHFFFFFFEEEEEDDDDDDDDDDDDDDDDDDDDDDD
GCCACCACGGAGCGAGAAGCCCAGATAGACGCCCCGGCGGCCCCGGGTCCTGGAGTCCCGCCGCCTGCTGCCCGGCCGAGGACCCCACCCCGCCTGCCGCCCTGTCCATGGCGGGCCCCACTGCAAGCATCGGGCGGCAGGCGGGGTGGGGTCCTCGGCCGGGCAGCAGGCGGCGGGACTCCAGGACCCGGGGCCGCCGG
[Ateeq@BIO-DT-415 NEW_NGS]$/usr/local/bin/fastx_quality_stats -i SRR1604991.fastq -o SRR1604991.txt-Q33 fastx_quality_stats: Invalid quality score value (char '#' ord 35 quality value -29) on line 4  -Ateeq Khaliq RNA-Seq • 2.0k views ADD COMMENT 1 Entering edit mode if you specify -Q33 it should work, since "#" is an accepted symbol for the phred33, as you can see from this page. I tried with your example and for me it works fine. I initially thought it was a typo, but did you check there's a space before -Q33 in your line? ADD REPLY 0 Entering edit mode I have no clue why its giving error :( its showing for 'J' also. please have a look :( [Ateeq@BIO-DT-415 NEW_NGS]$ /usr/local/bin/fastx_quality_stats -i SRR1604991.fastq -o SRR1604991.txt -Q33
fastx_quality_stats: Invalid quality score value (char 'J' ord 74 quality value 41) on line 4

Edit: I checked, space is making no difference :(

2
Entering edit mode
8.2 years ago
Martombo ★ 3.0k

ok then apparently it is using the old phred33 scores (sanger), that were defined only until 40, which is "I". while with the new illumina scheme a "J" is also possible. as I told you I tried running your command with your example reads and it's giving me no error. maybe you're using an old version of fastx? mine is 0.0.13.2. what's yours?

alternatively you could use another tool, like fastqc, which is giving a similar output

0
Entering edit mode

Thanks a lot, I resolved this, Its an version error !!! Thanks again, this helped me lot !!!