Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bio::Restriction::Analysis::positions does not return anything #2

Open
phlatphish opened this issue Nov 15, 2018 · 1 comment
Open

Comments

@phlatphish
Copy link

use strict;
use Bio::Seq;
use Bio::SeqIO;
use Bio::Restriction::Analysis;


my $fasta_file = "sample.fa";

my $inputstream = Bio::SeqIO->new(-file => $fasta_file,-format => 'Fasta');

my $all_collection = Bio::Restriction::EnzymeCollection->new();
my $apeki = $all_collection->get_enzyme('TseI');

my $seqobj = $inputstream->next_seq();
my $name = $seqobj->display_name();   

my $restriction_analysis = Bio::Restriction::Analysis->new(-seq=>$seqobj);

my @positions = $restriction_analysis->positions($apeki); 
print STDERR "$name: @positions\n";  ### Nothing

my @sizes = $restriction_analysis->sizes($apeki); 
print STDERR "$name: @sizes\n";  ### List of frag sizes

Like:

000000F: 
000000F: 2496 86 760 1050 1250 1219 1076 1942 2317 4755 1181 400 5746 523 80 692 656 2408 2663 871 75 797 240 375 300 1165 310 602 518 226 3679 2101 68 119 362 1819 1821 570 946 3568 5068 734 117 352 17 39 690 339 1687 44 1641 3265 145 108 1852 72 834 3041 290 1051 3724 1501 5038 523 51 108 762 484 9 1165 392 528 1063 1122 1675 638 323 143 23 1252 252 1365 244 56 2954 915 471

BioPerl conda installation:

perl-bioperl 1.7.2 pl526_8

file name : perl-bioperl-1.7.2-pl526_8.tar.bz2
name : perl-bioperl
version : 1.7.2
build string: pl526_8
build number: 8
channel : bioconda
size : 11 KB
arch : x86_64
has_prefix : False
license : perl_5
md5 : 6cc0ed7cb1f1b53a17f85e31f58d2f9d
noarch : None
platform : linux
requires : ()
subdir : linux-64
timestamp : 1533830894938
url : https://conda.anaconda.org/bioconda/linux-64/perl-bioperl-1.7.2-pl526_8.tar.bz2
dependencies:
perl >=5.26.2,<5.26.3.0a0
perl-bio-asn1-entrezgene
perl-bio-coordinate
perl-bio-featureio
perl-bio-samtools
perl-bio-tools-phylo-paml
perl-bio-tools-run-alignment-clustalw
perl-bio-tools-run-alignment-tcoffee
perl-bioperl-core 1.7.2
perl-bioperl-run

@carandraug
Copy link
Member

For what is worth, Bio::Restriction::Analysis has been moved to its own distribution https://github.com/bioperl/Bio-Restriction (which actually needs a release).

@carandraug carandraug transferred this issue from bioperl/bioperl-live Dec 5, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants