Yes, aside from possibly (but hopefully not) being
used for alignment, those TRs should be removed from
the analysis via 3dTcat.
It might be worth using afni_proc.py for the analysis,
and reviewing the steps it generates, even if you want
to write your own script.
- rick