Skip to content
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

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

Open
Prangejet opened this issue Nov 18, 2024 · 0 comments

Comments

@Prangejet
Copy link

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant