Fig. 2
From: LEMON: a method to construct the local strains at horizontal gene transfer sites in gut metagenomics
a Some paired-end reads cross HGT breakpoints s1 and s2. b The two sides of one junction pair are mapped to two different references. Bx and By are mapped positions of junction pair B on reference X and Y respectively. c Furthermore, all junction pairs mapped to the same two references are transformed to points in a coordinate system with their mapped positions as coordinates, e.g. junction pairs A and B are transformed to points A′(Ax,Ay) and B′(Bx,By). d We apply DBSCAN using Euclidean distance to cluster the breakpoint pair candidates, such as S1 and S2. The maximum diameter of the cluster circle is the insert size of pair-end reads