  Recent developments in (sub-)millimeter facilities have drastically changed
the amount of information obtained from extragalactic spectral scans. In this
paper, we present a feature extraction technique using principal component
analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc)
spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using
Band 3 of the Atacama Large Millimeter/submillimeter Array. We apply PCA to 16
well-detected molecular line intensity maps convolved to a common 150 pc
resolution. In addition, we include the [SIII]/[SII] line ratio and [CI]
$^3P_1$-$^3P_0$ maps in the literature, both of whose distributions show
remarkable resemblance with that of a kpc-scale biconical outflow from the
central AGN. We identify two prominent features: (1) central concentration at
the circumnuclear disk (CND) and (2) two peaks across the center that coincide
with the biconical outflow peaks. The concentrated molecular lines in the CND
are mostly high-dipole molecules (e.g., H$^{13}$CN, HC$_3$N, and HCN). Line
emissions from molecules known to be enhanced in irradiated interstellar
medium, CN, C$_2$H, and HNC, show similar concentrations and extended
components along the bicone, suggesting that molecule dissociation is a
dominant chemical effect of the cold molecular outflow of this galaxy. Although
further investigation should be made, this scenario is consistent with the
faintness or absence of the emission lines from CO isotopologues, CH$_3$OH, and
N$_2$H$^+$, in the outflow, which are easily destroyed by dissociating photons
and electrons.
