This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bioawk output in original fastq format

Hi, I was trying to use bioawk to extract some matched sequences from a fastq file. However, the output is somehow tabular format of fastx with $name, $seq, $qual and $comment. Did anyone know if there is an option from bioawk which can make the output keep the original fastq 4 line format? I know I can use script to process the output to make it similar as original fastq format, but if bioawk has an option for such type output, it will be very nice. Thank you.

Xiaokuan

bioawk fastq
• 3,426 views • link
updated by Stephane Plaisance Stephane Plaisance • written by weixiaokuan weixiaokuan

1 answer 

bioawk -c fastx '{print "@"$name" "$comment"\n"$seq"\n+\n"$qual}' input.fq > output.fq

this should output a duplcate of the input, add more code to filter or modify any

• 1 views • link
updated by ATpoint ATpoint • written by Stephane Plaisance Stephane Plaisance

Please use the formatting bar to highlight code.

bar

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this