下載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}'
計算CDS長度
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點，簡書系信息發(fā)布平臺，僅提供信息存儲服務(wù)。
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點，簡書系信息發(fā)布平臺，僅提供信息存儲服務(wù)。
相關(guān)閱讀更多精彩內(nèi)容
- 生信學習筆記 linux部分功能 查看文件夾 工具 選項 可以設(shè)置鼠標功能 可以設(shè)置右鍵粘貼 雙擊這個窗口可以再打...
- 要求 實現(xiàn)這個功能的軟件也很多，還是煩請大家先自己搜索幾個教程，入門請統(tǒng)一用htseq-count，對每個樣本都會...
- 原文鏈接：HOPTOP轉(zhuǎn)錄組入門(6)： reads計數(shù)-轉(zhuǎn)錄組-生信技能樹 要求 實現(xiàn)這個功能的軟件也很多，還是...
