Info

Cette question est clôturée. Rouvrir pour modifier ou répondre.

How to construct a for loop to split up a genomic sequence (~5 million nucleotides)?

1 vue (au cours des 30 derniers jours)
artfan100
artfan100 le 27 Juil 2013
Clôturé : MATLAB Answer Bot le 20 Août 2021
Hi! I'm constructing a script to basically find all promoter sequences in a bacterial genome. Problem is, I can't manage to split up the 5 million base nucleotide genome sequence into manageable segments so I can run seqshoworf command (it only works with max 1 million nucleotides).
here's the code i have so far:
putida = getgenbank('NC_009512.1','SequenceOnly',true);
putida1 = putida(1:100000);
putida2 = putida(100001:200000);
putida3 = putida(200001:300000); %etc
for i=1:1:60
putida'num2str(i)' = putida(100000*i-99999:100000*i);
end
Could someone correct the for loop so that I can work separately with each putida(n) segment?

Réponses (0)

Produits

Community Treasure Hunt

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

Start Hunting!

Translated by