hi

[indira mekala]

hi,


how do i take out the sequences of length greater than 100kb from a whole lot of sequences, please sent that programm.

[Jason Stajich]

Indira - 

What have you tried in terms of writing this yourself? Showing some of your code would be good. 

The bp_seq_length.pl does one aspect of this,  hopefully what you've learned about when it comes to Bio::SeqIO should let you figure out how to put it all together.

This one calculates total sequence length, can you modify it to just see if the length is > than your size.

https://github.com/bioperl/bioperl-live/blob/master/scripts/utilities/bp_seq_length.pl

Also see these tutorials

http://www.molecularevolution.org/resources/activities/bioperl

On Nov 22, 2011, at 12:45 AM, indira mekala wrote:

hi,


how do i take out the sequences of length greater than 100kb from a whole lot of sequences, please sent that programm.

-- 
Jason E Stajich, PhD
Assistant Professor, Plant Pathology & Microbiology

Institute for Integrative Genome Biology
University of California, Riverside
951.827.2363

http://lab.stajich.org http://fungalgenomes.org http://fungidb.org 

twitter: @stajichlab @hyphaltip @fungalgenomes @fungidb
http://plantpathology.ucr.edu http://genomics.ucr.edu