Mercurial > repos > edward-kirton > roche454_toolsuite
view sff_to_fastq_converter.pl @ 3:bf1f8bc4abe6 default tip
minor updates for v2.6
author | eskirton@lbl.gov |
---|---|
date | Wed, 21 Dec 2011 19:42:53 -0800 |
parents | 2d86d5b112e8 |
children |
line wrap: on
line source
#!/usr/bin/env perl use warnings; use strict; use Getopt::Long; use File::Basename; # VALIDATE INPUT die("Expected 3 args") unless @ARGV == 3; my ($sff, $extra_files_path, $fastq) = @ARGV; # DEFINE PATHS mkdir($extra_files_path) unless -d $extra_files_path; my $base = basename($sff); my $fasta = "$extra_files_path/$base.fasta"; my $qual = "$extra_files_path/$base.qual"; # GENERATE FASTA, QUAL, FASTQ my $outf; my $out; eval { $out=`sffinfo -seq $sff > $fasta` }; die("ERROR: $out") if $@; print $out; eval { $out=`sffinfo -qual $sff > $qual` }; die("ERROR: $out") if $@; print $out; eval { $out=`fasta_qual_to_fastq $fasta $qual $fastq` }; die("ERROR: $out") if $@; print $out; unlink($fasta, $qual); exit 0;