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?