Figure 3.
Increased m6A calling on long reads (>20 kb) via fibertools. (A) Percentage increase in fibertools m6A calls over the GMM model as a function of the minimum read length of the underlying sequencing data. The histogram below shows how many reads were used to calculate each percentage increase. (B) Comparison of fibertools and the subread model for m6A calling over CAGE-positive TSS in K562 cells across the genome, separated by read length. Reads are matched between fibertools and the subread model, and the number of Fiber-seq reads used in the calculation of each size range (n) is indicated.
