Eastern Indonesia is tectonically complex, formed by different plates and microplates interactions from different origins. This complexity gives geoscientists a challenge to solve the ’jigsaw’ of the complex interactions. The understanding of tectonic processes can lead to a breakthrough in both resource exploration and disaster risk reduction. We utilize teleseismic P wave coda for random coda from scattering and deterministic coda originated from the crust-mantle boundary (Moho) to derive the crustal properties, including thickness, Vp/Vs, and qualitative scattering characteristics. For the scattering properties, we apply Iterative Cross-Correlation and Stacking (ICCS) to align the waveform. At the same time, for the crust characteristic, we employ the Receiver Functions (RF) method alongside H-k stacking. The crustal thickness recovered from the RF and H-k stacking has a good correlation with the crustal origin, where the thickness in older and stable crust originated from Sundaland and Gondwana is thicker than a younger plate of the crust arc and subduction origin. The Vp/Vs is high in a region that is interpreted to be dominated by mafic lower crust originated from oceanic-oceanic subduction during Eocene, anisotropy, or by a magmatic anomaly. The P coda also correlated well with the subsurface magmatic anomaly by providing a unique pattern.