Extracting longest read from fastq file

Question

[KVC_bioinfo]

Hello All,

I want to extract the longest read from a fastq file. I am not aware of any way to do it. Could someone help me here?

TIA

[h.mon]

Try searching for Extracting longest read from fastq file. I've found some good results.

Answer 1

[Pierre Lindenbaum]

gunzip -c in.fq.gz | paste - - - - | awk -F '\t' '{L=length($2);if(L>M) {M=L;R=$0;}} END {print R;}' | tr "\t" "\n"