diamond icon indicating copy to clipboard operation
diamond copied to clipboard

How to make diamond output species name and proteins?

Open wanglu2014 opened this issue 1 year ago • 3 comments

For example, now the output is WP_217954729. We want to add the information response to WP_217954729: Alistipes timonensis, sulfatase into output. https://www.ncbi.nlm.nih.gov/protein/WP_217954729.1/ I do know it is common for diamond to directly output species and protein names, however, it not work for my case. Moreover, I do not know how to deal with microbial protein on conversion of protein from refseq accession.

This is my code #NR database ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz #building database diamond makedb --in nr --db nr #less nr 图片 #Diamond diamond blastp --db nr.dmnd -q X2test.fasta -o proteinadd_matches_fmt6.txt --outfmt 6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore testfile.zip Thank you for your time.

wanglu2014 avatar Dec 23 '22 07:12 wanglu2014