Hello
I use some pysam function to calculate the depth of sequencing and I check the result with samtools depth. Most of the time the result is identical but sometimes there is a shift in the between the samtools result and the pysam result.
I put my code in a zip file with the example and a readme that contains the command to execute and the result files
A better place to share code would be a GitHub Gist. I'm not going to download and extract a zip :-)
My first guess is that you see this difference because samtools has a default maximum depth of 8000. I believe if you set
-d 0
you disable the maximum (if you have a reasonable recent version of samtools).please use http://gist.github.com or paste the code here, instead of this temporary link.
I made the gist but I was unable to join the bam file that provokes the error
Hello,
some more general problems that might came up if one try to compare programs which calculate read depths.
fin swimmer