diamond
diamond copied to clipboard
Diamond blastp two sequences with themselves, there is 0 pairwise alignments reported
I used ./diamond makedb --in margin_seq.fasta -d margin_seq to build the Database.
Database report as follows:
diamond v2.1.9.163 (C) Max Planck Society for the Advancement of Science, Benjamin Buchfink, University of Tuebingen
Documentation, support and updates available at http://www.diamondsearch.org
Please cite: http://dx.doi.org/10.1038/s41592-021-01101-x Nature Methods (2021)
#CPU threads: 128
Scoring parameters: (Matrix=BLOSUM62 Lambda=0.267 K=0.041 Penalties=11/1)
Database input file: margin_seq.fasta
Opening the database file... [0s]
Loading sequences... [0s]
Masking sequences... [0.007s]
Writing sequences... [0s]
Hashing sequences... [0s]
Loading sequences... [0s]
Writing trailer... [0s]
Closing the input file... [0s]
Closing the database file... [0s]
Database sequences 2
Database letters 1547
Database hash b401acfc064057d03ec56f2e0cd4daa3
Total time 0.008000s
And then, ./diamond blastp -d margin_seq -q margin_seq.fasta -o margin_seq_diamond.csv -f 6 qseqid sseqid qlen slen pident length mismatch gapopen qstart qend sstart send evalue bitscore -p 108 --id 99 --query-cover 99 --ultra-sensitive -e 1e-5 was used to blastp. And report as follows:
Reported 0 pairwise alignments, 0 HSPs.
0 queries aligned.
The sequences file was supported as follows:. Github doesn't accept FASTA files and I supported the txt file instead. margin_seq.txt
This is due to the repeat masking, you need to use --masking 0.