PRINSEQ bug with -no_qual_header

26 views
Skip to first unread message

Emmanuelle Beyne

unread,
Jun 2, 2021, 4:25:40 AM6/2/21
to Edwards Lab Tools

Hello,

I have several FASTQ files (from Illumina sequencing) with reads as the following :
@ST-E00181:926:HFWF7CCX2:7:1108:17939:71770 1:N:0:TAGGCATG+CTCCTTAC
ACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACT
+
AAFFFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJAJJJJJJJJFJJJJJJJJJJJJJJJJJJJJJJJJJJJAFFJJJJJJJJFJJJJJJFAF7
@ST-E00181:926:HFWF7CCX2:7:1108:17959:71770 1:N:0:TAGGCATG+CTCCTTAC
...

If I specify the '-no_qual_header' option in my Prinseq command line, we obtain for this read:
@ST-E00181:926:HFWF7CCX2:7:1108:17939:71770 1:N:0:TAGGCATG+CTCCTTAC
ACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACT
+
AAFFFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJFJJJJJJJJJJJJJJ+
AAFFFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJ
@ST-E00181:926:HFWF7CCX2:7:1224:26230:11347 1:N:0:TAGGCATG+CTCCTTAC
...

with a supplementary and false line. And there are many errors like this one...
No problem if I don't use this option (but files have a bigger size).
I know the Prinseq  software is a "old" tool but some biologists like it (and use it )...

Best regards,
Emmanuelle


Reply all
Reply to author
Forward
0 new messages