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

Peter Cock p.j.a.cock at googlemail.com
Mon Sep 5 10:25:15 UTC 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?

Peter



More information about the Bioperl-l mailing list