Info

Questa domanda è chiusa. Riaprila per modificarla o per rispondere.

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

1 visualizzazione (ultimi 30 giorni)
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?

Risposte (0)

Prodotti

Community Treasure Hunt

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

Start Hunting!

Translated by