Alternative content
#look at GC across the entire CDS. gffread -x - -g <genome.fasta> <annotation.gff> | \ seqtk comp - | \ awk -v OFS="\t" '{ print $1, "0", $2, ($4 + $5) / $2 }'