Redundant species in nr database
I used blastp to search for the homologs of the serine threonine kinase protein (present in mycobacterium tuberculosis). I used the nr database, and all other default settings. My blast results have multiple hits for the same species. For eg, it has many hits for mycobacterium abscessus. Each hit has a different accession no. but the species is the same.
This leads to all my visible hits being of a few number of species only, and I cannot view more distantly related species in the list.
Why does the nr database have multiple entries of the same protein in the same species? How can I solve my issue?
• 1,737 views
•
link
updated
by
Biostar
•
written
by
kritika.tauras •
•
written
by
kritika.tauras •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Confused about which database to use - nt or nr?
written by DNAngel
-
Diamond and Taxonomic information
written by emiliomastriani •
-
Efetch from eutilities gives '400 Bad Request'; how to fix this?
written by Shraddha •
-
how to run blastp against nr databse locally
written by Arora •
-
Why does some hits that show up in search against smaller database, but does not show up in entire …
written by sajal •
-
How can I download a fasta file of all homologs to a given protein?
written by A248 •
-
Strore in separated fasta files fasta sequences with same locus
written by Amy •
-
Command line blast database problem
written by jolespin
-
similarety search for homologous alignment
written by shoxan-rahman •
-
Will the real kinases please stand up?
written by cdsouthan
You have kind of answered your own question
Those must be multiple submissions (possibly for different strains/from different submitters) from the same species. Instead of doing a simple blast you may want to start looking at psi- or delta-blast to get more distant homologs.
The same problem occurs with both psi- an delta-blast. The result still has duplicate entries.
While it is going to be extremely painful you can un-select most of those mycobacterial entries going into round 2 of Psi- or Delta-blast.
Only other way to do this may be by using a local copy of
nrsuitably cleansed of these multiple accessions.You may want to try the blast at Ensembl bacteria.
It has an option to specifically select
Distant homologies(undersearch sensitivity) as an explicit option.Edit: Moved this to a comment since there is a limit of a maximum 25 genomes to search against. Perhaps that may work for you.