muellan/metacache

Out file interpretation

bielasilva opened this issue · 2 comments

bielasilva commented

Hello,
I'm having trouble interpreting the per-read mapping report, specifically the taxid. I noticed there are some negative ones and when plotting with KronaTools most are not identified but some are classified as Eukaryota. Could you please clarify what they are?

Is there a recommended tool for plotting the results?

Thank you.

muellan commented

Hi,

since Metacache can classify down to the sequence level and individual sequences do not have an official NCBI taxid, we use negative taxids to identify the individual sequences in the current database. If you don't want/need that you can restrict the classification output with the option -lowest <rank name>, so e.g., -lowest species or -lowest subspecies.

bielasilva commented

Ah, got it! Thank you, I'll do that!