annotate pyPRADA_1.2/tools/bwa-0.5.7-mh/qualfa2fq.pl @ 0:acc2ca1a3ba4

Uploaded
author siyuan
date Thu, 20 Feb 2014 00:44:58 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
1 #!/usr/bin/perl -w
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
2
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
3 use strict;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
4 use warnings;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
5
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
6 die("Usage: qualfa2fq.pl <in.fasta> <in.qual>\n") if (@ARGV != 2);
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
7
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
8 my ($fhs, $fhq, $q);
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
9 open($fhs, ($ARGV[0] =~ /\.gz$/)? "gzip -dc $ARGV[0] |" : $ARGV[0]) || die;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
10 open($fhq, ($ARGV[1] =~ /\.gz$/)? "gzip -dc $ARGV[1] |" : $ARGV[1]) || die;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
11
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
12 $/ = ">"; <$fhs>; <$fhq>; $/ = "\n";
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
13 while (<$fhs>) {
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
14 $q = <$fhq>;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
15 print "\@$_";
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
16 $/ = ">";
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
17 $_ = <$fhs>; $q = <$fhq>;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
18 chomp; chomp($q);
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
19 $q =~ s/\s*(\d+)\s*/chr($1+33)/eg;
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
20 print $_, "+\n";
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
21 for (my $i = 0; $i < length($q); $i += 60) {
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
22 print substr($q, $i, 60), "\n";
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
23 }
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
24 $/ = "\n";
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
25 }
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
26
acc2ca1a3ba4 Uploaded
siyuan
parents:
diff changeset
27 close($fhs); close($fhq);