saf <- scrapR::parse_saf_pf("polyadb32.hg38.saf.gz")

mat <- scrapR::scraps_to_matrix("R2_counts.tsv.gz", pf = saf)

so <- scrapR::scraps_to_seurat("R2_counts.tsv.gz", pf = saf, so)

res <- PA_DEXSeq(mat, cell_ids1, cell_ids2)

motifs_prep (rna, id1, id2) # by default creates a file hs_motifs_05_oRNAament_rbp.csv.gz
mdb <- read_csv("hs_motifs_05_oRNAament_rbp.csv.gz")
DE_PA_motifs(mdb, res)