fmaumusINRA

Thanks a lot for your response. I tried reducing with "-k 5" which improved a bit and then got an error when trying "-k 4". But isn't it doing a...

My query is a reverse transcriptase (RT) domain. For instance >Athila41 LDAGVIYPISDSTWVSPVHCVPKKGGMTVVKNEKDELIPTRTITGHRMCIDYRKLNAASRKDHFPLPFIDQMLERLANHPYYCFLDGYSGFFQIPIHPNDQEKTTFTCPYGTFAYKRMPFGLCNAPATFQRCMTSIFSDLIEEMVEVFMDDFSVYGPSFSSCLLNLGRVLTRCEETNLVLNWEKCHFMVKEGIVLDHKISEKGIEVDKGKVEVMMQLQPPKTVKDIRSFLGHAGFYRRFIKDFA It is about 240 aa so I wonder if --min-length will help. My target is a plant genome....

OK, maybe I over simplified the question because I was actually counting merged hits with a minimum size of 520 nt (empirical, a bit less than RT size). I indeed...

Thanks a lot for your help, see you soon with another question, I am exploring still

Thanks for your reply, I have tried the following commands but got "Segmentation fault" right away whenn launching filterresult. Any idea what could cause this? mmseqs createdb FASTANUCLFILE FASTANUCLFILE.db mmseqs...

Thanks a lot, it's very kind of you to spare time trying to make this possible. Yea, nucleotides again! ;)

Hi, I was wondering if by chance you had any opportunity to progress on that. I am still interested! Thanks