Question: Is it necessery to perform GATK-BQSR before performing variant calling using freebayes.
0
gravatar for sktbanerjee1
3 months ago by
sktbanerjee120
sktbanerjee120 wrote:

I am trying to call variants from whole exome paired end sequences aligned to HG38p7 using the Freebayes variant caller. As it is a general recommendation to perform BQSR before any type of variant calling, I was wondering is it necessary to perform the BQSR before calling variants with freebayes?

freebayes bqsr gatk • 188 views
ADD COMMENTlink modified 3 months ago by finswimmer77190 • written 3 months ago by sktbanerjee120
1
gravatar for finswimmer77
3 months ago by
finswimmer77190
finswimmer77190 wrote:

Hello,

the authors of freebayes itself are saying, that neither BQSR nor VQSR is neccessary.

fin swimmer

ADD COMMENTlink written 3 months ago by finswimmer77190

Could you tell me where can I find the original that "the authors of freebayes itself are saying, that neither BQSR nor VQSR is neccessary." I want to find out why BQSR is unnecessary

ADD REPLYlink written 7 weeks ago by xuzhougeng0

Hello, just follow the link in my post and you will find it. A copy and paste from the website:

FreeBayes incorporates a number of features in order to reduce the complexity of variant detection for researchers and developers:

Indel realignment is accomplished internally using a read-independent method, and issues resulting from discordant alignments are dramatically reducedy through the direct detection of haplotypes.

The need for base quality recalibration is avoided through the direct detection of haplotypes. Sequencing platform errors tend to cluster (e.g. at the ends of reads), and generate unique, non-repeating haplotypes at a given locus.

Variant quality recalibration is avoided by incorporating a number of metrics, such as read placement bias and allele balance, directly into the Bayesian model. (Our upcoming publication will discuss this in more detail.)

ADD REPLYlink written 7 weeks ago by finswimmer77190
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: 415 users visited in the last hour