Info

This question is closed. Reopen it to edit or answer.

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

1 view (last 30 days)
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?

Answers (0)

Products

Community Treasure Hunt

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

Start Hunting!