Category: BioPerl

Rewriting perl code using bioperl

Code: #!/usr/bin/perl use strict; use warnings; die “Usage: $0 <fasta_file style=”box-sizing: border-box; color: rgb(102, 102, 102); font-family: ‘Open Sans’, sans-serif; font-size: 13px; font-style: italic; line-height: 20px; background-color: rgb(251, 251, 251);”><coord_file style=”box-sizing: border-box;”>n” unless @ARGV > 0; my ($fasta_file, $coord_file) = @ARGV; open FASTA, “<” . $fasta_file; my $seq_id; my $dna_string =…

Continue Reading Rewriting perl code using bioperl

Any tools converting Genbank format to GFF3 format?

Any tools converting Genbank format to GFF3 format? 4 Dear all, As my title describe, I am asking help to convert Genbank format to GFF format.  By GOOGLing, I found a perl script bp_genbank2gff3.pl), which has many bioperl dependencies.  Could anyone inform me other easy-to-use tools?  Any suggestions will be…

Continue Reading Any tools converting Genbank format to GFF3 format?

biopython extract sequence from fasta

My two questions are: What is the simplest way to do this? This unique book shows you how to program with Python, using code examples taken directly from bioinformatics. using python-bloom-filter, just replace the set with seen = BloomFilter(max_elements=10000, error_rate=0.001). This book is suitable for use as a classroom textbook,…

Continue Reading biopython extract sequence from fasta

Tools To Calculate Average Coverage For A Bam File?

Tools To Calculate Average Coverage For A Bam File? 12 I would like to get the average coverage of all the captured bases in a bam file. What would be the best way to do this? What I am looking is a simple one number like 40X. Given that there…

Continue Reading Tools To Calculate Average Coverage For A Bam File?

Line length limit on input FASTA file: 65,536 characters (limit imposed by bioperl)

Hello, I’m trying to run the following command: agat_sp_extract_sequences.pl -g JU2526_Y39G10AR.22.gff -f JU2526*_region.fa -p And it throws the following error: ————- EXCEPTION: Bio::Root::Exception ————- MSG: Each line of the file must be less than 65,536 characters. Line 2 is 67824 chars. STACK: Error::throw STACK: Bio::Root::Root::throw /home/lgs6452/.conda/envs/exonerate_env/lib/site_perl/5.26.2/Bio/Root/Root.pm:447 STACK: Bio::DB::IndexedBase::_check_linelength /home/lgs6452/.conda/envs/exonerate_env/lib/site_perl/5.26.2/Bio/DB/IndexedBase.pm:757 STACK:…

Continue Reading Line length limit on input FASTA file: 65,536 characters (limit imposed by bioperl)

Accepted bioperl-run 1.7.3-7 (source) into unstable

—–BEGIN PGP SIGNED MESSAGE—– Hash: SHA512 Format: 1.8 Date: Tue, 14 Sep 2021 18:52:51 +0200 Source: bioperl-run Architecture: source Version: 1.7.3-7 Distribution: unstable Urgency: medium Maintainer: Debian Med Packaging Team <debian-med-packag…@lists.alioth.debian.org> Changed-By: Étienne Mollier <emoll…@debian.org> Changes: bioperl-run (1.7.3-7) unstable; urgency=medium . * d/watch: fix broken link to github and adjust…

Continue Reading Accepted bioperl-run 1.7.3-7 (source) into unstable

Difficulty installing Bioperl for use with GUIDANCE

Difficulty installing Bioperl for use with GUIDANCE 0 Hi, I am trying to use the Guidance software to analyze a set of 427 sequences. I have the software downloaded and installed the modules indicated in the Guidance user guide (Perl/BioPerl/Ruby), but when I attempt to run Guidance with the following…

Continue Reading Difficulty installing Bioperl for use with GUIDANCE

Can’t locate Bio/SeqIO.pm – path error?

Can’t locate Bio/SeqIO.pm – path error? 0 Hi, I am trying to run snippy on my UBUNTU18 (inside a conda env that has this tool installed) but am getting an error: Can’t locate Bio/SeqIO.pm in @INC (you may need to install the Bio::SeqIO.pm module) I checked the posts here, and…

Continue Reading Can’t locate Bio/SeqIO.pm – path error?

Ubuntu Manpage: Bio::Procedural – Simple low-dependency procedural interfaces to BioPerl

Provided by: libbio-procedural-perl_1.7.4-2_all NAME Bio::Procedural – Simple low-dependency procedural interfaces to BioPerl VERSION version 1.7.4 SYNOPSIS See Bio::Perl for examples. DESCRIPTION The code in this distribution focuses on simple low-dependency procedural interfaces to BioPerl that don’t require knowledge of BioPerl classes and methods. FEEDBACK Mailing lists User feedback is an…

Continue Reading Ubuntu Manpage: Bio::Procedural – Simple low-dependency procedural interfaces to BioPerl

Bioperl SeqIO.pm cannot be found

Hello, I am trying to run this command line in the terminal from Mac OS, which is copied from a protocol: FEELnc_codpot.pl –i <unannotated_transcript.fa> -a <mRNA_sequence.fa> -g <reference_genome.fa> -m ’shuffle’ –o <unannotated_transcript> However, I get the following error: command not found, probably because I am not using perl, so I…

Continue Reading Bioperl SeqIO.pm cannot be found

Bioperl SeqIO.pm cannot be found

Hello, I am trying to run this command line in the terminal from Mac OS, which is copied from a protocol: FEELnc_codpot.pl –i <unannotated_transcript.fa> -a <mRNA_sequence.fa> -g <reference_genome.fa> -m ’shuffle’ –o <unannotated_transcript> However, I get the following error: command not found, probably because I am not using perl, so I…

Continue Reading Bioperl SeqIO.pm cannot be found

Bioperl SeqIO.pm cannot be found

Hello, I am trying to run this command line in the terminal from Mac OS, which is copied from a protocol: FEELnc_codpot.pl –i <unannotated_transcript.fa> -a <mRNA_sequence.fa> -g <reference_genome.fa> -m ’shuffle’ –o <unannotated_transcript> However, I get the following error: command not found, probably because I am not using perl, so I…

Continue Reading Bioperl SeqIO.pm cannot be found