Service notice: Several of our team members are on vacation so service will be slow through at least July 13th, possibly longer depending on how much backlog accumulates during that time. This means that for a while it may take us more time than usual to answer your questions. Thank you for your patience.

base re-calibration when I don't have a publicly available dbSNP

I want to perform base re-calibration using BaseRecalibrator. However, I work on a non-model organism and there is no available SNP database.

Would it be right way if I creat my custom dbSNP file using very stringent criteria (--min_base_quality_score 30) for SNP calling and re-calibrate my data using this file?
I also wonder if I will benefit from re-calibration when I use custom dbSNP.vcf, which is produced from the same samples I analyze?

I create dbSNP.vcf using UnifiedGenotyper to save time, while for genotyping after base re-calibration I use HaplotypeCaller. I have 24 samples in this analysis.

I would appreciate your assistance.

Best Answer

Answers

Sign In or Register to comment.