计算TMB

作者: 裁尘的人儿 | 来源:发表于2019-07-01 22:55 被阅读0次

id=TCGA.READ.muse.ec8ec3ad-f08d-46eb-9571-42806e304b37.DR-10.0.somatic.maf.gz

zcat $id|grep -v '^#' |cut -f 16 |sort -u |grep "^TCGA" >> READ.list

cat READ.list| grep "^TCGA" | while read line; do num=`zcat $id| grep -v '^#' | grep $line | wc -l`; echo "$line $num" >> TMB_READ.txt; done

相关文章

网友评论

    本文标题:计算TMB

    本文链接:https://www.haomeiwen.com/subject/jleocctx.html