Hello,
I was calling variants from RNAseq data using samtools mpileup|bcftools
call. I am aware that the BAQ (base alignment quality) strategy is invoked
by default in mpileup. BAQ is in itself derived from the original base
quality score. Now if the BAQ is lower, no variants will be called in that
base. What is the general range of BAQ? How much is high and how much is
low? Is the BAQ score reflected in the VCF file in any form?
I was also wondering for a way visualize an mpileup file? What can be the
best way to do it?
--
Regards,
Prasun Dutta, M.Sc. (Bioinformatics)
PhD Student, Developmental Biology Division (Hume Group)
The Roslin Institute, The University of Edinburgh, Scotland, UK
M: +447438743406
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help
