How to separate this DNA sequence?

5 vues (au cours des 30 derniers jours)
Nanda Lauzan
Nanda Lauzan le 17 Fév 2012
Modifié(e) : Image Analyst le 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? :)

Réponse acceptée

Ned Gulley
Ned Gulley le 17 Fév 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 commentaire
Nanda Lauzan
Nanda Lauzan le 18 Fév 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 ???

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Genomics and Next Generation Sequencing dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by