David,just check the .DOS output (use block ProjectedDensityOfStates) - the levels that are on both sides of EF would be HOMO and LUMO.
Keep in mind that HOMO-LUMO gap is strongly underestimated in DFT (only B3LYP can give reasonable results but it is not implemented in SIESTA).
