SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Task: True homologue? Britta Bioinformatics 0 03-18-2014 06:44 AM
How to do a discontinuous MEGABLAST using megablast? maimaiti2008 Bioinformatics 4 08-28-2013 11:07 AM
Interesting Behavior: The megablast standalone behaves as blastn.. Who knows why??? pbrand Bioinformatics 0 05-04-2012 07:45 AM
help on a simple (?) assembly task Kennels Bioinformatics 3 08-08-2011 09:08 PM
is bedTools useful for this task? dnusol Bioinformatics 1 06-27-2011 07:59 AM

Reply
 
Thread Tools
Old 06-09-2014, 05:24 PM   #1
augenblick
Junior Member
 
Location: San Francisco Bay Area

Join Date: Jun 2014
Posts: 1
Default indexing the nt database with makembindex to improve blastn (task megablast) runtime

I am trying to dust mask the NT database and then to index it using makembindex. However, when I run megablast on the indexed and dust masked nt database the runtime is slower (about 3 times slower) than when I run megablast on nt without indexing. Could anyone please check the set of commands I am using? Perhaps there is an error in the parameters. I got these commands off the ncbi website's readme for megablast. Is there any way to troubleshoot why the runtime with indexing is so much slower?

Thanks!
________________________________
Masking command:
dustmasker -in nt -infmt blastdb -parse_seqids -outfmt maskinfo_asn1_bin -out nt_dust.asnb

Merging mask data into nt:
makeblastdb -in nt -input_type blastdb -dbtype nucl -parse_seqids -mask_data nt_dust.asnb -out nt_dustmasked -title "nt"

Indexing command:
makembindex -input nt_dustmasked -iformat blastdb -old_style_index false -db_mask 11

Megablast command:
blastn -num_threads 8 -evalue 1e-30 -perc_identity 90 -word_size 45 -task megablast -outfmt 0 -show_gis -query contigs.fa -db nt_dustmasked -use_index true -db_soft_mask 11 > contigs.out 2> contigs.err
augenblick is offline   Reply With Quote
Reply

Tags
blastn, indexing, masking, megablast, reference

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:30 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO