Advertisement
Guest User

Untitled

a guest
Sep 20th, 2017
50
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.44 KB | None | 0 0
  1. #!/usr/bin/perl
  2. use strict;
  3. use warnings;
  4.  
  5. use Bio::DB::Fasta;
  6. use Bio::SeqIO;
  7.  
  8. my $db = Bio::DB::Fasta->new('dmel-all-chromosome-r6.13.fasta');
  9.  
  10. my $out = Bio::SeqIO->new(-file => "out.fasta", -format => 'fasta');
  11.  
  12. while(<DATA>) {
  13. chomp;
  14. my $seq = $db->get_Seq_by_id($_);
  15. # print $_, "\t", $seq->subseq(1,10) , "\n";
  16. print $out->write_seq($seq);
  17. }
  18.  
  19.  
  20. __DATA__
  21. 211000022278279
  22. 211000022278436
  23. 211000022278449
  24. 211000022278760
  25. 211000022279165
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement