From mpileup on bam files to TMB calculation, how is it done?
0
0
Entering edit mode
2.7 years ago
Pin.Bioinf ▴ 300

Hello, I want to calculate the Tumor mutational burden of my responders and non responders to treatment. I have many question regarding this matter:

  1. Should I calculate a TMB for responders and another TMB for non responders or only a global TMB for all samples?

  2. samtools mpileup , bcftools norm and bcftools view have to be applied to all bams together or to each bam of each sample? Right now I am doing it for all bams at the same time to get one vcf. Then, I would annotate the vcf file and count the number of snps that are non synonymous, and divide by region size to get the TMB.

  3. I am doing this with RNA-Seq data for orientative purpose and practice, I know it is not the ideal situation. But, what is the size of the coding region of the targeted territory by which I have to divide my mutation number in RNA-Seq? Is it the length of the reads ( 100b) or something like 3099750718 bp (got it from a coverage reference file when mapping), or 30 (length used when TMB is calculated for whole exon sequencing)?

Thank you so much for your help.

mutations TMB snps bam mpileup • 733 views
ADD COMMENT

Login before adding your answer.

Traffic: 3274 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6