SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
commandline blast kaps General 2 02-17-2015 01:17 AM
Count BLAST hits and 'no hits' fibar Bioinformatics 1 12-17-2014 12:28 AM
MEGAN - how many BLAST hits to keep for each? camilaz Bioinformatics 3 07-08-2014 07:30 AM
BLAST filtering out 'no hits' in output kga1978 Bioinformatics 16 06-01-2013 04:38 AM
Local BLAST no hits found PSchneeb Bioinformatics 6 04-18-2012 05:16 AM

Reply
 
Thread Tools
Old 02-17-2015, 01:28 AM   #1
kaps
Member
 
Location: Uganda

Join Date: Jan 2015
Posts: 71
Default settings for the best hits in commandline blast

Hi,

How can I set for the best 10 hits when blasting with command line blast?
kaps
kaps is offline   Reply With Quote
Old 02-17-2015, 01:33 AM   #2
sarvidsson
Senior Member
 
Location: Berlin, Germany

Join Date: Jan 2015
Posts: 137
Default

Did you try
Code:
-max_target_seqs 10
?

"Best" doesn't mean the same thing to everyone - but max_target_seqs should select the hits with the highest e-value and bit score.
sarvidsson is offline   Reply With Quote
Old 02-17-2015, 02:10 AM   #3
kaps
Member
 
Location: Uganda

Join Date: Jan 2015
Posts: 71
Default

yes,
but its like -max_target_seqs 10 does not give 10 maximum hits since i am getting many hits, more than 10!!
kaps is offline   Reply With Quote
Old 02-17-2015, 02:48 AM   #4
sarvidsson
Senior Member
 
Location: Berlin, Germany

Join Date: Jan 2015
Posts: 137
Default

Could you post the results? Often you get several (partial) alignments against one database sequence, which then counts as one "target sequence".
sarvidsson is offline   Reply With Quote
Old 02-17-2015, 03:53 AM   #5
kaps
Member
 
Location: Uganda

Join Date: Jan 2015
Posts: 71
Default

Attached is a txt file with the results.
I wait for your comments!
Attached Files
File Type: txt lib4_17_2_15.txt (16.1 KB, 9 views)
kaps is offline   Reply With Quote
Old 02-17-2015, 03:59 AM   #6
sarvidsson
Senior Member
 
Location: Berlin, Germany

Join Date: Jan 2015
Posts: 137
Default

As far as I can see, there are only one query with more than ten result lines (query 4_S4_L001_R1_001_(paired)_trimmed_(paired)_contig_23); in that case, one database target sequence has two alignments. This is the expected behaviour.
sarvidsson is offline   Reply With Quote
Old 08-25-2016, 07:33 AM   #7
ssully
Member
 
Location: NYC

Join Date: Aug 2010
Posts: 43
Default

You are getting the 10 best target sequences for each query. For some of the queries, there is more than one hit (hsp) to the target. Those are returned too, so sometimes there are more than 10 lines per query. This is good behavior, e.g. if your query is cDNA and the target is a gene with an intron, the query is split into two matches.
ssully is offline   Reply With Quote
Reply

Tags
best hits, blast, commandline

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 12:15 AM.


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