Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/env perl
- use strict;
- use warnings FATAL => 'all';
- use Data::Dumper::Concise;
- use Log::Log4perl ':easy';
- use Perl6::Slurp;
- use HTGT::DBFactory;
- Log::Log4perl->easy_init;
- my $file = shift || die "Usage: $0 <targ_vec_names>";
- my $conn = HTGT::DBFactory->connect('eucomm_vector');
- my @targ_vec_names = slurp $file, { chomp => 1 };
- my %mgi_accession_ids = ();
- for my $targ_vec_name (@targ_vec_names) {
- if ( $targ_vec_name =~ m/^(\w+)_(\w+)$/o ) {
- my $targ_vec_rs = $conn->resultset('WellSummaryByDI')->search_rs(
- { pgdgr_plate_name => $1,
- pgdgr_well_name => $2,
- pgdgr_distribute => 'yes'
- },
- { distinct => 1 }
- );
- while ( my $targ_vec = $targ_vec_rs->next ) {
- $mgi_accession_ids{$targ_vec_name}
- = $targ_vec->project->mgi_gene->mgi_accession_id;
- }
- }
- }
- INFO scalar keys %mgi_accession_ids, '/', scalar @targ_vec_names,
- ' are in fact distributable.';
- exit 0;
Add Comment
Please, Sign In to add comment