[Bioperl-l] Merging overlapping Seq Features
Dan
dan.halligan at gmail.com
Wed Feb 17 22:04:05 UTC 2010
Hi, I'm just trying to get to grips with Bioperl as was wondering if
anyone could help me accomplish this task: I want to merge together
any overlapping sequence features, or features within xbp of each
other, that are stored in a BED file.
I have got as far as reading the features in and iterating through
them, e.g.:
my $in = Bio::FeatureIO->new(-format => 'bed', -file => "$file");
while (my $feat = $in->next_feature) {
print "start: ".$feat->start." end: ".$feat->end."\n";
}
Does Bioperl provide anything that could accomplish this merging, or
would it be something I would have to write myself?
Any help much appreciated
More information about the Bioperl-l
mailing list