Saturday, July 18, 2020

mutations ~ dewpoint in spike region


There are 3 mutations in the Spike region.





msa_position p q ref_position ref_seq type gene start end protein_id feature_position feature_position_mod3 triplet_start triplet_end ref_triplet mut_triplet ref_position_inProtein ref_AA ref_AAPosition mut_AA ref_freq mut_freq mut_seq CHR
3654 36985 0.00453724672342917 0.0019392834243773 23730 C CDS S 21563 25384 YP_009724390.1 2169 0 23728 23730 ACC ACT 7910 T T723 T 40314 1381 T 1
3655 37185 0 0 23928 C CDS S 21563 25384 YP_009724390.1 2367 0 23926 23928 TAC TAT 7976 Y Y789 Y 41149 530 T 1
3653 36626 1.98774342659889E-68 4.81792614569765E-68 23402 A CDS S 21563 25384 YP_009724390.1 1841 2 23401 23403 GAT GGT 7801 D D614 G 11231 30413 G 1

ef_tripletmut_tripletref_position_inProteinref_AAref_AAPositionmut_AAref_freqmut_freq
GATGGT7801DD614G1123130413G
ACCACT7910TT723T403141381T
TACTAT7976YY789Y41149530

Are synonymous codon adaptation in play?
In bat, ACT > ACC, in human, ACC > ACT, so T723T from ACC -> ACT change from human major codon to human minor codon.

In human, TAT > TAC, so Y789Y change from human major codon to minor codon.
So, the answer seems to be NO.

Bat virus codon usage versus human codon bias usage

West Caucasian bat virus [gbvrl]: 4 CDS's (1478 codons)
http://www.kazusa.or.jp/codon/cgi-bin/showcodon.cgi?species=249584

West Caucasian bat virus [gbvrl]: 4 CDS's (1478 codons)

fields: [triplet] [frequency: per thousand] ([number])

UUU 18.9(    28)  UCU 29.8(    44)  UAU 22.3(    33)  UGU 12.9(    19)
UUC 21.0(    31)  UCC 18.9(    28)  UAC 18.3(    27)  UGC  8.8(    13)
UUA 10.8(    16)  UCA 18.9(    28)  UAA  0.7(     1)  UGA  1.4(     2)
UUG 20.3(    30)  UCG  4.1(     6)  UAG  0.7(     1)  UGG 12.9(    19)

CUU 14.2(    21)  CCU 18.3(    27)  CAU 12.9(    19)  CGU  3.4(     5)
CUC 12.9(    19)  CCC 13.5(    20)  CAC 16.2(    24)  CGC  1.4(     2)
CUA  8.8(    13)  CCA 14.2(    21)  CAA 16.9(    25)  CGA  2.7(     4)
CUG 16.2(    24)  CCG  2.7(     4)  CAG 10.1(    15)  CGG  1.4(     2)

AUU 16.2(    24)  ACU 15.6(    23)  AAU 26.4(    39)  AGU 16.2(    24)
AUC 27.1(    40)  ACC 12.2(    18)  AAC 17.6(    26)  AGC  6.8(    10)
AUA 18.3(    27)  ACA 22.3(    33)  AAA 23.7(    35)  AGA 29.8(    44)
AUG 23.7(    35)  ACG  3.4(     5)  AAG 35.2(    52)  AGG 17.6(    26)

GUU 17.6(    26)  GCU 15.6(    23)  GAU 37.9(    56)  GGU 10.1(    15)
GUC 18.9(    28)  GCC 13.5(    20)  GAC 23.0(    34)  GGC  6.1(     9)
GUA  9.5(    14)  GCA 14.2(    21)  GAA 38.6(    57)  GGA 21.7(    32)
GUG 21.0(    31)  GCG  2.0(     3)  GAG 32.5(    48)  GGG 19.6(    29)

Coding GC 43.96% 1st letter GC 46.75% 2nd letter GC 39.17% 3rd letter GC 45.94%
homo sapiens
http://www.kazusa.or.jp/codon/cgi-bin/showcodon.cgi?species=9606
Homo sapiens [gbpri]: 93487 CDS's (40662582 codons)

fields: [triplet] [frequency: per thousand] ([number])

UUU 17.6(714298)  UCU 15.2(618711)  UAU 12.2(495699)  UGU 10.6(430311)
UUC 20.3(824692)  UCC 17.7(718892)  UAC 15.3(622407)  UGC 12.6(513028)
UUA  7.7(311881)  UCA 12.2(496448)  UAA  1.0( 40285)  UGA  1.6( 63237)
UUG 12.9(525688)  UCG  4.4(179419)  UAG  0.8( 32109)  UGG 13.2(535595)

CUU 13.2(536515)  CCU 17.5(713233)  CAU 10.9(441711)  CGU  4.5(184609)
CUC 19.6(796638)  CCC 19.8(804620)  CAC 15.1(613713)  CGC 10.4(423516)
CUA  7.2(290751)  CCA 16.9(688038)  CAA 12.3(501911)  CGA  6.2(250760)
CUG 39.6(1611801)  CCG  6.9(281570)  CAG 34.2(1391973)  CGG 11.4(464485)

AUU 16.0(650473)  ACU 13.1(533609)  AAU 17.0(689701)  AGU 12.1(493429)
AUC 20.8(846466)  ACC 18.9(768147)  AAC 19.1(776603)  AGC 19.5(791383)
AUA  7.5(304565)  ACA 15.1(614523)  AAA 24.4(993621)  AGA 12.2(494682)
AUG 22.0(896005)  ACG  6.1(246105)  AAG 31.9(1295568)  AGG 12.0(486463)

GUU 11.0(448607)  GCU 18.4(750096)  GAU 21.8(885429)  GGU 10.8(437126)
GUC 14.5(588138)  GCC 27.7(1127679)  GAC 25.1(1020595)  GGC 22.2(903565)
GUA  7.1(287712)  GCA 15.8(643471)  GAA 29.0(1177632)  GGA 16.5(669873)
GUG 28.1(1143534)  GCG  7.4(299495)  GAG 39.6(1609975)  GGG 16.5(669768)

Coding GC 52.27% 1st letter GC 55.72% 2nd letter GC 42.54% 3rd letter GC 58.55%

No comments:

Post a Comment