1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
open(FASTAQ,$ARGV[0]) or die;
my $entry="";
my $i=0;
while(my $line=<FASTAQ>){
#We need to process a new entry
if($i==4){
my ($top,$sequence,$top2,$qualsequence) = $entry =~ /^
\@?(.+?)\n
([^\@]*?)\n
\+?(.+?)\n
(.*)\n
/xs;
print ">$top\n";
print $sequence."\n";
$entry=$line;
$i=0;
}else{
$entry.=$line;
}
$i++;
}
#parsing the last secuence
my ($top,$sequence,$top2,$qualsequence) = $entry =~ /^
\@?(.+?)\n
([^\@]*?)\n
\+?(.+?)\n
(.*)\n
/xs;
print ">$top\n";
print $sequence."\n";
|
