MATLAB: How to separate this DNA sequence

bioinformaticsdna

Hey everyone, I've the problem for separating this DNA sequence. for example :
 sequence = 'AAATTTATGTGACAGTAG';
i've tried like this :
 [one, two] = strtok(sequence)
but i've a result like this
 one =
  AAATTTATGTGACAGTAG
 two =
   Empty string: 1-by-0
i just want the result
 one = AAA
 two = TTT
 three = ATG
 four = TGA
 five = CAG
 six = TAG
and so on. can you all fix this? 🙂

Best Answer

Here's a simple way to chunk your DNA up into codons. Then you can just pull the rows off to get each triplet.
 >> sequence = 'AAATTTATGTGACAGTAG';
 >> codons = sequence;
 >> codons = reshape(codons(:),3,length(codons)/3)'
 codons =
 AAA
 TTT
 ATG
 TGA
 CAG
 TAG
Related Question