other/pipeline/trunk/readlength_averager.awk

Code
Comments
Other
Rev Date Author Line
6419 23 Sep 21 nicklas 1 #!/bin/awk -f
6419 23 Sep 21 nicklas 2
6419 23 Sep 21 nicklas 3 # Nicklas Nordborg, September 2021
6419 23 Sep 21 nicklas 4 # Calculates average read length in FASTQ files
6419 23 Sep 21 nicklas 5
6419 23 Sep 21 nicklas 6 BEGIN {
6419 23 Sep 21 nicklas 7     mean=0;
6419 23 Sep 21 nicklas 8     n=0;
6419 23 Sep 21 nicklas 9 }
6419 23 Sep 21 nicklas 10
6419 23 Sep 21 nicklas 11 {
6419 23 Sep 21 nicklas 12     if (NR%4==2) {
6419 23 Sep 21 nicklas 13       n++;
6419 23 Sep 21 nicklas 14       mean+=length
6419 23 Sep 21 nicklas 15     }
6419 23 Sep 21 nicklas 16 }
6419 23 Sep 21 nicklas 17
6419 23 Sep 21 nicklas 18 END {
6419 23 Sep 21 nicklas 19     mean/=n;
6419 23 Sep 21 nicklas 20     printf("%ld\t%g\n",n,mean);
6419 23 Sep 21 nicklas 21 }