计算CDS长度
2019-07-28 本文已影响0人
裁尘的人儿
-
下载TCGA 基因注释文件：gencodev22
https://gdc.cancer.gov/about-data/data-harmonization-and-generation/gdc-reference-files -
安装bedtools
conda install -c bioconda bedtools
orconda 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}'
