How to count the number of reads in each chromosome in a bam file? The bam file is already sorted by the chromosome names. If the bam file is indexed, you may quickly get these info from the index: samtools idxstats in.bam | awk '{print $1" "$3}' If the bam file is not indexed, you […]
The post How to count the number of reads in each chromosome in a bam file? appeared first on SysTutorials.