how to find out mate reads?

[hujiacun]

Hi,
When I got the BAM file of alignments in the region of interest, how to retrieve the mate's data? I want to know how those mates are distributed throught the whole genome.

[student-t]

Can you just visualize them on a viewer such as IGV?

[hujiacun]

Can you just visualize them on a viewer such as IGV?

student-t, thanks for your reply.
I use "samtools view originalfile.bam Chr4:100k-200k" to get a regional BAM file. In IGV, I can visualize the reads among my interested region. Because of Structure Viaration, the mates of those reads may not locate in the same region, and are not included in the BAM file.
How can I get a file including data of those mates?

[maubp]

Do you want to do this by hand/eye in a viewer? If so look at the full BAM file and the viewer should be able to show you the partner read (try right clicking etc). This is possible with Tablet via the tool tip, and you can right click to jump to the partner etc:

https://ics.hutton.ac.uk/tablet/

Or, you can do this programatically by parsing the SAM/BAM file (e.g. in Python using pysam). There are specific columns used to report the mapping position of the partner read.