From Code to Community: Sponsoring The Perl and Raku Conference 2025 Learn more

LRR / Polloc-v1.5.3 / scripts / polloc_gel.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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96

              
#!/usr/bin/perl
use strict;
use Bio::Polloc::LocusIO 1.0501;
use Bio::Polloc::Genome;
use Bio::Polloc::TypingI;
use Pod::Usage;
# ------------------------------------------------- INPUT
my $gff_in =  shift @ARGV;
my $out    =  shift @ARGV;
my @names  = split /:/, shift @ARGV;
pod2usage(1) unless $gff_in and $out and $#names > -1;
Bio::Polloc::Polloc::Root->DEBUGLOG(-file=>">$out.log");
$Bio::Polloc::Polloc::Root::VERBOSITY = 4;
# ------------------------------------------------- READ INPUT
my $genomes = [];
for my $G (0 .. $#names){
   push @$genomes, Bio::Polloc::Genome->new(-name=>$names[$G], -id=>$G) }
my $LocusIO = Bio::Polloc::LocusIO->new(-file=>$gff_in);
my $inloci = $LocusIO->read_loci(-genomes=>$genomes);
# ------------------------------------------------- TYPING
my $typing = Bio::Polloc::TypingI->new(-type=>'bandingPattern::amplification');
open IMG, ">", "$out.png" or die "I can not open '$out.png': $!\n";
binmode IMG;
print IMG $typing->graph(-locigroup=>$inloci)->png;
close IMG;
__END__
=pod
=head1 AUTHOR
Luis M. Rodriguez-R < lmrodriguezr at gmail dot com >
=head1 DESCRIPTION
polloc_gel.pl - Draws the expected gel for a given set of amplicons.
=head1 LICENSE
This script is distributed under the terms of
I<The Artistic License>.  See LICENSE.txt for details.
=head1 SYNOPSIS
C<perl polloc_gel.pl> B<arguments>
The arguments must be in the following order:
=over
=item Input gff
GFF3 file containing the amplicons.
Example: C</tmp/polloc-primers.out.amplif.1.gff>
=item Output
Path to the base of the output files.
Example: C</tmp/polloc-gel.out>
=item Names
The names of the genomes separated by colons (C<:>). Alternatively, it
can be an empty string (C<''>) to assign genome names from files.
Example: C<Xci3:Xeu8:XamC>
=back
Run C<perl polloc_gel.pl> without arguments to see the help
message.
=head1 SEE ALSO
=over
=item *
L<Bio::Polloc::Locus::amplicon>
=item *
L<Bio::Polloc::Typing::bandingPattern>
=back
=cut