Question: bcftools variant calling query
1
gravatar for f1978hp
9 months ago by
f1978hp30
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
0
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.

Kevin

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.

Help
Access

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