[Bioperl-l] Fwd: blast result not matching.

Peter Cock p.j.a.cock at googlemail.com
Mon Sep 5 06:25:15 EDT 2011

Please CC the mailing list.

On Mon, Sep 5, 2011 at 11:19 AM, Manju Rawat <manju.rawat2 at gmail.com> wrote:
> Hi,
> Thanks for the reply...
> but when i am blasting after getting sequence of any gene (from NCBI using
> bioperl see below)..it showing me same result as shown in NCBI..
> #!usr/bin/perl -w
> use Bio::Perl;
> $seq_object = get_sequence('NCBI',"NM_181451");
> $blast_result = blast_sequence($seq);
> write_blast(">roa1.blast",$blast_report);
> I dnt know why its not working when i am blasting my own sequence..

Maybe you need give the sequence as a FASTA entry
rather than a plain string?


More information about the Bioperl-l mailing list