Identify similarly mapped reads in BAM/SAM?

Question

[Eliveri]

I have a BAM (or SAM) file of reads mapped to a reference genome.

I would like to identify which reads are mapped similarly/to approximately the same location, which reads might be "orthologs".

For example, when visualizing the BAM file in IGV, the two reads circled in the orange are similar.

How can I parse this information from a BAM or SAM file? Using Python or R?

Answer 1

[Pierre Lindenbaum]

 bedtools bamtobed -i in.bam | datamash -g 1,2,3  collapse 4

[Eliveri]

Thank you! This is a great and straightforward solution.