How to separate this DNA sequence?

5 visualizzazioni (ultimi 30 giorni)
Nanda Lauzan
Nanda Lauzan il 17 Feb 2012
Modificato: Image Analyst il 16 Ott 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? :)

Risposta accettata

Ned Gulley
Ned Gulley il 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 Commento
Nanda Lauzan
Nanda Lauzan il 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 ???

Accedi per commentare.

Più risposte (0)

Categorie

Scopri di più su Genomics and Next Generation Sequencing in Help Center e File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by