search

lh3 / seqtk

Toolkit for processing sequences in FASTA/Q formats
MIT License
1.37k stars 308 forks source link

A bug about seqtk #66

Closed WinterLi1993 closed 8 years ago

WinterLi1993 commented 8 years ago

Hi I use seqtk to convert my fastq : seqtk seq -V -Q 64 $1.fq > $1.s.fq and Use fastqQc , but I meet such this error :
"FATAL ERROR (file fastq_qc_stat.c line 80): quality (-1) not in [0, 45]"

What happend about seqtk ? Thank you !

lh3 commented 8 years ago

Could you show the violating line before and after "-VQ64"?

WinterLi1993 commented 8 years ago

I think If I convert to 64phred ,should not appear "-29 quality" , but I use fastQc tell something wrong happend on fastq phred. seqtk seq -V -Q 64 $1.fq > $1.s.fq seqtk seq -V -Q 64 $2.fq > $2.s.fq

"fastq_qc_stat -a $1.fq $2.fq -f R02_N -q 45 -Q 64

Error Information： “FATAL ERROR (file fastq_qc_stat.c line 80): quality (-29) not in [0, 45]”

seqtk is work ??

lh3 commented 8 years ago

-VQ64 converts Q64 quality to Q33 quality, not the other way around.