Decomposição da cadeia de caracteres
Examine as frequências relativas de códons (grupo de 3 nucleotídeos consecutivos) na lista de nucleotídeos de um gene.
Obtenha a sequência de DNA do gene humano "SCNN1A".
In[1]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_9.png)
dnasequence = GenomeData["SCNN1A", "FullSequence"];
Use StringPartition para criar a lista correspondente de códons .
In[2]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_10.png)
codons = StringPartition[dnasequence, 3];
In[3]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_11.png)
Take[codons, 10]
Out[3]=
![](assets.pt-br/string-decomposition/O_8.png)
Calcule a frequência relativa de cada códon neste gene.
In[4]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_12.png)
frequencies = N[Counts[codons]/Length[codons]];
Existem 64 possíveis códons formados a partir dos nucleotídeos A, C, G, T, e todos aparecem no gene selecionado.
In[5]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_13.png)
frequencies // Length
Out[5]=
![](assets.pt-br/string-decomposition/O_9.png)
Ache os 3 códons com as frequências mais altas.
In[6]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_14.png)
TakeLargest[frequencies, 3]
Out[6]=
![](assets.pt-br/string-decomposition/O_10.png)
Ache os 3 códons com as frequências mais baixas.
In[7]:=
![Click for copyable input](assets.pt-br/string-decomposition/In_15.png)
TakeSmallest[frequencies, 3]
Out[7]=
![](assets.pt-br/string-decomposition/O_11.png)
Visualize todas as frequências relativas em um Grid.
mostre o input completo da Wolfram Language
Out[8]=
![](assets.pt-br/string-decomposition/O_12.png)