I imported a BLAST file with taxa names on the header (and also tried helping that with a "synonyms" file). It seems that MEGAN is able to detect the taxa for each read, however, I find plenty of reads that only aligned to one specific species but MEGAN has them on a LCA (last common ancestor). For example:
- Have many reads displayed as Capreolus capreolus. Many hundred reads with 1 match only have been assigned to the LCA infraorder Pecora, many with high scores (some higher scores than on some reads assigned to species B).
Why are reads within Pecora not assigned to the specific species when the score is 119 with a 100% identity (e.g. Bos taurs)? Some reads assigned to Capreolus capreolus have lower scores and identities below 90%. Can I change some parameters to achieve better separation, or are the default ones for "naive LCA" accurate or possibly conservative enough?
Thanks in advance!