This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Fastx quality filtering error

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 -p 75 -i sra_data.fastq -o sample_filtered.fastq
fastq_quality_filter: Invalid quality score value (char '.' ord 46 quality value -18) on line 4

I need a moderate threshold of average quality>25. please, can anyone suggest a possible solution for this?

software error fastx
• 2,474 views • link
updated by finswimmer finswimmer • written by athira112athira • athira112athira

Output of head -n 4 sra_data.fastq ?

• 1 views • link
written by ATpoint ATpoint

This is the output for head -n 4 sra_data.fastq

dcb@dcb-P6-2355IL:~$ head -n 4 sra_data.fastq

@SRR6031690.1.1 I2SY8KZ01A1O0T length=510  
CGCAACATTACCGCGGCTGCTGGCACGTAGTTAGCCGTCCTCTTGGTAAGATACCGTCACAGTGTGAACTTTCCACTCTCACACTCATTCTTCTCTTACAACAGAGCTTTACGATCCGAAACCTTCTTCACTCACGCGGCGTTGCTCGGTCAGACTTCCGTCCATTGCCGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCCTAGGTCGGCTATGTATCGTCGCCTTGGTGAGCCGTTACCCCACCAACTAGCTAATACAACGCAGGTCCATCTGGTAGTGATGCAATTGCCACCTTTAAGTAAATGTCATGCAACATCTACTATTATGCGGTATTAGCTATCGTTTCCAATAGTTATCCCGCCGTACCAGGCAGGTTACCTACGCGTTACTCACCCGTTCGCAACTCATCCAGAGAAGCAAGCTCCTCCTTCAGCGTTCTACTTGCATGTATTAGGCACGCCGCCAGCGTT
+SRR6031690.1.1 I2SY8KZ01A1O0T length=510
IIIIIHHIIIHHD888;FH@===IIIIIH@@@CDIIII544===>>>;;;>HADDEIIIIIIIIIIIII88888GFFFHIIIIIIIIIIIIIIIIEEEHIIIIIIHIHHAA@A@?@@@12266IHHFFHHIIIIIIIIIIIIIIIIIIIIIIIIHHHHHHF::;;;EHEGHHGHHH777CDF@@666>96111CCC@@@CE=:;;CCHFFFFCCA=;222?AFCF???>?>>>>A9;6666;666>BEEHFHEEGIIIIHIEEEEHHHHD@@@FEE@@@DFFHHHHHHHHHHEEEFFFFFFFFHHHFEEEEIIHHHEEEHDB=<;@=944465...77?=//06<BEHHHHHHHHHHHHEEEFFHFE>>>HHDDEHHHHH@:555FDCBB????72---0///000017===>>CDEDFEHA??>DDCC@@@CCCDHHHHHHHHHHHC@===?DD@AB@EDB===>?EDDEHHIDC>==DHHHECCBCCDCA@888<@<:9;555;=?..
• 1 views • link
updated by Bastien Hervé Bastien Hervé • written by athira112athira • athira112athira

Hello athira112athira,

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

Thank you!

• 1 views • link
written by Bastien Hervé Bastien Hervé

Thank you so much, I will take care about it from the very next post.

• 1 views • link
written by athira112athira • athira112athira

1 answer

Looks like a normal fastq file. In any case, I recommend going away from fastx-toolkit. It is very old and newer (and faster) alternatives exist. Options are among others: fastp, skewer, cutadapt etc. I checked with skewer and it does not throw any error. Check which tools you feel comfortable with, they all have a lot of options to customize the trimming.

• 1 views • link
written by ATpoint ATpoint

Thank you for your suggestion. I will try anyone of it as possible.

• 1 views • link
written by athira112athira • athira112athira

Log in to answer this question.

More posts like this