GenePrimp error: Illegal division by zero at BPlite.pm line 283, <$__ANONIO__> line 400764. · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

GenePrimp error: Illegal division by zero at BPlite.pm line 283, <$__ANONIO__> line 400764.

Hello, Did anybody has the error while the GenePrimp (http://www.nature.com/nmeth/journal/v7/n6/full/nmeth.1457.html) parsing the Blastp output at the first stage of analysis. I have the same issue by analysis both the gbk file from ncbi and the gbk file produced by Prokka. I also tested the different versions of Perl 5.8.5, 5.16, 5.24. I can't understand what is the problem. As I remember, a year ago everything was OK. I also found that analysis of fragmented prokaryotic genome with annotation was successful, but using complete genome (~1 Mb) fails.

Thank you in advance, for any suggestions.

software error gene genome

• 2,345 views • link

written 9.4 years ago by konstantinkul

I've just surfed my previous posts in BioStar and decided to explain how I solved this issue a year ago.

The problem was due to incorrect formatting of last releases of protein (nr) blast databases. When I extracted all sequence from NCBI blast database and again created the blast database but with old version makeblastdb. After this reformatting GenePrimp works fine.

• 1 views • link

written 7.8 years ago by konstantinkul

0 answers

No answers yet.

Log in to answer this question.

More posts like this

sc-WBGS library GC content, low bismark mapping rate, high CHH%
written 9.0 years ago by Jingyue •

Hi, I have met some problem in analysis sc-WGBS dataset. The single cell WGBS experimental protocol was followed the nature protocol: [Smallwood et.al][1]. This is …
Tool: MAGERI: a software tool for calling rare variants and detecting circulating tumor DNA from UM…
written 9.3 years ago by mikhail.shugay

Dear Colleagues, I would like to announce our recently published software tool called [MAGERI][1] that is designed to facilitate the detection of ultra-rare variants from …
For only a few genes, Cuffdiff report more than one FPKM values per locus. Why?
written 9.4 years ago by rice.researcher

I followed similar steps mentioned in the procedure part of [Trapnell et al, 2012][1] for an RNASeq analysis of oryza sativa datasets. The problem I …
(Modern, mid-2016) RNA-seq software pipeline
written 10.2 years ago by jp

I am trying to follow the RNA-seq protocol described in (this excellent paper): Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and …
How to get the alignment summary of HISAT2
written 10.2 years ago by Wet&DryImmunology

Hi, I was trying to HISAT2, which is suppose to be the improved version of tophat or tophat2, [HISAT publication][1] , I applied it to …
Tophat Protocol - nature paper
written 11.7 years ago by David_emir

Hi All, I was doing differential gene and transcript expression analysis by following the Tophat protocol ([published in Nature protocol in March 2012][1]) on the …
AHRD functional annotation program (error)
written 11.7 years ago by jomaco

Hi, I am carrying out functional annotation of predicted genes in our genome using [AHRD][1]. [1]: https://github.com/groupschoof/AHRD [2]: https://www.biostars.org/p/80714/ I have the same problem as …
Tool: Dgidb - Mining The Druggable Genome
written 12.8 years ago by Obi Griffith

The druggable genome can be defined as the genes or gene products that are known or predicted to interact with drugs, ideally with a therapeutic …
Tool: Dgidb - Mining The Druggable Genome
written 12.8 years ago by Obi Griffith

The druggable genome can be defined as the genes or gene products that are known or predicted to interact with drugs, ideally with a therapeutic …
Tutorial: Rna-Seq Review Papers
written 13.9 years ago by Istvan Albert

## Review papers on the topic of RNA-seq ### General Guides - [Computational methods for transcriptome annotation and quantification using RNA-seq (May 2011)][1] in Nature …