Mercurial > repos > edward-kirton > roche454_toolsuite
diff sff_to_fastq_converter.pl @ 2:2d86d5b112e8
Uploaded
author | edward-kirton |
---|---|
date | Thu, 14 Jul 2011 22:14:07 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/sff_to_fastq_converter.pl Thu Jul 14 22:14:07 2011 -0400 @@ -0,0 +1,31 @@ +#!/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;