We study a photonic implementation of a modified Fano-Anderson model - a waveguide array with two additional waveguides and by using the coupled mode theory we calculate its spectral and scattering properties. We classify eigenmodes according to vertical symmetry of the structure given by self-coupling coefficients of the additional waveguides and establish the conditions for bound states in the continuum (BIC) existence. The main predictions drawn from the theoretical model are verified by rigorous full-wave simulations of realistic structures. We use the Weierstrass factorization theorem and interpret the scattering spectra of the systems with broken symmetry in terms of the eigenmodes. The Fano resonance related with excitation of quasi-BIC is explained as arising from the interference between this mode and another leaky mode.