[ViewVC] Annotation of: gclib/scripts/fq2fa

#!/usr/bin/perl
use strict;
while (<>) {
 if (m/^@(\S+)/) {
   my $seqname=$1;
   my $seq=<>;
   $seq=~tr/\n\r//d;
   my $qh=<>;
   #chomp($qh);
   $qh=~tr/\n\r//d;
   die("Error: 3rd line of $seqname must start with '+' !\n")  
      unless ($qh eq '+') || ($qh eq '+'.$seqname);
   my $quals=<>;
   print ">$seqname\n$seq\n";
   }
}
Revision:	23
Committed:	Tue Jul 26 21:44:38 2011 UTC (13 years, 2 months ago) by gpertea
Original Path:	*ann_bin/fq2fa*
File size:	335 byte(s)
Log Message:	adding misc scripts
Line	User	Rev	File contents
1	gpertea	23	#!/usr/bin/perl
2			use strict;
3			while (<>) {
4			if (m/^@(\S+)/) {
5			my $seqname=$1;
6			my $seq=<>;
7			$seq=~tr/\n\r//d;
8			my $qh=<>;
9			#chomp($qh);
10			$qh=~tr/\n\r//d;
11			die("Error: 3rd line of $seqname must start with '+' !\n")
12			unless ($qh eq '+') \|\| ($qh eq '+'.$seqname);
13			my $quals=<>;
14			print ">$seqname\n$seq\n";
15			}
16			}