diamond
diamond copied to clipboard
How to make diamond output species name and proteins?
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.