Question: From mpileup on bam files to TMB calculation, how is it done?
0
gravatar for Pin.Bioinf
8 weeks ago by
Pin.Bioinf240
Malaga
Pin.Bioinf240 wrote:

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.

tmb mpileup mutations bam snps • 115 views
ADD COMMENTlink written 8 weeks ago by Pin.Bioinf240
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1062 users visited in the last hour