Hi, It would be very useful to be able to export readName_to_topBlastHitPercentIdentity. i.e. for each read print the read name and the percent identity obtained from the top megan-considered blast hit. This would help to keep a lot of reads with poorer top blast hits, while also being aware that the taxonomic assignment for these hits is less trustworthy.
While this suggestion would be useful for reads assigned to species-level bins, it would not necessarily make sense for higher taxonomic levels. For these I would suggest that it would also be useful to be able to export the average percent identity of all megan-considered hits for a read, which would be more informative when a read is assigned to a bin at a higher taxonomic level (readName_to_averageBlastHitPercentIdentity).