Question: copmaring mpileup parameters: -B vs -BQ0
0
gravatar for -_-
6.1 years ago by
-_-880
Canada
-_-880 wrote:

Hi there, I have seen the following difference which I cannot explain with samtools mpileup.

First, only -B is specified to disable BAQ calculation, it finds two reads with insertions. I am not sure if this matters, but if you look closely, the two inserted sequences are NOT exactly the same.

$ samtools mpileup -B -f ${REF} ${BAM_FILE} | grep 175200931
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-file depth to 8000
chr1    175200931    a    20    ....,...+30CATGAATATATACACACGTATATATACATA..+30CATGTATATATACACACGTATATATACATA..........    Jm_sEHIDFsEHIdJ>JLJD

Second, I rerun the experiment with -BQ0 specified, a third insertion appears.

$ samtools mpileup -BQ0 -f ${REF} ${BAM_FILE} | grep 175200931

[mpileup] 1 samples in 1 input files
<mpileup> Set max per-file depth to 8000
chr1    175200931    a    27    ...,.,...+30CATGAATATATACACACGTATATATACATA..+30CATGTATATATACACACGTATATATACATA...,+30catgtatatatacacacgtatatatacata......,,,,..    Jm_!sEHIDFsEHI!dJ>JLJ!!!!,D

From my understanding of -Q0, it means skip bases with BAQ smaller than 0. In other words, since BAQ cannot be samller than 0, it keep all reads whatever their BAQ are.

However, I don't understand if -B is speicified to disable BAQ calculation already, why would -Q0 make a difference? Where comes the BAQ data then?

mpileup samtools baq • 1.4k views
ADD COMMENTlink modified 6.1 years ago by Ram32k • written 6.1 years ago by -_-880
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: 1555 users visited in the last hour
_