RNA modifications are key regulators for RNA processes. tRNA-derived RNAs are small RNAs with size between 15 and 50 bases long that are processed from mature or precursor tRNAs. Despite their more recent discovery, tRNA-derived RNAs have been found to play regulatory roles in many cellular processes including gene silencing, protein synthesis, stress response, and transgenerational inheritance. Furthermore, tRNA-derived RNAs are highly abundant in bodily fluids, posing as potential biomarkers. A unique feature of tRNA-derived RNAs is that they are rich in RNA modifications. Many of the RNA modifications on tRNA-derived RNAs disrupt Watson-Crick base pairing and will thus stall reverse transcriptase, such as N(1)-methyladenosine (m(1)A), N(1)-methylguanosine (m(1)G) and N(2), N(2)-dimethylguanosine (m(2)(2)G). These RNA modifications add another layer of regulation onto tRNA-derived RNAs' functions and are of interests for future research. However, these RNA modifications could also lead to lower detection of modification-containing RNAs in genome-wide small RNA sequencing analysis due to reverse transcriptase stall. To circumvent this bias, TGIRT (Thermostable Group II Intron Reverse Transcriptase) has been used to readthrough RNA modifications inserting mismatches. These mismatch signatures can then be used to precisely map the modification sites at base resolution. Here we describe the step-by-step experimental protocol to start with purified RNAs from cells or tissues and use TGIRT to make small RNA sequencing library for Illumina sequencing to profile the abundance of tRNA-derived RNAs and the associated RNA modifications.