How to retrieve the number of reads mapped for each chromosome in the BAM file of methylation data

Question

[sata72]

I want to retrieve the number of reads mapped for each chromosome in the BAM file of methylation data. Anyone know how I can get this? Thank you in advance.

Answer 1

[bk11]

You can do this-

samtools idxstats your_data.bam | cut -f 1,3