You changed the structure and symmetries, so you can not expect the list of irreducible q and perturbations to remain the same: your input copies those for Al FCC, but the number and coordinates of the q will change.
Please check the tutorial and see how to generate the list of irred q. The simplest is to run a dummy GS calculation with kptopt 1 and ngkpt = your qpt grid. This will give you the list of irred q in the kpt variable (no need to converge or anything).
Check the new eph tutorial as well: It explains you can use qptopt 1 and iqpt with the index of the q you want.
Finally, have a look at the abipy tutorials. There are workflows which do all of this indexing for you and avoid forgetting perturbations etc.