[Kwant] outputing Hamiltonian with superconductivity

If I initialize a finite system with superconductivity using 2.6.2, finalize it, and output the Hamiltonian using sys.finalized().hamiltonian_submatrix() I have no clue how the sites are ordered. It is neither grouped with particle/hole nor with sites. For example with Delta=1, mu=2, t=5, L=3, the output is [[ 3.+0.j 0.+0.j 5.-0.j 0.+0.j 1.-0.j 0.+0.j] [ 0.+0.j -3.+0.j 0.+0.j 1.-0.j 0.+0.j -5.-0.j] [ 5.+0.j 0.+0.j 3.+0.j 5.-0.j 0.+0.j 1.-0.j] [ 0.+0.j 1.+0.j 5.+0.j 3.+0.j 0.+0.j 0.+0.j] [ 1.+0.j 0.+0.j 0.+0.j 0.+0.j -3.+0.j -5.-0.j] [ 0.+0.j -5.+0.j 1.+0.j 0.+0.j -5.+0.j -3.+0.j]] Is there a way to know how is the basis ordered?