|
Hi,
I am one of the newbee venturingout bioperl for my research purposes. I have a whole genome sequence of a pathogen. I am trying to break them into non-overlapping 1000bps subsequences. For example if my whole genome sequence is 400000 bps length, then I should be beak them into 4000 subsequences of each 1000 bps and they should be non-overlapping but at the same time continous. To be precise, my first substring would be from 1 to 1000 bps, second substing would be from 1001 to 2000 etcc.. Could anyone help me.
I tried with the following code but it gives me only the first substring and rest are not! I would appreciate very much if someone could help me!
|