hh-suite icon indicating copy to clipboard operation
hh-suite copied to clipboard

Error while creating an MSA for custom database creation

Open hites77 opened this issue 4 years ago • 2 comments

I ran the following code: First I ran the command to to build ffindex for my sequences kept in fasta file named as pgpb and then i ran command to generate MSA for each sequences:

mpirun -np 24 hhblits_mpi -i pgpb -d /scratch/users/nus/eriht/databases/uniclust30_2018_08 /uniclust30_2018_08 -oa3m pgpb_a3m_wo_ss -n 2 -cpu 24 -v 0

Got the following error

  • 17:46:54.836 ERROR: Error in /home/users/nus/eriht/scratch/hh-suite/src/hhalignment.cpp:1267: Compress:

  • 17:46:55.425 ERROR: no match states found in 1!

  • 17:46:55.434 ERROR: Better use '-M first' option or reduce gap percentage threshold for match states

Even put the -M 50 paramater but the error remain the same. Does for a sequence in my fasta file no match was found in uniprot database ?

hites77 avatar Feb 04 '20 07:02 hites77

It looks like one or your input MSAs has no column with more than 50% gaps. Is this possible? Does it also crash with -M first?

martin-steinegger avatar Feb 04 '20 18:02 martin-steinegger

MSA are to be build using this command. The same error occur even if I put -M first and even whne change the threshold to -M 10.

hites77 avatar Feb 05 '20 02:02 hites77