[Bioperl-l] Fwd: How to run this "bp_search2gff.pl"
Jason Stajich
jason at bioperl.org
Mon Mar 16 14:35:34 UTC 2009
redirecting message again.
Begin forwarded message:
> From: Dan Bolser <dan.bolser at gmail.com>
> Date: March 16, 2009 6:05:19 AM PDT
> To: jason at bioperl.org
> Subject: Fwd: [Bioperl-l] How to run this "bp_search2gff.pl"
>
> Ajay?
>
>
> ---------- Forwarded message ----------
> From: Dan Bolser <dan.bolser at gmail.com>
> Date: 2009/3/16
> Subject: Re: [Bioperl-l] How to run this "bp_search2gff.pl"
> To: bioperl-l at lists.open-bio.org
>
>
> 2009/3/5 ajay <jason at bioperl.org>:
>> Respected sir,
>>
>>
>> good morning! sir my self ajay a resercher in the
>> bioinformatics ,sir please
>> guide me how to run this bioperl scripts
>>
>> sir i have compleated the local BLAST P serch of 255 fasta files
>> and the
>> blast result was saved in the (.txt) format sir i need to
>> convert the
>> blast result to GFF3 format using this script "bp_search2gff.pl"
>>
>>
>> but i can not understand how to input my BLAST result file and get
>> the
>> output file in the GFF3 format
>>
>> so please tell me the steps to input the blast result file
>>
>> what changes i have to made in the program
>
> If you run "bp_search2gff.pl -h" you get the following output, which
> should help answer your question:
>
> SYNOPSIS
> Usage:
> search2gff [-o outputfile] [-f reportformat] [-i
> inputfilename] OR
> file1 file2 ..
>
> DESCRIPTION
> This script will turn a protein Search report (BLASTP, FASTP,
> SSEARCH,
> AXT, WABA) into a GFF File.
>
> The options are:
>
> -i infilename - (optional) inputfilename, will read
> either ARGV files or from STDIN
> -o filename - the output filename [default STDOUT]
> -f format - search result format (blast,
> fasta,waba,axt)
> (ssearch is fasta format). default is
> blast.
> -t/--type seqtype - if you want to see query or hit
> information
> in the GFF report
> -s/--source - specify the source (will be algorithm
> name
> otherwise like BLASTN)
> --method - the method tag (primary_tag) of the
> features
> (default is similarity)
> --scorefunc - a string or a file that when parsed
> evaluates
> to a closure which will be passed a
> feature
> object and that returns the score to
> be printed
> --locfunc - a string or a file that when parsed
> evaluates
> to a closure which will be passed two
> features, query and hit, and returns the
> location (Bio::LocationI compliant)
> for the
> GFF3 feature created for each HSP; the
> closure
> may use the clone_loc() and create_loc()
> functions for convenience, see their
> PODs
> --onehsp - only print the first HSP feature for
> each hit
> -p/--parent - the parent to which HSP features
> should refer
> if not the name of the hit or query
> (depending
> on --type)
> --target/--notarget - whether to always add the Target tag
> or not
> -h - this help menu
> --version - GFF version to use (put a 3 here to
> use gff 3)
> --component - generate GFF component fields
> (chromosome)
> -m/--match - generate a ’match’ line which is a
> container
> of all the similarity HSPs
> --addid - add ID tag in the absence of --match
> -c/--cutoff - specify an evalue cutoff
>
> Additionally specify the filenames you want to process on the
> com-
> mand-line. If no files are specified then STDIN input is
> assumed. You
> specify this by doing: search2gff < file1 file2 file3
>
> AUTHOR
> Jason Stajich, jason-at-bioperl-dot-org
>
> Contributors
> Hilmar Lapp, hlapp-at-gmx-dot-net
>
> clone_loc
>
> Title : clone_loc
> Usage : my $l = clone_loc($feature->location);
> Function: Helper function to simplify the task of cloning
> locations
> for --locfunc closures.
>
> Presently simply implemented using
> Storable::dclone().
> Example :
> Returns : A L<Bio::LocationI> object of the same type and
> with the
> same properties as the argument, but physically
> different.
> All structured properties will be cloned as well.
> Args : A L<Bio::LocationI> compliant object
>
> create_loc
>
> Title : create_loc
> Usage : my $l = create_loc("10..12");
> Function: Helper function to simplify the task of creating
> locations
> for --locfunc closures. Creates a location from a
> feature-
> table formatted string.
>
> Example :
> Returns : A L<Bio::LocationI> object representing the
> location given
> as formatted string.
> Args : A GenBank feature-table formatted string.
>
> perl v5.8.8 2009-01-14
> BP_SEARCH2GFF(1)
>
>
> Dan.
>
>> sir plesae help me
>> i sahll be very thankfull to you
>>
>>
>> waiting for your kind responce
>>
>> sincerly your!s
>> Ajay
>>
>>
>>
>>
>>
>> Ajay Kumar Mahato
>> (Research Associate)
>> National Research Centre for Plant Biotechnology (IARI)
>> (http://www.nrcpb.org)
>>
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>
Jason Stajich
jason at bioperl.org
More information about the Bioperl-l
mailing list