Dear Cheila,
In this case the -G is irrelevant. NielsenYang.bf uses codon models of evolution which include their own way to correct for rate variation across sites (this is actually the basis for detecting codon-level positive selection). The -G part only applies if you want to use nucleotide models to analyze your data. 012345 would be the appropriate input (using, for example, the MG94 custom model).
If you'd like, you can also perform model selection and analysis for positive selection using Multimedia File Viewing and Clickable Links are available for Registered Members only!! You need to
, which implements three different methods for finding positively and negatively selected codons. The original Nielsen-Yang can falsely identify hypervariable sites (e.g. relaxed functional constraint, along with elevated mutation rates) as positively selected. We recommend the use of models which allow both synonymous and non-synonymous rates to vary from codon to codon (NY assumes that synonymous rates are constant = 1, and non-synonymous rates vary from codon to codon).
If you are interested, or just want to learn more about how modeling assumptions can affect selection analyses, I would suggest that you take a look at Multimedia File Viewing and Clickable Links are available for Registered Members only!! You need to
and Multimedia File Viewing and Clickable Links are available for Registered Members only!! You need to
Cheers,
Sergei