Thiago Venancio | 27 Feb 21:47 2015

MSG: Start/stop out of range | Bio::Search::Tiling::MapTiling

Hi all,

By running the following code in the BLAST report attached I get the following error:

------------- EXCEPTION: Bio::Root::Exception -------------
MSG: Start/stop out of range [824, 1141]
STACK: Error::throw
STACK: Bio::Root::Root::throw /usr/local/share/perl/5.20.1/Bio/Root/
STACK: Bio::Search::HSP::HSPI::matches_MT /usr/local/share/perl/5.20.1/Bio/Search/Tiling/
STACK: Bio::Search::Tiling::MapTiling::_calc_stats /usr/local/share/perl/5.20.1/Bio/Search/Tiling/
STACK: Bio::Search::Tiling::MapTiling::length /usr/local/share/perl/5.20.1/Bio/Search/Tiling/
STACK: Bio::Search::Tiling::MapTiling::num_aligned /usr/local/share/perl/5.20.1/Bio/Search/Tiling/
STACK: Bio::Search::Tiling::MapTiling::frac_aligned /usr/local/share/perl/5.20.1/Bio/Search/Tiling/

I used the following code:

use strict;
use warnings;
use Bio::SearchIO;
use Bio::Search::Tiling::MapTiling;

my $searchio = Bio::SearchIO->new( -format => 'blast',
  -file   => shift );

my $gene;

while ( my $result = $searchio->next_result() ) {

    $result->query_name =~ /(\S+_\S+)_i/;
    my $g = $1;

    while( my $hit = $result->next_hit ) {
my $tiling = Bio::Search::Tiling::MapTiling->new($hit);
my $hsp    = $hit->next_hsp;

my $st    = $hsp->query->strand == "1" ? "p" : "m";
my $frame = $hsp->query->frame;
my $context = $st.$frame;

print $tiling->frac_aligned('query', 'exact', $context)."\n";

Any clues on how to deal with this error ?

Thanks in advance.

Bioperl-l mailing list
Bioperl-l <at>
Thiago Venancio | 27 Feb 19:21 2015

frac_identical | Bio::Search::Tiling::MapTiling

Hi all,

Bio::Search::Tiling::MapTiling frac_identical() is calculating it over the total query length instead of the aligned length in a a BLASTN. Below is the code I am using:

while( my $hit = $result->next_hit ) {
     my $tiling = Bio::Search::Tiling::MapTiling->new($hit);
     my $id = $tiling->frac_identical(-type => 'query', -action => 'exact');

Does anyone know how to change this behavior ?


! iv>
Bioperl-l mailing list
Bioperl-l <at>
Raoul J P Bonnal | 23 Feb 09:51 2015

OBF GSoC Registration

Please all possible mentors,

register yourself

and try to connect to the organization Open Bioinformatics Foundation
id: obf

in case you are aware that some mentor can not be reached by this 
message feel free to forward it.

Ra & Fra

On 2/20/15 8:51 PM, Raoul Bonnal wrote:
> Dear All,
> I have created a special section
> I think that would be very useful to add all the references where GSoC has been used for doing science.
> please go ahead and add your project/paper/or whatever you think that google should know about our work.
> --
> Ra
Fields, Christopher J | 14 Feb 05:38 2015

Google Summer of Code 2015


OBF and The Perl Foundation are both working on reapplying for the Google Summer of Code.  Note that we will
gladly accept Perl 6 related projects as well!  

Here is the spot for OBF-specific project ideas:

Perl-related ones can also be added here:

Weigang Qiu | 6 Feb 21:37 2015

Re: Bb-utils

Dear BP developers,

We have packaged the initial release (1.0) of bp-utils, which contains two bioperl-based utilities: "bioseq" and "bioaln" (another two utilities, "biopop" and "biotree", are in pre-release cleaning-up)

For project description, rationale, and call for developers, please see:

Thanks Chris and other for setting up the github access. Please feel free to email me for any feedback and bugs.


Weigang Qiu

On Fri, Dec 12, 2014 at 9:59 PM, Fields, Christopher J <cjfields <at>> wrote:
Hi Weigang,

I wonder whether it would be better to have a separate bputils repo in the BioPerl space.  This would allow development to continue w/o tying it directly to a release, and I think would solve the exposure problem much more so than having it included in the main bioperl-live repo.  We could also feasibly include it as part of the main CPAN bioperl release, maybe by simply linking to it as a git submodule and packaging it up.


On Dec 9, 2014, at 3:01 PM, Weigang Qiu <weigangq <at>> wrote:

Dear BioPerl developers,

I intend to share & (eventually) publish a suite of bioperl-based command-line utilities my lab has developed & found very useful in the past 10 years:

Hilmar encouraged me to get some feedback & advice from the bioperl-L on how to package and release it properly.

My specific questions are:

1. While it goes naturally with bioperl releases, but I would like it to be more exposed than hidden inside the bp-scripts folder. Do you think these utilities are useful enough to be housed in a separate folder (e.g., "bputils") by itself? I have a developer's account but I haven't commit anything for years.

2. How to make it sustainable (and attract new developers) since we are constantly revising and adding methods. It would be great if become a part of bioperl-live.

Any advice and facilitation will be appreciated!

Weigang Qiu, Ph.D.
Department of Biological Sciences
Hunter College of the City University of New York
695 Park Avenue, New York, NY 10065
Office: 1-212-772-5296 (Room 839 Hunter North Building)
Lab: 1-212-772-5721 (Room 830 Hunter North Building)

Bioperl-l mailing list
Bioperl-l <at>

Weigang Qiu
Bioperl-l mailing list
Bioperl-l <at>
Mark A. Jensen | 26 Jan 15:23 2015

Conversion of Phred 33 -> Phred 64 quality

Hi folks,
I know I could RTFM, but maybe someone knows off the top of their head: 
I understand that Illumina at one time made a switch in the constant 
added to quality scores to generate the FASTQ that comes off their 
instruments. This leads to a certain incomparability of data before and 
after that switch. This is about all I know of the issue; does anyone 
here have experience with this? Are there any BP modules that do this 
much appreciated-
Shane McCoy | 14 Jan 00:23 2015

Properly installing Bioperl - Ubuntu - Git

Not too familiar w/ Bioperl and want to make sure i've properly installed on my VM ubuntu 14.4 server w/ git. 

sudo apt-get update
sudo apt-get upgrade
sudo apt-get check
sudo apt-get autoclean

sudo apt-get install build-essential git gcc
sudo apt-get install libexpat1 expat
sudo apt-get install libgd-dev  (unable to locate libgd recommended on Dependencies Outside perl)
sudo apt-get install libssl-dev
sudo apt-get install libpq-dev
sudo apt-get install libdb-dev libperl-dev
sudo apt-get install libgd-gd2-perl
sudo apt-get install libxml2
sudo apt-get install libxml2-dev

Installed local::lib manually for perl modules. Set environment path 
echo '[ $SHLVL -eq 1 ] && eval "$(perl -I$HOME/perl5/lib/perl5 -Mlocal::lib)"' >>~/.bashrc

cpan> install Bundle::CPAN (archive::zip failed? pff!) 
cpan> install Bundle::BioPerl
Failed during this command:
 LDS/AcePerl-1.92.tar.gz                      : make_test NO
 CJFIELDS/Bio-ASN1-EntrezGene-1.70.tar.gz     : make_test NO
 CAPTTOFU/DBD-mysql-4.029.tar.gz              : writemakefile NO '/usr/bin/perl Makefile.PL 
and several more. missing packages etc and successfully installed except the above 3. Don't think i need those.

sudo apt-get install libgd-gd2-perl

cd bioperl-live 
perl Build.PL 
--fails - needs more modules

Configuring C/CJ/CJFIELDS/BioPerl-1.6.924.tar.gz with Build.PL
Checking prerequisites...

HTML::TableExtract is not installed
PostScript::TextBlock is not installed
Net::SSLeay Crypt::SSLeay IO::Socket::SSL
(all above installed via CPAN)

Checking optional features...
    ! Bio::ASN1::EntrezGene is not installed
 - starts installing bioperl 
i let it run and it installed Bioperl fine via CPAN.  build OK. 
Installed bio::asn1::entrezgene & bio::root::build
perl Build.PL
now needs Inline::C...
cpan>install Inline::C
 cd bioperl-live 

 perl Build.PL 
 -- fails. missing Manifest?
Build manifest
perl Build.PL --still fails
perl Build -- works
./Build test 
sudo ./Build install 

Installed / passed tests OK. 
export PERL5LIB="$HOME/bioperl-live:$PERL5LIB"
to .bashrc

So i've installed via CPAN which i didn't really want, but couldn't run Build.PL w/o prereq Bio::ASN1::EntrezGene which starts installing bioperl. What should i have done differently? 
i still did Build w/ Git clone, but when i check Version i get 1.006924. I thought the bioperl-live version was 1.006925? i set the environment variable in bashrc for bioperl-live. 
Ok, thanks for your time! Just want to make sure i'm doing this correctly. 
Bioperl-l mailing list
Bioperl-l <at>
Demian Riccardi | 13 Jan 15:52 2015

HackaMol 0.01 released

Hello BioPerl community!

HackaMol 0.01 has been released.  The library, very much inspired by BioPerl, targets structural analysis and multiscale molecular modeling. 

While this is an early release, I've pushed it as far as I could over the past year. There are surely areas that need work.  Contributions are welcome!  



Bioperl-l mailing list
Bioperl-l <at>
Carnë Draug | 9 Jan 01:04 2015

bioperl executables files and dist-zilla


when installing new bioperl modules with dist zilla, the executable files
are not being installed.  The reason is that our dist zilla plugin bundle [1]
is searching for them in the default directory "bin/" while bioperl seems
to have them in a "scripts/" directory.  This is done by the ExecFiles plugin.

I can change this easily or I can just move the binaries to "bin/" (which
I think is more common).  Also, dropping the file extensions would make
things simpler.

One last thing, I was thinking of adding the SetScriptShebang plugin [3]
to the bundle. Opinions?



Bioperl-l mailing list
Bioperl-l <at>
Zakariyya Mughal | 24 Dec 06:01 2014

A fast and natural interface to R from Perl

Hi everyone,

I have uploaded modules for working with the R interpreter from Perl.
The CPAN links are below, but to get a taste of what the API
looks like, check out my blog post <>.

- Statistics::NiceR <>
- Data::Frame <>

I'd love to have feedback on how to improve them.

Regards and happy hacking,
- Zaki Mughal
Thiago M. Venancio | 17 Dec 19:04 2014

MSG: Undefined sub-sequence

Hi all,
I am parsing a BLASTX report and got the following error for some sequences.
It always happen with sub-sequences of 1 base. I found some posts on this issue, but they are old.
I am use the latest version installed through cpan.
$ perl -MBio::Perl -le 'print Bio::Perl->VERSION;'

Any help is appreciated.

------------- EXCEPTION: Bio::Root::Exception -------------
MSG: Undefined sub-sequence (241,242). Valid range = 96 - 242
STACK: Error::throw
STACK: Bio::Root::Root::throw /usr/share/perl5/Bio/Root/
STACK: Bio::Search::HSP::HSPI::matches /usr/share/perl5/Bio/Search/HSP/
STACK: Bio::Search::SearchUtils::_adjust_contigs /usr/share/perl5/Bio/Search/
STACK: Bio::Search::SearchUtils::tile_hsps /usr/share/perl5/Bio/Search/
STACK: Bio::Search::Hit::GenericHit::frac_aligned_hit /usr/share/perl5/Bio/Search/Hit/

Thiago Motta Venancio, M.Sc., PhD
Bioperl-l mailing list
Bioperl-l <at>