Dear all,

I have multiple sequence aligned data and I want to calculate position wise KaKs. Previously I was using SELECTON/2.4 (http://selecton.tau.ac.il/versions.html), but the current version is out dated as it was released 12 years in 2011.

SELECTON 2.4 can only be compiled with GCC 4.9.2 but not with the newer versions of GCC.

I am looking for an alternative solution, your feedback will be highly appreciated.

Thank you.

The type of data (for example) is:

>ENSVURG00010003966

-------------------------------------------------------MASPS

RRLQTKPVISCLKSVLLTYTFIFWFTGVILLTVGIWGKVSLEVYFSLLNEKATNVPFVLI

GP-GTIIILLGTFGCFATSRASAWMLKLYAMFLTLIFFVELVAAIVGFVFRHEIRNSFEM

NYEASLKQYNDTA-IQRS--EAVDTIQKTLHCCGVQSAVDWETTAYYAE--KGFPRSCCR

FHSC--SWADRKDLKKQI-----QEGCFVKVTTVIQSEMGVVAGISFGIACFQ-------

----------LFGIFLACCLSRSITNNQYEIV----------------------------

------------------------------------------------------------

-------------------------------

>ENSPCIG00000003889

-------------------------------------------------------MASPS

RRLQTKPVITCLKSVLLTYTFIFWFTGVILLTVGIWGKVSLEVYFSLLNEKATNVPFVLI

GT-GTVIILLGTFGCFATCRASAWMLKLYAMFLTLIFLVELVAAIVGFVFRHEIKNSFKM

NYEASLKQYNDTG-IQRS--EAVDTIQKTLHCCGVQSATDWETTAYYAE--KGFPRSCCR

FHSC--TWADRKDLKKQS-----QEGCFVKVTTVIQSEMGVVAGISFGVACFQ-------

----------LIGIFLACCLSRSITNNQYEIV----------------------------

------------------------------------------------------------

-------------------------------

>ENSSHAG00000022159

-------------------------------------------------------MASPS

RRLQTKPVITCLKSVLLTYTFIFWFTGVILLTVGIWGKVSLEVYFSLLNEKATNIPFVLI

AT-GTVIILLGTFGCFATCRASAWMLKLYAMFLTLIFLVELVAAVVGFVFRHEIKDSFKM

NYEASLKQYNDTG-VQRS--EAVDTIQKTLHCCGVLSATDWESTAYYVE--KGFPRSCCR

YHSC--SWADRKDLKKQN-----QEGCFLRVTTIIQSEMGVVAGISFGVACFQ-------

----------LIGIFLACCLSRSITNNQYEIV----------------------------

------------------------------------------------------------

-------------------------------

>ENSG00000000003

-------------------------------------------------------MASPS

RRLQTKPVITCFKSVLLIYTFIFWITGVILLAVGIWGKVSLENYFSLLNEKATNVPFVLI

AT-GTVIILLGTFGCFATCRASAWMLKLYAMFLTLVFLVELVAAIVGFVFRHEIKNSFKN

NYEKALKQYNSTG-DYRS--HAVDKIQNTLHCCGVTDYRDWTDTNYYSE--KGFPKSCCK

LEDC--TPQRDADKVNN-------EGCFIKVMTIIESEMGVVAGISFGVACFQ-------

----------LIGIFLAYCLSRAITNNQYEIV----------------------------

------------------------------------------------------------

-------------------------------

and the expected outcomes, that I am looking for

Selecton Bayesian Ka/Ks Results

Displayed on sequence ENSVURG00010003966 including gaps

==========================================================================================================================================

POS AMINO Ka/Ks [Confidence Interval](* if lower bound > 1) POSTERIOR PROBABILITIES

w = 0.0058 0.028 0.062 0.11 0.17 0.25 0.36 0.59 1

==========================================================================================================================================

1 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

2 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

3 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

4 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

5 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

6 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

7 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

8 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

9 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

10 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

11 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

12 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

13 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

14 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

15 - 0.22 [0.0058,1] 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.028

16 - 0.39 [0.062,1] 0.00044 0.0086 0.033 0.075 0.13 0.18 0.23 0.28 0.067

17 - 0.4 [0.062,1] 0.00031 0.0065 0.027 0.066 0.12 0.18 0.24 0.29 0.064

18 - 0.44 [0.11,1] 0.0002 0.0042 0.018 0.048 0.095 0.16 0.24 0.34 0.087

19 - 0.45 [0.11,1] 1.2e-05 0.0011 0.0096 0.036 0.089 0.17 0.26 0.35 0.085

20 - 0.57 [0.17,1] 1.2e-07 5.4e-05 0.0011 0.0072 0.029 0.086 0.21 0.48 0.19

21 - 0.56 [0.17,1] 4.2e-06 0.0004 0.0035 0.014 0.041 0.095 0.2 0.45 0.19

22 - 0.5 [0.11,1] 6.3e-06 0.00062 0.0056 0.023 0.063 0.14 0.25 0.41 0.12

23 - 0.5 [0.11,1] 6.5e-06 0.00063 0.0057 0.023 0.063 0.13 0.24 0.41 0.13

24 - 0.29 [0.028,1] 0.015 0.06 0.1 0.13 0.15 0.16 0.17 0.17 0.039

25 - 0.55 [0.11,1] 5.8e-05 0.0014 0.0067 0.02 0.049 0.1 0.2 0.43 0.19

26 - 0.4 [0.062,1] 0.00028 0.0058 0.025 0.062 0.12 0.18 0.25 0.3 0.067

27 - 0.52 [0.17,1] 4.5e-06 0.00045 0.0042 0.018 0.052 0.12 0.23 0.43 0.14

28 - 0.43 [0.062,1] 0.00027 0.0056 0.023 0.057 0.11 0.17 0.24 0.32 0.085

29 - 0.48 [0.11,1] 6.9e-06 0.00068 0.0061 0.025 0.068 0.14 0.25 0.4 0.11

30 - 0.34 [0.062,1] 0.00055 0.011 0.043 0.095 0.16 0.21 0.24 0.21 0.035

31 - 0.54 [0.17,1] 2.9e-06 0.00031 0.0032 0.015 0.046 0.11 0.23 0.44 0.16

32 - 0.21 [0.028,0.59] 0.02 0.082 0.14 0.18 0.19 0.17 0.13 0.078 0.011

33 - 0.44 [0.062,1] 0.00026 0.0054 0.023 0.055 0.1 0.16 0.23 0.32 0.1

34 - 0.47 [0.11,1] 8.8e-06 0.00085 0.0074 0.029 0.076 0.15 0.26 0.38 0.096

35 - 0.52 [0.11,1] 6.7e-06 0.00064 0.0055 0.022 0.057 0.12 0.23 0.42 0.14