To identify the functional 3′UTR SNPs, i.e. those whose ratio between reference and alternative allele was significantly different between the vector-expressed RNA and the plasmid DNA extracted from the same transfected cells, a generalized linear model was used. Similar to ASE analysis, a negative binomial distribution was used to account for the over-dispersion effects on the UMI count data.