Say I have a read with two supplementary alignments. How do I infer the gap length of the two aligned segments relative to the read?
In the example above, the gap length is the size of the insertion (orange thin line). I'm confused as heck as how to calculate the gap length from reads in a BAM file. Would the CIGAR string be informative, or do I have to go back to the original FASTQ?