Hi,
I used DIAMOND blastp to blast a series of proteins to the nr database and outputs a .daa files.
I then meganizae my .daa file and explore it with MEGAN. Works very smooth.
However, when I inspect a taxa and the proteins hits, there are a lots of matches that don’t have a taxonomy assigned (lines with the following format: “?; score=250.0”). I attach an example below.
It worries me as usually the matches with no taxonomy have a higher score that the match that has been retained (match having a taxonomy with the highest score). Seems like it would be that lots of sequences in the nr database don’t have a taxonomy assigned in the mapping file (I am using megan-map-Jan2021.db).
Is there a way to fix this issue? How do you guys make the mapping file?
Thank you very much
Best,
Thomas