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