I wanted to build an awesome place for people to discuss module specific issues, but I don't have any more time for this, and there are much better places to discuss Perl-related issues. I'd recommend asking your question on Stack Overflow or on Perl Monks.
If you are looking for a Perl tutorial or Perl-related news, I hope these links will serve you well.
Posted on 2007-04-10 08:54:03-07 by gopu
extract nonoverlapping subsequences from a whole genome
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!
Direct Responses: 5081 | Write a response