comparison parse_primersearch.pl @ 5:b321e0517be3 draft

Uploaded
author ben-warren
date Thu, 22 May 2014 20:30:19 -0400
parents 21053f7f9ed1
children
comparison
equal deleted inserted replaced
4:be070a68521e 5:b321e0517be3
1 #!/usr/bin/perl 1 #!/usr/bin/perl
2 #parse_primersearch.pl 2 #parse_primersearch.pl
3 #reformat EMBOSS primersearch output into columnar Galaxy interval format 3 #reformat EMBOSS primersearch output into columnar Galaxy interval format
4 4 #usage perl parse_primersearch.pl <input primersearch file> <tabular output file>
5 #Copyright 2012 John McCallum 5 #Copyright 2012 John McCallum
6 #New Zealand Institute for Plant and Food Research 6 #New Zealand Institute for Plant and Food Research
7 #This program is free software: you can redistribute it and/or modify 7 #This program is free software: you can redistribute it and/or modify
8 # it under the terms of the GNU General Public License as published by 8 # it under the terms of the GNU General Public License as published by
9 # the Free Software Foundation, either version 3 of the License, or 9 # the Free Software Foundation, either version 3 of the License, or
18 # along with this program. If not, see <http://www.gnu.org/licenses/>. 18 # along with this program. If not, see <http://www.gnu.org/licenses/>.
19 19
20 open (IN, "<$ARGV[0]"); 20 open (IN, "<$ARGV[0]");
21 open (OUT, ">$ARGV[1]"); 21 open (OUT, ">$ARGV[1]");
22 22
23 #print OUT "primerset_id","\t","sequence_id","\t","hit_start","\","mismatches","\t","amplimer_size",\n"; 23
24 24
25 25
26 26
27 while (<IN>) { 27 while (<IN>) {
28 /^Primer name (\S+)/ && ($name = $1); # get primer set name 28 /^Primer name (\S+)/ && ($name = $1); # get primer set name