soedinglab / MMseqs2

MMseqs2: ultra fast and sensitive search and clustering suite
https://mmseqs.com
GNU General Public License v3.0
1.38k stars 195 forks source link

mmseqs search return no results #780

Open Leo-T-Zang opened 11 months ago

Leo-T-Zang commented 11 months ago

Expected Behavior

Expect search return identity value for query sequences

Current Behavior

No values (identity) has been returned.

Steps to Reproduce (for bugs)

Data can be found here: link

# Create MMseqs2 databases
mmseqs createdb gen.fasta gen_db
mmseqs createdb train.fasta train_db

# Perform the search
mmseqs search gen_db train_db result tmp -a -s 7.5
mmseqs convertalis gen_db train_db result result.m8

Use python code to read back results

mmseqs_result = pd.read_csv('./mmseq/result.m8', sep='\t', names=[
    'query', 'target', 'identity', 'alignment_length', 
    'mismatches', 'gap_openings', 'qstart', 'qend', 
    'tstart', 'tend', 'evalue', 'bit_score'
])

The dataframe is empty. Does this mean there is no identity (0 identity) for my query sequences in target set?

Leo-T-Zang commented 6 months ago

Hi MMseqs Team,

Can anyone help me with this problem?

Thanks a lot!