-
Notifications
You must be signed in to change notification settings - Fork 0
/
fasta_to_fastq.pl
46 lines (41 loc) · 1.85 KB
/
fasta_to_fastq.pl
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
34
35
36
37
38
39
40
41
42
43
44
45
46
#Copyright (c) 2010 LUQMAN HAKIM BIN ABDUL HADI ([email protected])
#
#Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files
#(the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify,
#merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is
#furnished to do so, subject to the following conditions:
#The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.
#THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
#OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
#LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR
#IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
#!/usr/bin/perl
use strict;
my $file = $ARGV[0];
open FILE, $file;
my ($header, $sequence, $sequence_length, $sequence_quality);
while(<FILE>) {
chomp $_;
if ($_ =~ /^>(.+)/) {
if($header ne "") {
print "\@".$header."\n";
print $sequence."\n";
print "+"."\n";
print $sequence_quality."\n";
}
$header = $1;
$sequence = "";
$sequence_length = "";
$sequence_quality = "";
}
else {
$sequence .= $_;
$sequence_length = length($_);
for(my $i=0; $i<$sequence_length; $i++) {$sequence_quality .= "I"}
}
}
close FILE;
print "\@".$header."\n";
print $sequence."\n";
print "+"."\n";
print $sequence_quality."\n";