Tumor mutation burden (TMB) and microsatellite instability (MSI) score
First, dysfunctional mutations in exonic region and at splicing junction sites with variant frequency not smaller than 0.05 were counted. Then, the TMB value was calculated as the panel size normalized number of counted mutations. TMB was determined by dividing the total number of identified mutations by the panel size (2MB). MSI score was calculated using MSIsensor-pro (v1.0.2) based on BAM file recording the aligned sequencing reads