diamond icon indicating copy to clipboard operation
diamond copied to clipboard

Diamond blastp two sequences with themselves, there is 0 pairwise alignments reported

Open Prangejet opened this issue 1 year ago • 1 comments

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

Prangejet avatar Nov 18 '24 13:11 Prangejet

This is due to the repeat masking, you need to use --masking 0.

bbuchfink avatar Dec 23 '24 19:12 bbuchfink