Hi i'm a begginer in Perl and my script don't return a correct value, please help me. Example:The id 16 should return me to the species " Homo sapiens chromosome 16 , GRCh38.p2 Primary Assembly " and nucleotide sequence.
#!/usr/bin/perl use strict; use warnings; use Bio::DB::GenBank; use Bio::SeqIO; my $output = Bio::SeqIO->new(-format=>"fasta",-file=>">output.fasta"); my $gb = Bio::DB::GenBank->new(); while(<STDIN>){ chomp; my $id = $_; my $seq = $gb->get_Seq_by_id($id); $output->write_seq($seq); }
In reply to BioPerl sequences by wirlleyd
| For: | Use: | ||
| & | & | ||
| < | < | ||
| > | > | ||
| [ | [ | ||
| ] | ] |