Guest User

Untitled

a guest
Jun 18th, 2018
73
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.96 KB | None | 0 0
  1. #! /usr/bin/env perl
  2. use strict;
  3. use warnings FATAL => 'all';
  4. use BioUtils 'get_all';
  5. use Data::Dumper::Concise;
  6. use File::Basename;
  7. use List::MoreUtils 'uniq';
  8.  
  9. my %cassette_vector_name_file_lookup = (
  10. L1L2_6XOspnEnh_Bact_P =>
  11. "pL1L2_GTIRES_6xostpn_BetactP_FLFL_Map_Verified.gbk",
  12. L1L2_Bact_EM7 => "FINAL_pL1L2_GTIRES_BetactP_FLFL_EM7_neo.gbk",
  13. L1L2_Bact_P => "pL1L2_GTIRES_BetactP_FLFL_MapVerified.gbk",
  14. L1L2_Del_BactPneo_FFL => "pL1L2_Del_BactPneo_FFL_PG.gbk",
  15. L1L2_NorCOMM => "L1L2_NorCOMM.gbk",
  16. L1L2_Pgk_P => "pL1L2_GTIRES_PgkP_FLFL_Verified_Map.gbk",
  17. L1L2_Pgk_PM => "pL1L2_GTIRES_PgkP_FLFL_Neo_Mutant_Verified_Map.gbk",
  18. L1L2_gt0 => "pL1L2_GT0_EUCOMM_MapVerified.gbk",
  19. L1L2_gt1 => "pL1L2_GT1_EUCOMM_MapVerified.gbk",
  20. L1L2_gt2 => "pL1L2_GT2_EUCOMM_MapVerified.gbk",
  21. L1L2_gtk => "pL1L2_GTK_EUCOMM_MapVerified.gbk",
  22. L1L2_hubi_P => "pL1L2_GTIRES_hubiqui_P_FLFL_Map_Verified.gbk",
  23. L1L2_st0 => "pL1L2_ST0_EUCOMM_MapVerified.gbk",
  24. L1L2_st1 => "pL1L2_ST1_EUCOMM_MapVerified.gbk",
  25. L1L2_st2 => "pL1L2_ST2_EUCOMM_MapVerified.gbk",
  26. Ty1_EGFP => "Ty1_EGFP_tag.gbk",
  27. V5_Flag_biotin => "V5_Flag_biotin_tag.gbk",
  28. 'L1L2_GOHANU' => 'pGOHANU_Promoter_L1L2_for_Sanger.gbk',
  29. 'L1L2_NTARU-0' => 'NorCOMM_L1L2_frame_0.gbk',
  30. 'L1L2_NTARU-1' => 'NorCOMM_L1L2_Frame_1.gbk',
  31. 'L1L2_NTARU-2' => 'NorCOMM_L1L2_frame_2.gbk',
  32. 'L1L2_NTARU-K' => 'NorCOMM_L1L2_Frame_K.gbk',
  33. );
  34.  
  35. my %filenames = reverse %cassette_vector_name_file_lookup;
  36.  
  37. print Dumper [
  38. map eval($_),
  39. uniq
  40. map Dumper($_),
  41. sort {
  42. $a->{key} cmp $b->{key}
  43. || $a->{qualifier} cmp $b->{qualifier}
  44. || $a->{sequence} cmp $b->{sequence}
  45. }
  46. map {
  47. my ( $k, $v ) = %$_;
  48. { key => $k, qualifier => $v->{tags}{note}, sequence => $v->{seq} }
  49. }
  50. get_all( "objects", grep exists $filenames{ basename($_) }, @ARGV )
  51. ];
  52.  
  53. exit 0;
Add Comment
Please, Sign In to add comment