This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Trouble with transcriptome data

Hi everyone

I have 14 Gb of transcriptome data of a plant sequences using Ion Proton. Initially I ran a quality assessment tool for checking it and there are a lot of overrepresented sequence probably rRNA. When I try to use fastq_quality_filter to remove low quality bases it gives me an error "Invalid quality score on line 210650840 [quality tok >:::;998;;::;5;;5<<18:;<<<=?7<;;;;;;<<==?;;;5;;299;::"

I tried searching for similar posts. Now m thinking whether to remove the whole read or not. Please suggest what can be done

rna-seq next-gen ion proton
• 2,213 views • link
updated by Benn Benn • written by gayachit gayachit

Please chose a more informative title and add relevant tags, such as ion proton

• 0 views • link
written by WouterDeCoster WouterDeCoster

2 answers

Hi gayachit,

This is a naive answer. But here it goes :

It will not affect your data if remove a single read.

and generally the first few bases of most sequencing reads have poor qualities and they need to be chopped off.

so, you can try chopping of first few bases.

• 0 views • link
written by gangireddy gangireddy

Thanks gangreddy Those bases were at the end of my file. So i removed them

• 0 views • link
written by gayachit gayachit

You'll probably have to set the phred score parameter to -Q33.

• 0 views • link
written by Benn Benn

-Q33 didn't work. I tried that earlier... But thanks anyway

• 0 views • link
written by gayachit gayachit

Okay, it would be helpful if you explained already what you have tried. Or present some code, just an idea.

• 0 views • link
written by Benn Benn

Log in to answer this question.

More posts like this