extensions/net.sf.basedb.reggie/trunk/src/net/sf/basedb/reggie/grid/scripts/readlength_averager.awk

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