[Bioperl-l] seqIO.pm

Smithies, Russell Russell.Smithies at agresearch.co.nz
Mon Nov 21 20:04:01 UTC 2011


Or you could use the builtin script bp_sreformat.pl

--Russell

> -----Original Message-----
> From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-
> bounces at lists.open-bio.org] On Behalf Of Rondon Neto
> Sent: Tuesday, 22 November 2011 6:31 a.m.
> Cc: bioperl-l at bioperl.org
> Subject: Re: [Bioperl-l] seqIO.pm
> 
> Hi! try this script:
> 
> #!/usr/bin/perl
> use warnings;
> use strict;
> use Bio::SeqIO;
> 
> if (!$ARGV[0]) { die "\n### USAGE::: perl $0 [file.fastq]\n"; }
> 
> my $fastq = $ARGV[0];
> 
> my $in = Bio::SeqIO->new( -file => $fastq,
>                           -format => 'fastq' );
> 
> my $out = Bio::SeqIO->new ( -file => ">out.phd",
>                             -format=> 'phd');
> 
> while (my $seq = $in->next_seq()) {
>       $out->write_seq($seq);
> }
> 
> exit;
> 
> 
> Best wishes,
> Rondon, a brazilian friend.
> 
> 
> 
> 
> 
> 
> ________________________________
>  De: Peter Cock <p.j.a.cock at googlemail.com>
> Para: Nathan Taylor <zntayl at gmail.com>
> Cc: bioperl-l at bioperl.org
> Enviadas: Segunda-feira, 21 de Novembro de 2011 14:15
> Assunto: Re: [Bioperl-l] seqIO.pm
> 
> On Thu, Nov 17, 2011 at 1:07 AM, Nathan Taylor <zntayl at gmail.com> wrote:
> > Hello,
> >
> >   Can SeqIO.pm convert a file of fastq reads into .phd files. Or,
> > barring that, a file of fastas and file of quals into .phd files?
> >
> > Many thanks,
> > Nathan
> 
> In principle that is possible (e.g. Biopython can do fastq to phd).
> Have you tried using BioPerl's SeqIO to do this? Was there an error message?
> 
> Peter
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
=======================================================================
Attention: The information contained in this message and/or attachments
from AgResearch Limited is intended only for the persons or entities
to which it is addressed and may contain confidential and/or privileged
material. Any review, retransmission, dissemination or other use of, or
taking of any action in reliance upon, this information by persons or
entities other than the intended recipients is prohibited by AgResearch
Limited. If you have received this message in error, please notify the
sender immediately.
=======================================================================




More information about the Bioperl-l mailing list