--
View this message in context: http://old.nabble.com/Convert-fastq-to-fasta-tp31492543p31492543.html
Sent from the Perl - Bioperl-L mailing list archive at Nabble.com.
_______________________________________________
Bioperl-l mailing list
Biop...@lists.open-bio.org
http://lists.open-bio.org/mailman/listinfo/bioperl-l
try fastq_to_fasta part of http://hannonlab.cshl.edu/fastx_toolkit/
usage: fastq_to_fasta [-h] [-r] [-n] [-v] [-z] [-i INFILE] [-o OUTFILE]
Malcolm Cook
Stowers Institute for Medical Research - Bioinformatics
Kansas City, Missouri USA
------------- transform.pl ---------------
use Bio::SeqIO;
my ($file1,$file2)=@ARGV;
my $seqin = Bio::SeqIO -> new (-format => 'fastq',-file => $file1);
my $seqout = Bio::SeqIO -> new (-format => 'fasta',-file => ">$file2");
while (my $seq_obj = $seqin -> next_seq)
{
$seqout -> write_seq($seq_obj);
}
-------------------------------------------
run as this:
perl transfrom.pl sequence.fastq sequence.fasta
> Message: 7
> Date: Wed, 27 Apr 2011 17:26:49 -0700 (PDT)
> From: perlbio007 <perlb...@gmail.com>
> Subject: [Bioperl-l] Convert fastq to fasta
> To: Biop...@lists.open-bio.org
> Message-ID: <314925...@talk.nabble.com>
> Content-Type: text/plain; charset=us-ascii
>
>
> Iam new to Bioperl. Pls help.
> I have a zip folder of sequences which is in fastq format. I need to
> convert
> it in fasta format?
> How I do that using bioperl?What module do I need?
>
--
Tao Zhu, College of Life Sciences, Beijing Normal University, Beijing
100875, China
Email: tz...@mail.bnu.edu.cn
Website: http://bnuzt.org (mainly written in Chinese)
On 28/04/11 10:26, perlbio007 wrote:
> Iam new to Bioperl. Pls help.
> I have a zip folder of sequences which is in fastq format. I need to convert
> it in fasta format?
> How I do that using bioperl?What module do I need?
_______________________________________________