[Bioperl-l] Newbie question on Bio::SeqIO

Chris Fields cjfields at illinois.edu
Wed Jan 26 23:07:51 UTC 2011


Also, FASTQ parsing only works correctly if you are using BioPerl 1.6.1.  Not sure what version fink currently has.

chris

On Jan 26, 2011, at 4:30 PM, Adam Sjøgren wrote:

> On Wed, 26 Jan 2011 14:05:08 -0800 (PST), JayPea wrote:
> 
>> #!/usr/bin/perl
>> use Bio::Perl;
>> use Bio::SeqIO;
> 
>> my $seqio_obj = Bio::SeqIO->new(-file => "dna.fa",
>>                                -format => 'Fasta');
>> my $output = $seqio_obj->subseq(1,20);
> 
> Try replacing that line by these two:
> 
> my $seq=$seqio_obj->next_seq;
> my $output=$seq->subseq(1,20);
> 
>> print "$output\n";
> 
> See "perldoc Bio::SeqIO" :-)
> 
> 
> Also: Never forget to start your scripts with "use strict; use
> warnings;" - a great way to avoid typos and stuff.
> 
> 
>  Best regards,
> 
>    Adam
> 
> -- 
> "Här kommer rädslan, gamle vän                               Adam Sjøgren
>  När alla fjärilar i magen vaknar upp                   asjo at koldfront.dk
>  Viskar välkommen hem"
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l





More information about the Bioperl-l mailing list