How to separate this DNA sequence?
5 views (last 30 days)
Show older comments
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? :)
0 Comments
Accepted Answer
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
More Answers (0)
See Also
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!