[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:	24
Committed:	Tue Jul 26 21:46:39 2011 UTC (13 years, 6 months ago) by gpertea
File size:	335 byte(s)
Log Message:
Line	File contents
1	#!/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	}