Hi, Shinyoung Jung-

Typically, one wouldn't directly align a statistical dataset to a template, because they statistics map doesn't have clear anatomical features (the kinds of data contained there are effect estimate and statistics blobs that might extend over structural boundaries).

For FMRI, during processing prior to statistical mapping, we might align the subject EPI to their anatomical dataset (creating the EPI-anatomical transformation), and then that anatomical dataset to the MNI template (creating the anatomical-template transformation), and combine/concatenate those transformations procedures into a single one (EPI-to-template transformation) to apply to the EPI so that it overlays well on the template. Then, when GLM modeling is performed, the results are already well aligned to the template, having used structural information in the EPI and anatomical datasets.

This is discussed, along with other alignment procedures that are part of standard processing and *also* concatenated into that single EPI-to-template transform, here:
https://www.youtube.com/watch?v=PaZinetFKGY&list=PL_CD549H9kgqJ1GDXAs1BWkgEimAHZeNX
And many other education videos are also on that AFNI Academy channel:
https://www.youtube.com/c/afnibootcamp

And actually, it is generally much more convenient to combine the alignment, statistical modeling, and other FMRI processing steps into a single pipeline that can be adjusted with options to do exactly what you want. In AFNI, this would be done with afni_proc.py, which is described here:
https://www.youtube.com/watch?v=XNPNAgpmjxI&list=PL_CD549H9kgpGgC8xWFexw-8R95HZokAr
... and there are several starter examples in the afni_proc.py help, too:
https://afni.nimh.nih.gov/pub/dist/doc/htmldoc/programs/alpha/afni_proc.py_sphx.html#ahelp-afni-proc-py

--pt