Hi, there
I compare the mpileup outputs between samtools 1.2 and 0.1.18. And I find
some discordance of the base quality between the old and new version.
I extract the pileup file with -B option to disable BAQ of samtools firstly. And I get a position like this:
gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA BBBB?@@ABBB@=@?C>ABBCB;
Then I perform the same command except the -B option, but for 0.1.18, I get this:
gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA <<<<>>:=;<<5=@;:>==;::'
for 1.2:
gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA BBBB?@@ABBB@=@?C>ABBCB'
There is a huge difference of the base quality between these two version. The old version seems change base quality a lot.
I failed to find such information through google. I am confused about this. Did I miss something?
Thanks
I compare the mpileup outputs between samtools 1.2 and 0.1.18. And I find
some discordance of the base quality between the old and new version.
I extract the pileup file with -B option to disable BAQ of samtools firstly. And I get a position like this:
gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA BBBB?@@ABBB@=@?C>ABBCB;
Then I perform the same command except the -B option, but for 0.1.18, I get this:
gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA <<<<>>:=;<<5=@;:>==;::'
for 1.2:
gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA BBBB?@@ABBB@=@?C>ABBCB'
There is a huge difference of the base quality between these two version. The old version seems change base quality a lot.
I failed to find such information through google. I am confused about this. Did I miss something?
Thanks
Comment