Hello I am running gmap with paired ends fq.gz files. after a while I get this error message. Length 4783 of quality score differs from length 539 of nucleotides in sequence ��<���.���*D#P:/���7"A@��]�e.�m��>�rGӕ�3x��L�=�}��+|�'�nb{��5h�GXU�F�Bw/D�n�TN�Cm��
Is there any way to tell gmap to skip this sequence or should I try and find the offending line? Thank you
Do you know which line have this issue ?
If so you can extract this part from fastq file to see how it look like
assuming that the issue in line 4
sed -n '1,4p' yourFile.fastq
then you can extract it by this command
sed -n '1,4d' yourFile.fastq > newFile.fastq
BUT
If you have pair ended read you need to remove it from both files
in case if you do not know the line number you can use the awk code in this post to find it out Fastq Quality Read And Score Length Check
Thank you medhat for your reply. I do not have the line numbers, I will have to look for them as you suggested. However, it seems strange that the length of the lines is so long . I have looked at one of the fastq files and the lines there are as expected not longer than 120 chr. In the message I get its say: "Length 6925 of quality score differs from length 3471 of nucleotides".
Log in to answer this question.