Question: bcftools variant calling query
gravatar for f1978hp
9 months ago by
f1978hp30 wrote:

"During variant calling bcftools compares the average mapping quality of reads supporting a variant with those they don't. If there is a significant difference it will skip this variant. " Are there someone that can explain me how to create a BAM file with the reads that didn't supported the variant calling process in this case?

ADD COMMENTlink modified 9 months ago by Kevin Blighe59k • written 9 months ago by f1978hp30
gravatar for Kevin Blighe
9 months ago by
Kevin Blighe59k
Kevin Blighe59k wrote:

Can you provide the source of the statement?

Anyway, not a direct answer but you can negate this filtering by just filtering out all reads below a certain MAPQ threshold, e.g., MAPQ 60:

samtools view -bq 60 Aligned.bam > Aligned_FiltMAPQ.bam

Then, there cannot be a statistically significant difference (based on MAPQ) between those reads that support a call versus those that don't.


ADD COMMENTlink written 9 months ago by Kevin Blighe59k

Hi Kevin!

Thanks for your interesting!

The source of the statement is another answer to my query that I was propose later. See the link Command line explain There you will undertand the problem better.

ADD REPLYlink written 9 months ago by f1978hp30
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1513 users visited in the last hour