The recently introduced quasi-extremal energy principle for incremental non-potential problems in rate-independent plasticity is applied to select the deformation pattern and active slip systems in single crystals. The standard crystal plasticity framework with a non-symmetric slip-system interaction matrix at finite deformation is used. The incremental work criterion for the formation of deformation bands is combined with the quasi-extremal energy principle for determining the active slip systems and slip increments in the bands. In this way, the incremental energy minimization approach has been extended to the non-potential problem of deformation banding in metal single crystals. It is shown that fulfilment of the mathematical criterion for incipient deformation banding in a homogeneous crystal in the multiple-slip case under certain conditions requires non-positive determinant of the hardening moduli matrix. Numerical examples of energetically preferable patterns of deformation bands are presented for Cu and Ni single crystals.