How to separate this DNA sequence?

5 views (last 30 days)
Nanda Lauzan
Nanda Lauzan on 17 Feb 2012
Edited: Image Analyst on 16 Oct 2013
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? :)

Accepted Answer

Ned Gulley
Ned Gulley on 17 Feb 2012
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
  1 Comment
Nanda Lauzan
Nanda Lauzan on 18 Feb 2012
thanks Gulley for the answer :)
but if codons =
AAA
TTT
ATG
TGA
CAG
TAG
and change into this :
codons1 = AAA
codons2 = TTT
codons3 = ATG
codons4 = TGA
codons5 = CAG
codons6 = TAG ???

Sign in to comment.

More Answers (0)

Categories

Find more on Genomics and Next Generation Sequencing in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!