This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!
Base Quality Recalibration With Low Frequency Variants
Hi I'm working on calling low frequency heteroplasmic variants in the mitochondrial genome, using WGS data. My data is high coverage (~2000x), and heteroplasmic variants might be represented on as few as 10 reads. Am I correct that BaseRecalibrator will treat low frequency haplotypes as sequencing error in this case and reduce Phred scores for these reads? If so, is BQSR not appropriate for this data?