AutoAug.pl running with PASA+GMAP

Discussions about training AUGUSTUS from various sources of evidence. Not discussed here: BRAKER1 and WebAUGUSTUS!

Moderator: bioinf

Post Reply
User avatar
katharina
Site Admin
Posts: 531
Joined: Wed Nov 18, 2015 6:14 pm
Location: Greifswald
Contact:

AutoAug.pl running with PASA+GMAP

Post by katharina »

Originally posted in the old forum by Vickie on 18.01.2013 - 13:25
Hi,
I'm writing this post because the AutoAug.pl README file is a little bit out of date. An exactly that part which advises to delete -S option from process_gmap_alignments.pl . Because of that, new PASA+GMAP starts complaining on duplicate entry, and asks to rerun pipline at index 4.
The exact error message is:

Code: Select all

DBD::mysql::db do failed: Duplicate entry 'gmap-comp0_c0_seq1_Sample1' for key 2 at /home/user/species/PASA_r2012-06-25/PerlLib//Mysql_connect.pm line 155, <BTAB> line 3. 
Undefined subroutine &Scalar::Util::blessed called at /usr/lib/perl5/5.8.8/overload.pm line 89, <BTAB> line 3.
And error appears because the gmap.tophits.btab file contains duplicates. Example:

Code: Select all

contig642325 gmap comp0_c0_seq1_Sample1 294 546 
1 253 99 1 1 
contig642325 gmap comp0_c0_seq1_Sample1 421 540 
254 373 99 1 2 
scaffold25491 gmap comp100019_c0_seq1_Sample1 8575 
6355 62 2284 99 12 2 
scaffold25491 gmap comp100019_c0_seq1_Sample1 3901 
3815 2285 2371 100 12 3
So when I returned -S flag back, and the line again looked like:

Code: Select all

my $cmd = "gmap_setup -S -D gmap_db_dir -d $genome_db $genome_db"; 
it started work properly without any errors.
I used PASA_r2012-06-25 and GMAP-2012-04-07 (Don't use GMAP beta ver. because it causes segmentation fault).
And one more thing, when using MySQL without root priveleges, use should first create database for pasa (for example:dbname-PASA , user priv> read and write), and write it into conf.txt (as said in PASA readme) and then create a new DB with read and write priveleges with name PASAspeciesname, where species name is the same name you use in --species=speciesname. Then there would be no poblems.
Hope this information would be somehow useful.
Vickie
User avatar
katharina
Site Admin
Posts: 531
Joined: Wed Nov 18, 2015 6:14 pm
Location: Greifswald
Contact:

Re: AutoAug.pl running with PASA+GMAP

Post by katharina »

by katharina on 21.01.2013 - 10:30
Hi Vickie,
thanks for telling us. I have included your comments in the README of autoAug.pl.
It might take a little while until it takes effect in the sources that are available for download since Mario has to make a new AUGUSTUS release.
Best,
Katharina
Post Reply