how to split dna sequence into three letters each
显示 更早的评论
how to split a dna sequence into three letters each
3 个评论
KSSV
2017-5-12
how is dna sequence? what format? One example is needed to help.
Pick from:
reshape
reshape and num2cell
mat2cell
regexp
...
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
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.
类别
在 帮助中心 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息
2017-5-12
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
