Thiago Venancio | 27 Feb 21:47 2015
Picon

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/Root.pm:449
STACK: Bio::Search::HSP::HSPI::matches_MT /usr/local/share/perl/5.20.1/Bio/Search/Tiling/MapTileUtils.pm:583
STACK: Bio::Search::Tiling::MapTiling::_calc_stats /usr/local/share/perl/5.20.1/Bio/Search/Tiling/MapTiling.pm:1163
STACK: Bio::Search::Tiling::MapTiling::length /usr/local/share/perl/5.20.1/Bio/Search/Tiling/MapTiling.pm:508
STACK: Bio::Search::Tiling::MapTiling::num_aligned /usr/local/share/perl/5.20.1/Bio/Search/Tiling/MapTiling.pm:694
STACK: Bio::Search::Tiling::MapTiling::frac_aligned /usr/local/share/perl/5.20.1/Bio/Search/Tiling/MapTiling.pm:668
STACK: blParse.pl:40
-----------------------------------------------------------

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.
Thiago

_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l
Thiago Venancio | 27 Feb 19:21 2015
Picon

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 ?

Best,
Thiago




! iv>
_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l
Raoul J P Bonnal | 23 Feb 09:51 2015

OBF GSoC Registration

Please all possible mentors,

register yourself

https://www.google-melange.com/gsoc/homepage/google/gsoc2015

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.

Best,
--
Ra & Fra

On 2/20/15 8:51 PM, Raoul Bonnal wrote:
> Dear All,
> I have created a special section
>
>
> http://www.open-bio.org/wiki/Google_Summer_of_Code#Scientific_Achievements
>
> 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

All,

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:

 http://www.open-bio.org/wiki/Google_Summer_of_Code

Perl-related ones can also be added here:

 http://blogs.perl.org/users/shadowcat_mdk/2015/02/gsoc-2015.html
 http://wiki2.enlightenedperl.org/index.php5?title=2015_Ideas_Page

chris
Weigang Qiu | 6 Feb 21:37 2015
Picon

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)

https://github.com/bioperl/bp-utils/blob/master/bp-utils-current-release.tar.gz

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

https://github.com/bioperl/bp-utils/blob/master/README.md

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

cheers,

--
Weigang Qiu


On Fri, Dec 12, 2014 at 9:59 PM, Fields, Christopher J <cjfields <at> illinois.edu> 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.

chris

On Dec 9, 2014, at 3:01 PM, Weigang Qiu <weigangq <at> gmail.com> 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:

http://diverge.hunter.cuny.edu/labwiki/Bioutils

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!

Thanks,
--
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)
Web:http://diverge.hunter.cuny.edu/labwiki/

_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l




--
Weigang Qiu
_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l
Mark A. Jensen | 26 Jan 15:23 2015
Picon

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 
translation?
much appreciated-
MAJ
Shane McCoy | 14 Jan 00:23 2015
Picon

Properly installing Bioperl - Ubuntu - Git

Hello! 
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 gmod.org/wiki/BioPerl 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...
  recommends:

Algorithm::Munkres 
Array::Compare 
Bio::Phylo 
Convert::Binary::C 
GraphViz
HTML::TableExtract is not installed
PostScript::TextBlock is not installed
SOAP::Lite 
Net::SSLeay Crypt::SSLeay IO::Socket::SSL
Sort::Naturally 
XML::LibXML 
DBD::Pg 
DBD::SQLite
(all above installed via CPAN)

Checking optional features...
EntrezGene............disabled
  requires:
    ! 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 
 --PASS
sudo ./Build install 

Installed / passed tests OK. 
Added;
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. 
Shane
_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l
Demian Riccardi | 13 Jan 15:52 2015
Picon

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!  

  

Demian


_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l
Carnë Draug | 9 Jan 01:04 2015
Picon

bioperl executables files and dist-zilla

Hi

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?

Carnë

[1] https://metacpan.org/pod/Pod::Weaver::PluginBundle::BioPerl
[2] https://metacpan.org/pod/Dist::Zilla::Role::ExecFiles
[2] https://metacpan.org/pod/Dist::Zilla::Plugin::SetScriptShebang

_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l
Zakariyya Mughal | 24 Dec 06:01 2014
Picon

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 <http://enetdown.org/dot-plan/posts/2014/12/24/a_fast_and_natural_interface_to_R_from_Perl/>.

- Statistics::NiceR <http://p3rl.org/Statistics::NiceR>
- Data::Frame <http://p3rl.org/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
Picon

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;'
1.006924

Any help is appreciated.
Thiago

------------- 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/Root.pm:472
STACK: Bio::Search::HSP::HSPI::matches /usr/share/perl5/Bio/Search/HSP/HSPI.pm:716
STACK: Bio::Search::SearchUtils::_adjust_contigs /usr/share/perl5/Bio/Search/SearchUtils.pm:431
STACK: Bio::Search::SearchUtils::tile_hsps /usr/share/perl5/Bio/Search/SearchUtils.pm:201
STACK: Bio::Search::Hit::GenericHit::frac_aligned_hit /usr/share/perl5/Bio/Search/Hit/GenericHit.pm:1319

--
=================================
Thiago Motta Venancio, M.Sc., PhD
http://venancio.openwetware.org/
=================================
_______________________________________________
Bioperl-l mailing list
Bioperl-l <at> mailman.open-bio.org
http://mailman.open-bio.org/mailman/listinfo/bioperl-l

Gmane