Base recalibration in normal vs. tumor somatic variant calling in WXS data?
I have a tumor and a normal BAM file and am preparing to run base recalibration.
I was planning on calling variants on the normal and using that, in addition to dbSNP, as input for recalibration of tumor BAM(s), e.g.:
gatk BaseRecalibrator -I tumor.bam -R hg38.fasta --known-sites normal.vcf --known-sites dbSNP_hg38.vcf -O tumor_recal.table
Before producing the normal VCF however, it’s not clear to me whether I should run base recalibration on the normal BAM. If this is advised, I had planned using dbSNP as the known sites (for normal), e.g.:
gatk BaseRecalibrator -I normal.bam -R hg38.fasta --known-sites dbSNP_hg38.vcf -O normal_recal.table
Alternatively, I could keep things simple and run base recalibration on both tumor and normal using dbSNP only.
Is one of these workflows more preferable?
Any clarity here would be much appreciated. Thanks!
• 20 views
Read more here: Source link