Dear Samtools team,
The samtools depth command sets a max depth of 8000 to the aligned reads.
Previously the max threshold of 8000 was replaced by adding the
bam_mplp_set_maxcnt(mplp, 1000000)
line in the "bam2depth.c" file.
In the 1.3.1 release of Samtools this has been corrected, and "bam2depth.c"
already contains the following lines:
if (0 < max_depth)
bam_mplp_set_maxcnt(mplp,max_depth);
The max_depth is set with the -d option of the mpileup command.
Is there also a way to specify the max_depth for the depth command?
Thanks,
Lilit
--
Lilit Nersisyan
Junior Researcher, MSc
Institute of Molecular Biology NAS RA
7, Hasratyan str., 0014, Yerevan, Armenia
Tel.: +37494 601703, +37410 282622
------------------------------------------------------------------------------
What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
patterns at an interface-level. Reveals which users, apps, and protocols are
consuming the most bandwidth. Provides multi-vendor support for NetFlow,
J-Flow, sFlow and other flows. Make informed decisions using capacity
planning reports. https://ad.doubleclick.net/ddm/clk/305295220;132659582;e
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help