/usr/share/bwa/qualfa2fq.pl is in bwa 0.7.12-5.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
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 | #!/usr/bin/perl -w
use strict;
use warnings;
die("Usage: qualfa2fq.pl <in.fasta> <in.qual>\n") if (@ARGV != 2);
my ($fhs, $fhq, $q);
open($fhs, ($ARGV[0] =~ /\.gz$/)? "gzip -dc $ARGV[0] |" : $ARGV[0]) || die;
open($fhq, ($ARGV[1] =~ /\.gz$/)? "gzip -dc $ARGV[1] |" : $ARGV[1]) || die;
$/ = ">"; <$fhs>; <$fhq>; $/ = "\n";
while (<$fhs>) {
$q = <$fhq>;
print "\@$_";
$/ = ">";
$_ = <$fhs>; $q = <$fhq>;
chomp; chomp($q);
$q =~ s/\s*(\d+)\s*/chr($1+33)/eg;
print $_, "+\n";
for (my $i = 0; $i < length($q); $i += 60) {
print substr($q, $i, 60), "\n";
}
$/ = "\n";
}
close($fhs); close($fhq);
|