Originally posted by Siva
View Post
Awk is truly powerful!!
Can you modify it further to add a column that gives the number of exons per transcript?
Code:
awk '$3=="exon" && $9~/^Parent/{sub(/;.*/,"",$9);ID=substr($9,8);L[ID]+=$5-$4+1;N[ID]++}END{for(i in L){print i,L[i],N[i]}}' ZmB73_4a.53_WGS_chr1_21.gff
s.
Comment