[Bioperl-l] From Blast hits to Taxanomy lineage for Short DNA Sequences (reads)
shalabh.sharma7 at gmail.com
Wed Mar 9 12:34:57 EST 2011
The other way to deal with this that you can download
the gi_taxaid file from ncbi.
Convert all your GI's to taxaid and use Bio::DB:Taxanomy.
I think there are lot of other options too, if you are interested you can
search for the thread which i started long time back.
Hope this helps.
Scientific Computing Professional Associate (Bioinformatics Specialist)
Department of Marine Sciences
University of Georgia
Athens, GA 30602-3636
On Wed, Mar 9, 2011 at 4:20 AM, Miguel Pignatelli
<miguel.pignatelli at uv.es>wrote:
> Hi Abhishek,
> For a non bioperl related solution, take a look at Bio::LITE::Taxonomy.
> It has been design to deal with great number of sequences (it is fast and
> You may also find interesting the Blast2lca tool,
> It currently works with the best hits for each query (calculates the lower
> common ancestor), but if you want to use only the best hit, please drop me a
> Please, let me know if you need further help with any of these,
> On 08/03/11 22:42, Abhishek Pratap wrote:
>> Hi All
>> I have results from different megablast of short reads(DNA sequences)
>> and after extracting the tophit for each read I want to bin them by
>> their lineage creating a tree.
>> For example.
>> If blast query hits the reference ->
>> I want to get the lineage for this specie.
>> The final goal is to do the above mapping as efficiently as possible.
>> Any pointers will be appreciated.
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
More information about the Bioperl-l