Info
La pregunta está cerrada. Vuélvala a abrir para editarla o responderla.
How to construct a for loop to split up a genomic sequence (~5 million nucleotides)?
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
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?
0 comentarios
Respuestas (0)
Ver también
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!