Guest User

Untitled

a guest
Jul 16th, 2018
104
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.01 KB | None | 0 0
  1. #!/usr/bin/env perl
  2. use strict;
  3. use warnings FATAL => 'all';
  4. use Data::Dumper::Concise;
  5. use Log::Log4perl ':easy';
  6. use Perl6::Slurp;
  7. use HTGT::DBFactory;
  8.  
  9. Log::Log4perl->easy_init;
  10.  
  11. my $file = shift || die "Usage: $0 <targ_vec_names>";
  12. my $conn = HTGT::DBFactory->connect('eucomm_vector');
  13. my @targ_vec_names = slurp $file, { chomp => 1 };
  14. my %mgi_accession_ids = ();
  15.  
  16. for my $targ_vec_name (@targ_vec_names) {
  17. if ( $targ_vec_name =~ m/^(\w+)_(\w+)$/o ) {
  18. my $targ_vec_rs = $conn->resultset('WellSummaryByDI')->search_rs(
  19. { pgdgr_plate_name => $1,
  20. pgdgr_well_name => $2,
  21. pgdgr_distribute => 'yes'
  22. },
  23. { distinct => 1 }
  24. );
  25. while ( my $targ_vec = $targ_vec_rs->next ) {
  26. $mgi_accession_ids{$targ_vec_name}
  27. = $targ_vec->project->mgi_gene->mgi_accession_id;
  28. }
  29. }
  30. }
  31.  
  32. INFO scalar keys %mgi_accession_ids, '/', scalar @targ_vec_names,
  33. ' are in fact distributable.';
  34.  
  35. exit 0;
Add Comment
Please, Sign In to add comment