UPDATE (Fri 16 Apr 2021 06:58:42 PDT): Use goseq instead of a customized script.
goseq
pathoAD BP
pathoAD MF
pathoAD CC
apoe.e4 BP
apoe.e4 MF
apoe.e4 CC
msex BP
msex MF
msex CC
cog BP
cog MF
cog CC