Skip to content

Instantly share code, notes, and snippets.

@l-modolo

l-modolo/fastq_by_size.sh

Last active August 29, 2015 13:57
Show Gist options
  • Save l-modolo/9550357 to your computer and use it in GitHub Desktop.
Save l-modolo/9550357 to your computer and use it in GitHub Desktop.
Download ZIP
Only read with more than 20bp from a fastq
Raw
fastq_by_size.sh
awk 'BEGIN {OFS = "\n"} {header = $0 ; getline seq ; getline qheader ; getline qseq ; if (length(seq) >= 20){print header, seq, qheader, qseq}}' < input.fastq > output.fastq
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment