Only One Blast Result Per Query (Lowest E-Value)
I'm BLASTING to get gene names, and am only interested in the lowest E-value BLAST result for each query. From there, I'm looking to extract the gene name. How can this be done? I'm using the online server version of BLAST, but can easily transition to Linux.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
E value > 1
written by mmj.g •
-
blastn difference in command line version vs webpage version
written by yinglin1202 •
-
the relationship between query sequence length and e-value/bit score in NCBI BLAST+
written by bill-zt •
-
FASTA to Gene Names (Is BLAST the best option)
written by gtasource
-
how to write a script to grep the lowest e-value for each gene?
written by Anny •
-
How to Determine the Most Conserved Sequence from BLAST Result?
written by Bahrani •
-
Extract fasta from BLAST hits
written by biotech
-
Transcriptome Analysis with only a fasta file
written by coreyhowe99 •
-
Blast E Value Calculation
written by darxsys
-
Retrieve Total Length Of Query Sequence From Text Blast Output Via Bioperl
written by Pukeko •
In command line, local Blast :
On Blast API this should be the
Algorithm parameters>General Parameters>Max matches in a query rangeQuestion. I recently read this article: https://academic.oup.com/bioinformatics/advance-article-abstract/doi/10.1093/bioinformatics/bty833/5106166?redirectedFrom=fulltext
Does max_target give the lowest e-value? Or just any result that matches reaches the threshold.
Hum interesting article, I would say, try the parameter on a well known dataset, if you get what you had expected, go ahead. If not, you can also filter the output yourself using bash command, python, perl etc...
have a look at these ongoing threads :
I couldn't reproduce the problem of max_target_seqs
Misunderstood parameter of NCBI BLAST
that one you thus most likely should not set. Just ask any number and then filter the output afterwards to select the one with the lowest e-value