how to split dna sequence into three letters each

11 次查看（过去 30 天）
Aarsha mv
Aarsha mv 2017-5-12
回答： Luuk van Oosten 2017-5-12
how to split a dna sequence into three letters each
  3 个评论
显示 1更早的评论
Stephen23
Stephen23 2017-5-12
编辑：Stephen23 2017-5-12
Pick from:
reshape
reshape and num2cell
mat2cell
regexp
...
John D'Errico
John D'Errico 2017-5-12
Get used to working in MATLAB where you think about variables in terms of their shape, their size, as complete arrays. Think of things in terms of how you can transform those arrays into what will be useful to you. Then do some quick reading in the help of the tools in matfun.
help elmat
There are many useful tools in elmat. Some of them, like reshape, might be useful to you.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Luuk van Oosten
Luuk van Oosten 2017-5-12
Dear Aarsha mv,
I assume you want to get the codons from your piece of DNA, but correct me if I am wrong.
Let us take the example of a piece of DNA coding for the protein insulin (I took a part of this piece of DNA).
Define a string of nucleotides as being your DNA sequence:
your_DNA = 'CTCGAGGGGCCTAGACATTGCCCTCCAGAGAGAGCACCCAACACCCTCCAGGCTTGACCGGCCAGGGTG';
To get the codons you can use 'reshape', as was suggested by Stephen Cobeldick and John D'Errico.
Use reshape as follows:
codons = reshape(your_DNA(:),3,length(your_DNA)/3)'
In a 'real life' scenario you would probably need a workaround if length(your_DNA) is not a number that can be divided by 3.
Besides having a look at the tips that Cobeldick and d'Errico already gave you, I suggest you have a look at the Nucleotide Sequence Analysis overview page. There might be some functions that will aid you in your DNA analysis.

类别

在 Help Center 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by