计算CDS长度

下载TCGA 基因注释文件：gencodev22
https://gdc.cancer.gov/about-data/data-harmonization-and-generation/gdc-reference-files

安装bedtools
conda install -c bioconda bedtools
or conda install -c bioconda/label/cf201901 bedtools

计算
awk -F "\t" '{if($3=="CDS"){print}}' gencode.v22.annotation.gtf|bedtools sort -i -|bedtools merge -s -i -|awk -F "\t" '{print $3-$2+1}'| awk '{sum += $1};END {print sum}'