I would like to compute a Lorenz curve of per-base coverage depth for a BAM file on a certain BED file. I can find a way to do it by myself, but it would be great if you could point me to some existing software.
Thank you in advance!
I would like to compute a Lorenz curve of per-base coverage depth for a BAM file on a certain BED file. I can find a way to do it by myself, but it would be great if you could point me to some existing software.
Thank you in advance!
Hi Igor,
I wrote a python tool to do so:
https://github.com/yhoogstrate/bam-lorenz-coverage
all best,
Youri
htSeqTools (with giniCoverage function) seems the only ready-to-use tool I have found for this task. However, it does not allow to restrict the analysis on a specific region (for example, by using a BED file).
https://blog.addgene.org/ngs-quality-control-for-pooled-libraries you can check the comments there is a comment by will including a link which is plsa.py
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.