Open SilasK opened 1 year ago
In my opinion, yes, but this is safer in theory, and so it is sorta a pick of convenience for scientific defensibility. There is no option to switch it out but, you can change that column name to any column that exits and that lambda to any lambda that will parse it into a list. You can also add subfamily EC from DBcan. That is probably the easiest way to get what you want.
Oh, and we only match to the dbCAN family, not the subfamily.
https://github.com/WrightonLabCSU/DRAM/blob/ac83ba7c6ed081c68abca6290be55a6ba7d28865/mag_annotator/summarize_genomes.py#L45