summary_blast <- function(x){
uniq <- c( UniProt = uniqueN(x$uniprot), genes= uniqueN(x$gene), transcripts = uniqueN(x$transcript), proteins = uniqueN(x$protein) )
y <- x[,.( genes=uniqueN(gene), transcripts=uniqueN(transcript), proteins = uniqueN(protein) , total =.N ), by=.(uniprot, domain, genus, name) ]
y <- y[order(-y$genes, tolower(y$name)),]
annot <- colSums(y[, .(total, genes, transcripts, proteins) ])
z <- cbind(unique=uniq, annotations= annot)
attr(y, "counts") <- z
message(nrow(y), " rows")
y
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.