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? :)

 採用された回答

Ned Gulley
Ned Gulley 2012 年 2 月 17 日

2 投票

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 件のコメント

Nanda Lauzan
Nanda Lauzan 2012 年 2 月 18 日
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 ???

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

ヘルプ センター および File ExchangeGenomics and Next Generation Sequencing についてさらに検索

質問済み:

2012 年 2 月 17 日

編集済み:

2013 年 10 月 16 日

Community Treasure Hunt

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

Start Hunting!

Translated by