BEGIN{FS=OFS="\t"; getline; getline; ns=NF; for (i=1;i<=ns;i++){ sampID[i]=$i; } }{ for (i=7;i<=ns;i++){ x=$i; nr[i]+=x;fl[NR]=$6; r[NR,i]=x;l[NR]=$0; } }END{ nl=NR; print "Sample","Total_mapped_reads" for (i=7;i<=ns;i++){ print sampID[i],int(nr[i]) } }