Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/perl
- use strict;
- use warnings;
- use Bio::DB::Fasta;
- use Bio::SeqIO;
- my $db = Bio::DB::Fasta->new('dmel-all-chromosome-r6.13.fasta');
- my $out = Bio::SeqIO->new(-file => "out.fasta", -format => 'fasta');
- while(<DATA>) {
- chomp;
- my $seq = $db->get_Seq_by_id($_);
- # print $_, "\t", $seq->subseq(1,10) , "\n";
- print $out->write_seq($seq);
- }
- __DATA__
- 211000022278279
- 211000022278436
- 211000022278449
- 211000022278760
- 211000022279165
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement