I'm trying to get the number of reads that start at a given position. I'm using bbmap pileup.sh. This is the script:
./pileup.sh in=sorted.bam basecov=coverage.txt startcov=t
How can I get coordinate information. For example, chrM 1033: 10 (so 10 reads start at that position).
Thank you