BQSR using GATK without dbsnp
Hi, I have run the base quality score recalibration (BQSR) tool from GATK but for the known sites parameter I used only two databases:
Mills_and_1000g_gold_standard_indels1000G_phase1.indels
I realise I should possibly also have included the dbsnp_138 database in this list. Is it necessary to re-run with all three databases?
1 answer
Log in to answer this question.
More posts like this
-
Genome mapping mm39 via GTK (BaseRecalibration) - Where to get known sites of polymorphisms?
written by Rashid •
-
Help with Base quality Score Recalibration (BQSR) using GATK
written by sainavyav22 •
-
Apply BSQR for Targeted Sequencing
written by Manuel Sokolov Ravasqueira •
-
dbSNP and indels
written by wonde2000 •
-
Base recalibration in normal vs. tumor somatic variant calling in WXS data?
written by rebeliscu •
-
So many variants detected.
written by Mehmet
-
GATK quality recalibration to get recalibrated BAM files before SNP calling
written by alicia.poppy
-
BQSR vs VQSR
written by James Reeve
-
BQSR on one sample vs. all samples in a run
written by ari.nazarian •
-
Advice On Pre-Processing/Filtering Ngs Data
written by DonJoe •
Ideally, it would be better if you can re-run the BQSR analysis using both SNPs and Indels. But it may or may not affect the downstream variant results depending on how good your sequencing data is. Honestly speaking, I have never found a big difference in the number and quality of variants called before or after BQSR but it could be solely due to genuinely good quality of reads that we get from our sequencing core. The other steps such as indel-realignment do have a big effect on the quality of variants.