计算TMB

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