On 27 Jul 2016, at 15:34, Adam Witney <[email protected]> wrote:
> I am upgrading an analysis from samtools version 0.1.19-44428cd to version
> 1.3.1 (using htslib 1.3.1). I am trying to understand why the pileup changes
> so much in this example
>
> samtools mpileup -s -f genomes/NC_000962.fna alignments/sample-1.bam | head
> [mpileup] 1 samples in 1 input files
> <mpileup> Set max per-file depth to 8000
> NC_000962.3 1 T 13
> ^],^],^],^],^],^].^].^].^].^].^].^],^], ;GF/GGGGGBCDH ]]]]]]]]]]]]]
> NC_000962.3 2 T 13 ,,,,,......,, ?GFF?GDGGBGGA
> ]]]]]]]]]]]]]
> [...]
> ./samtools-1.3.1/samtools mpileup -s -f genomes/NC_000962.fna
> alignments/sample-1.bam | head
> [mpileup] 1 samples in 1 input files
> <mpileup> Set max per-file depth to 8000
> NC_000962.3 1 T 0
> NC_000962.3 2 T 0
> [...]
> Is this change a change in the default settings somewhere such that the reads
> are filtered out?
Does this difference persist throughout the genome, or just in these first few
tens of positions in NC_000962.3 and other chromosomes? This sort of
difference in mpileup read filtering often comes down to BAQ calculations... If
you rerun with samtools-0.1.19 -B and samtools-1.3.1 -B to disable BAQ
recalculation, do they come back with more similar results?
If they do and if the differences you have noticed occur at the beginning of
chromosomes but the two samtools versions produce more similar results further
on, I suspect what you are seeing is this bug fix in samtools 1.3:
• The mpileup command now applies BAQ calculations at all base positions,
regardless of which ‑l or ‑r options are used (previously with -l it was not
applied to the first few tens of bases of each chromosome, leading to different
mpileup results with -l vs. -r; #79, #125, #286, #407).
John
--
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.
------------------------------------------------------------------------------
_______________________________________________
Samtools-help mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/samtools-help
