Transipedia / dekupl-run

Identify differentially expressed k-mers between RNA-Seq datasets
MIT License
11 stars 11 forks source link

Bug in reverse complementation of FASTQ files #73

Open dcansugurer opened 6 months ago

dcansugurer commented 6 months ago

Hello,

In revCompFastq.pl , there is a bug when reversing the quality score line of the FASTQ file. "Else if" condition in line 39 is:

elsif($i % 3 == 3)

But something modulo 3 cannot be equal to 3 so this condition is simply skipped everytime, give no errors and quality scores cannot be reversed while read sequence is reverse complemented.

The problem is solved when we modify it as:

elsif($i % 4 == 3)

Best wishes. Cansu.