Question: Is it necessery to perform GATK-BQSR before performing variant calling using freebayes.
0
gravatar for sktbanerjee1
14 months ago by
sktbanerjee130
sktbanerjee130 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 • 568 views
ADD COMMENTlink modified 14 months ago by finswimmer2.1k • written 14 months ago by sktbanerjee130
2
gravatar for finswimmer
14 months ago by
finswimmer2.1k
Germany
finswimmer2.1k wrote:

Hello,

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

fin swimmer

ADD COMMENTlink written 14 months ago by finswimmer2.1k

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 13 months 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 12 months ago by finswimmer2.1k
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: 764 users visited in the last hour