This is a test version of Biostars. For the public version, visit https://www.biostars.org.
tophat: qual length error after trimming

I did trimming with Trimmomatics v0.36 with parameter SLIDINGWINDOW:10:20 MINLEN:60, and I ran tophat and got this error:

prep_reads v2.1.1 (13b3420)
---------------------------

gzip: /lustre/beagle2/xyang/NB/Target/RNAseq/FastQ_trim/SRR1559168_1.fastq.gz: unexpected end of file
Error: qual length (13) differs from seq length (101) for fastq record SRR1559168.54082770!

gzip: stdout: Broken pipe

Is there a tool to fix this?

rna-seq rna-seq software error
• 1,919 views • link
updated by GenoMax GenoMax • written by zz105 • zz105

Looks like you ended up with a corrupted file. There probably was an error generated during the trimmomatic run. Go back and look for it. Best option would be to redo the trimming. If your original sequence file is corrupt then you would need to redownload it.

• 0 views • link
written by GenoMax GenoMax

Thanks again, this is another file that passed the fastqc step, so there is definitely something wrong in the trim process

• 0 views • link
written by zz105 • zz105

0 answers

No answers yet.

Log in to answer this question.

More posts like this