The transcriptional programs activated by p53 in B cells in vivo following exposure to ionizing radiation were studied through the integrated analysis of various types of next-generation sequencing data: genome-wide profiling of p53 binding sites, mapping of histone marks and open chromatin regions and quantification of gene expression. Moreover, the binding of p53 was associated to a series of specific motifs on the DNA, which were directly inferred from the data. Here, we describe in detail the computational analysis of the datasets associated with our study (Tonelli et al., Oncotarget 6 (2015), 24611-26), deposited in the GEO archive (accession code GSE71180), and we provide the R scripts needed to generated the figures of the paper.