New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
colnames and EC number #776
Comments
The columns are explained here: https://github.com/bbuchfink/diamond/wiki/1.-Tutorial
All options are explained in the Wiki.
These mappings can be downloaded at sites e.g. Uniprot. I'm to aware of a tool to do this though. |
@bbuchfink Thank you I didn't notice that. |
There are other tools that can do that, but diamond works only on proteins. |
@bbuchfink thank you |
@bbuchfink sorry again by in this command : diamond cluster -d INPUT_FILE -o OUTPUT_FILE --approx-id 30 -M 64G |
Supported formats are FASTA and DIAMOND (.dmnd), so you can provide after -d also your FASTA file. |
Yes of course you can use the blastx mode of diamond. |
hello,
I want to thank you for this amazing tool.
I tried to use it and it went really fast
I used these commands as you suggested :
downloading the tool
wget http://github.com/bbuchfink/diamond/releases/download/v2.1.8/diamond-linux64.tar.gz
tar xzf diamond-linux64.tar.gz
creating a diamond-formatted database file
./diamond makedb --in reference.fasta -d reference
running a search in blastp mode
./diamond blastp -d reference -q queries.fasta -o matches.tsv
running a search in blastx mode
./diamond blastx -d reference -q reads.fasta -o matches.tsv
Now I'm wondering what are the column names for my output data file
PNEG_00003T0 sp|Q09895|YAI8_SCHPO 45.7 387 195 7 1 373 1 386 7.08e-114 340
PNEG_00003T0 sp|Q9Y282|ERGI3_HUMAN 38.2 387 215 8 8 383 8 381 4.22e-83 261
second question : how to parameter the tool in term of e-value, qcov_hsp_perc etc
final question how I could determine EC number from this output file !!
thank you in advance !!
The text was updated successfully, but these errors were encountered: