Advertisement
Guest User

Untitled

a guest
Dec 2nd, 2016
81
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.21 KB | None | 0 0
  1. #!/software/bin/perl
  2.  
  3. use strict;
  4. use warnings;
  5.  
  6. use Bio::EnsEMBL::Registry;
  7.  
  8. use Getopt::Long;
  9.  
  10. my ($host, $port, $user, $pass, $species, $version, $registry);
  11.  
  12. GetOptions(
  13. 'host=s' => \$host,
  14. 'user=s' => \$user,
  15. 'pass=s' => \$pass,
  16. 'port=i' => \$port,
  17. 'species=s' => \$species,
  18. 'version=i' => \$version,
  19. 'registry=s' => \$registry,
  20. );
  21.  
  22.  
  23. my $reg = 'Bio::EnsEMBL::Registry';
  24.  
  25. if(defined($host) && $host =~ /staging|variation|livemirror|genebuild|ensdb/) {
  26. $port ||= 3306;
  27. $user ||= 'ensro';
  28. $pass ||= '';
  29. }
  30.  
  31. else {
  32. $host ||= 'ensembldb.ensembl.org';
  33. $port ||= 5306;
  34. $user ||= 'anonymous';
  35. $pass ||= '';
  36. }
  37.  
  38. $species ||= 'homo_sapiens';
  39.  
  40.  
  41. if(defined($registry)) {
  42. $reg->load_all($registry);
  43. }
  44. else {
  45. $reg->load_registry_from_db(
  46. -host => $host,
  47. -user => $user,
  48. -pass => $pass,
  49. -port => $port,
  50. -db_version => $version,
  51. );
  52. }
  53.  
  54. $DB::single = 1;
  55.  
  56. my $va = $reg->get_adaptor($species,'variation','variation');
  57. my $v = $va->fetch_by_name('rs699');
  58.  
  59. my $vf = shift @{$v->get_all_VariationFeatures};
  60.  
  61. $DB::single = 1;
  62.  
  63. my $sa = $reg->get_adaptor($species, 'core', 'slice');
  64. my $s = $sa->fetch_by_region('chromosome', 1, 230845294, 230846294);
  65.  
  66. $DB::single = 1;
  67.  
  68. 1;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement