On Wed, Jul 12, 2017 at 04:41:17PM +0200, Kristian Ullrich wrote: > whatever I am doing with both 'samtools mpileup' and 'bcftools > mpileup', in both cases not all sites are reported, even if I > specify the '-a -a' options e.g.:
My standard first try at explaining anything odd with mpileup is to test it with (default) and without (-B) the BAQ calculation. Also consider -A too. Does adding -B change things for you? (If it does, explaining why and deciding if it is an improvement is a totally different and more complex issue!) James -- James Bonfield (j...@sanger.ac.uk) | Hora aderat briligi. Nunc et Slythia Tova | Plurima gyrabant gymbolitare vabo; A Staden Package developer: | Et Borogovorum mimzebant undique formae, https://sf.net/projects/staden/ | Momiferique omnes exgrabure Rathi. -- The Wellcome Trust Sanger Institute is operated by Genome Research Limited, a charity registered in England with number 1021457 and a company registered in England with number 2742969, whose registered office is 215 Euston Road, London, NW1 2BE. ------------------------------------------------------------------------------ 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