Towards realistic codon models: among site variability and dependency of synonymous and non-synonymous rates