MMseqs2 icon indicating copy to clipboard operation
MMseqs2 copied to clipboard

extract sequences after mmseqs search

Open lzhangUT opened this issue 2 years ago • 0 comments

Hi Thanks for the great contributions here. I have a fasta file containing 36 sequences, I am trying to do the homology search against uniprot, Here is what I did: mmseqs databases UniProtKB db/uniprot tmp mmseqs createdb input_fasta queryDB mmseqs search queryDB db/uniprot resultDB tmp

it looks like everything is running fine, I got the resultDB.* files. Also I did mmseqs convertalis queryDB targetDB resultDB resultDB.m6 But my real goal here is to get the sequences for each fasta, not just the ids and scores, as my next step is to deduplicate the sequences and then do clustering. I did mmseqs convert2fasta resultsDB resultsDB.fasta

it seems like it is not converting: MMseqs Version: 7ebd2e0441e5c3bdec585317c2b1c3cdbf943568 Use header DB false Verbosity 3

Database resultsDB needs header information

Any help or suggestion is greatly appreciated. Thank you

lzhangUT avatar Aug 12 '22 02:08 lzhangUT