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!
•
written
by
igor.fobia •
3 answers
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).
•
written
by
igor.fobia •
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
Log in to answer this question.
-
Gene performance/coverage calculation methodology in .bam files using samtools
written by slzr_ •
-
Using Bedtools coverage for a single alignment file?
written by caro-ca •
-
How to obtain zero-based coordinates read depth using bedtools coverage for a specific region?
written by vladimiralan •
-
samtools depth -r
written by lz1903 •
-
Merging and then unmerging small contigs
written by Philipp Bayer
-
Gini Coefficient and Lorenz curve to measure uniformity of coverage
written by am •
-
GC Bias in Sequence Capture - Analysis & Quantification
written by jnowacki
-
Get genes and coverage/depth per gene from sam/bam
written by prishly •
-
Median coverage from BAM file
written by firestar
-
Depth of Coverage - Calculated Per Probe / Bed Region
written by jnowacki
